Changeset 3ab391 in git

Timestamp:

Jan 16, 2009, 3:21:09 AM (15 years ago)

Author:

Motsak Oleksandr <motsak@…>

Branches:

(u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')

Children:

1378915d904a163ade5fd8be0a96eb3e73bac75c

Parents:

2653d3c2809f0b96d6da699e54066fdb07623b92

Message:

*motsak: jacob for modules


git-svn-id: file:///usr/local/Singular/svn/trunk@11323 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• Singular/iparith.cc (modified) (3 diffs)

Singular/iparith.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: iparith.cc,v 1.491 2009-01-15 17:22:02 Singular Exp $ */
+/* $Id: iparith.cc,v 1.492 2009-01-16 02:21:09 motsak Exp $ */
 
 /*
@@ -4045 +4045 @@
   return FALSE;
 }
+/*2
+ * compute Jacobi matrix of a module/matrix
+ * Jacobi(M) := ( diff(Mt,var(1))|, ... ,| diff(Mt,var(pVariables))  ),
+ * where Mt := transpose(M)
+ * Note that this is consistent with the current conventions for jacob in Singular,
+ * whereas M2 computes its transposed.
+ */
+static BOOLEAN jjJACOB_M(leftv res, leftv a)
+{
+
+        ideal id = (ideal)a->Data();
+        id = idTransp(id);
+        int W = IDELEMS(id);
+
+        ideal result = idInit(W * pVariables, id->rank); 
+        poly *p = result->m;
+        
+        for( int v = 1; v <= pVariables; v++ )
+  {
+    poly* q = id->m;
+          for( int i = 0; i < W; i++, p++, q++ ) 
+                  *p = pDiff( *q, v );
+  }
+        idDelete(&id);
+
+        res->data = (char *)result;
+  return FALSE;
+}
+
+
 static BOOLEAN jjKBASE(leftv res, leftv v)
 {
@@ -5084 +5114 @@
 ,{jjJACOB_P,    JACOB_CMD,       IDEAL_CMD,      POLY_CMD       ALLOW_PLURAL}
 ,{mpJacobi,     JACOB_CMD,       MATRIX_CMD,     IDEAL_CMD      ALLOW_PLURAL}
+,{jjJACOB_M,    JACOB_CMD,       MODUL_CMD,      MODUL_CMD      ALLOW_PLURAL}
 ,{jjJanetBasis, JANET_CMD,       IDEAL_CMD,      IDEAL_CMD      ALLOW_PLURAL}
 ,{jjKBASE,      KBASE_CMD,       IDEAL_CMD,      IDEAL_CMD      ALLOW_PLURAL}