jiri/LA_library
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

t.cc example for profiling

Browse Source
This commit is contained in:
Jiri Pittner
2025-11-19 15:45:57 +01:00
parent c67d07ace0
commit e0f5485b94
1 changed files with 32 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

33
t.cc
View File

@@ -4346,7 +4346,7 @@ cout <<"xx = "<<xx.shape<< " "<<xx.names<<endl;
cout <<"Error = "<<(xx-xxx).norm()<<endl; cout <<"Error = "<<(xx-xxx).norm()<<endl;
} }
if(1) if(0)
{ {
int nn; int nn;
cin>>nn; cin>>nn;
@@ -4392,5 +4392,36 @@ cout <<"Error = "<<(z-zzz).norm()<<endl;
} }
if(1)
{
//for profiling
int nn;
cin>>nn;
int r=4;
NRVec<INDEXGROUP> s(r);
for(int i=0; i<r; ++i)
{
s[i].number=1;
s[i].symmetry=0;
s[i].range=nn;
s[i].offset=0;
}
Tensor<double> x(s); x.randomize(1.);
INDEXNAME xlist[] = {"i","j","k","l"};
x.names=NRVec<INDEXNAME>(xlist);
Tensor<double> y(s); y.randomize(1.);
INDEXNAME ylist[] = {"n","m","j","i"};
y.names=NRVec<INDEXNAME>(ylist);
Tensor<double>z(s); z.clear();
INDEXNAME zlist[] = {"k","l","m","n"};
z.names=NRVec<INDEXNAME>(zlist);
z.add_permuted_contractions("k/l,m",x,y,1,0,false,false);
cout <<z.norm()<<endl;
}
}//main }//main