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

fourindex: by default ignore assignments to zero elements due to antisymmetry

Browse Source
This commit is contained in:
Jiri Pittner
2025-11-27 14:11:54 +01:00
parent 831404d08d
commit 76d00b6b2b
1 changed files with 12 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
fourindex.h
View File

@@ -1385,7 +1385,12 @@ if(i<j) elem = -elem;
if(k<l) elem = -elem; if(k<l) elem = -elem;
int I = ASMat_index_1(i,j); int I = ASMat_index_1(i,j);
int J = ASMat_index_1(k,l); int J = ASMat_index_1(k,l);
if (I<0 || J<0) laerror("assignment to nonexisting element"); if (I<0 || J<0)
#ifdef FORBID_NONEXISTENT_ASSIGNMENTS
laerror("assignment to nonexisting element");
#else
return;
#endif
#ifdef DEBUG #ifdef DEBUG
if (I>=NRSMat<T>::nn || J>=NRSMat<T>::nn) laerror("index out of range"); if (I>=NRSMat<T>::nn || J>=NRSMat<T>::nn) laerror("index out of range");
if (!NRSMat<T>::v) laerror("access to unallocated fourindex_dense"); if (!NRSMat<T>::v) laerror("access to unallocated fourindex_dense");
@@ -1401,7 +1406,12 @@ if(i<j) elem = -elem;
if(k<l) elem = -elem; if(k<l) elem = -elem;
int I = ASMat_index_1(i,j); int I = ASMat_index_1(i,j);
int J = ASMat_index_1(k,l); int J = ASMat_index_1(k,l);
if (I<0 || J<0) laerror("assignment to nonexisting element"); if (I<0 || J<0)
#ifdef FORBID_NONEXISTENT_ASSIGNMENTS
laerror("assignment to nonexisting element");
#else
return;
#endif
#ifdef DEBUG #ifdef DEBUG
if (I>=NRSMat<T>::nn || J>=NRSMat<T>::nn) laerror("index out of range"); if (I>=NRSMat<T>::nn || J>=NRSMat<T>::nn) laerror("index out of range");
if (!NRSMat<T>::v) laerror("access to unallocated fourindex_dense"); if (!NRSMat<T>::v) laerror("access to unallocated fourindex_dense");