From f13eacbc8e5ab714dd3544adab8189c313382c77 Mon Sep 17 00:00:00 2001
From: Ian Jauslin <ian.jauslin@roma1.infn.it>
Date: Wed, 22 Jul 2015 13:55:29 +0000
Subject: Support for non-commuting fields

---
 src/parse_file.c | 24 ++++++++++++++++++------
 1 file changed, 18 insertions(+), 6 deletions(-)

(limited to 'src/parse_file.c')

diff --git a/src/parse_file.c b/src/parse_file.c
index 6054372..19b91c6 100644
--- a/src/parse_file.c
+++ b/src/parse_file.c
@@ -44,16 +44,17 @@ limitations under the License.
 #define PP_EXTERNAL_MODE 8
 #define PP_INTERNAL_MODE 9
 #define PP_FERMIONS_MODE 10
+#define PP_NONCOMMUTING_MODE 11
 // indices
-#define PP_INDEX_MODE 11
+#define PP_INDEX_MODE 12
 // factors or monomials
-#define PP_BRACKET_MODE 12
+#define PP_BRACKET_MODE 13
 // labels
-#define PP_LABEL_MODE 13
+#define PP_LABEL_MODE 14
 // polynomial
-#define PP_POLYNOMIAL_MODE 14
+#define PP_POLYNOMIAL_MODE 15
 // group
-#define PP_GROUP_MODE 15
+#define PP_GROUP_MODE 16
 
 
 // parse fields list
@@ -101,6 +102,11 @@ int parse_input_fields(Char_Array str_fields, Fields_Table* fields){
 	  mode=PP_FERMIONS_MODE;
 	}
 	break;
+      case 'a':
+	if(mode==PP_NULL_MODE){
+	  mode=PP_NONCOMMUTING_MODE;
+	}
+	break;
 
       // reset buffer
       case ':':
@@ -123,6 +129,9 @@ int parse_input_fields(Char_Array str_fields, Fields_Table* fields){
 	else if(mode==PP_FERMIONS_MODE){
 	  int_array_append(i,&((*fields).fermions));
 	}
+	else if(mode==PP_NONCOMMUTING_MODE){
+	  int_array_append(i,&((*fields).noncommuting));
+	}
 	buffer_ptr=buffer;
 	*buffer_ptr='\0';
 	break;
@@ -142,6 +151,9 @@ int parse_input_fields(Char_Array str_fields, Fields_Table* fields){
 	else if(mode==PP_FERMIONS_MODE){
 	  int_array_append(i,&((*fields).fermions));
 	}
+	else if(mode==PP_NONCOMMUTING_MODE){
+	  int_array_append(i,&((*fields).noncommuting));
+	}
 	mode=PP_NULL_MODE;
 	break;
 
@@ -417,7 +429,7 @@ int parse_input_identities(Char_Array str_identities, Fields_Table* fields){
   (*fields).ids.length=0;
   for(i=0;i<tmp;i++){
     sign=1;
-    monomial_sort((*fields).ids.lhs[i], 0, (*fields).ids.lhs[i].length-1, *fields, &sign);
+    monomial_sort((*fields).ids.lhs[i], *fields, &sign);
     polynomial_simplify((*fields).ids.rhs+i, *fields);
     polynomial_multiply_Qscalar((*fields).ids.rhs[i],quot(sign,1));
   }
-- 
cgit v1.2.3-70-g09d2