From 76d00b6b2b094023495123d18286455b3e259698 Mon Sep 17 00:00:00 2001
From: Jiri Pittner <jiri@pittnerovi.com>
Date: Thu, 27 Nov 2025 14:11:54 +0100
Subject: [PATCH] fourindex: by default ignore assignments to zero elements due
 to antisymmetry

---
 fourindex.h | 14 ++++++++++++--
 1 file changed, 12 insertions(+), 2 deletions(-)

diff --git a/fourindex.h b/fourindex.h
index 4001401..06cb7c0 100644
--- a/fourindex.h
+++ b/fourindex.h
@@ -1385,7 +1385,12 @@ if(i<j) elem = -elem;
 if(k<l) elem = -elem;
 int I = ASMat_index_1(i,j);
 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
 if (I>=NRSMat<T>::nn || J>=NRSMat<T>::nn) laerror("index out of range");
 if (!NRSMat<T>::v) laerror("access to unallocated fourindex_dense");
@@ -1401,7 +1406,12 @@ if(i<j) elem = -elem;
 if(k<l) elem = -elem;
 int I = ASMat_index_1(i,j);
 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
 if (I>=NRSMat<T>::nn || J>=NRSMat<T>::nn) laerror("index out of range");
 if (!NRSMat<T>::v) laerror("access to unallocated fourindex_dense");