Changeset 13185 for issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp

Timestamp:

08/30/12 06:59:13 (13 years ago)

Author:

utke

Message:

NEW : register externally differentiated function for ADOL-C ; introduce a generic parameter class to hold a structure of such registry items

File:

• issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp

@@ -45 +45 @@
 }/*}}}*/
 #ifdef _HAVE_ADOLC_
+
+int EDF_for_solverx(int n, IssmPDouble *x, int m, IssmPDouble *y) {
+    if(m*(m+1)!=n)_error_("Stiffness matrix should be square!");
+    SolverxSeq(&y,x, x+m*m, m);
+    return 0;
+}
+
 void SolverxSeq(IssmDouble** pX,IssmDouble* A,IssmDouble* B,int n){//{{{
         /* if we use Adol-C then the IssmDouble will be an adouble
@@ -61 +68 @@
 
         /*First, transfer from IssmDouble to double all our matrices and vectors: */
-        pdoubleA=xNew<double>(n*n);
-        pdoubleB=xNew<double>(n);
+        pdoubleA=xNew<IssmPDouble>(n*n);
+        pdoubleB=xNew<IssmPDouble>(n);
         for(i=0;i<n*n;i++)pdoubleA[i]=reCast<IssmPDouble>(A[i]);
         for(i=0;i<n;i++)pdoubleB[i]=reCast<IssmPDouble>(B[i]);