source: issm/oecreview/Archive/12321-12677/ISSM-12577-12578.diff

u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/modules/Solverx/SolverxGsl.cpp

@@ -42 +42 @@
         gsl_vector_view  b;
         gsl_vector      *x = NULL;
         gsl_permutation *p = NULL;
-
+#ifdef _HAVE_ADOLC_
+        // if we use Adol-C then the IssmDouble will be an adouble
+        // and the calls to gsl_... will not work
+        // and we should call  a suitable wrapped solve instead
+        _error2_("SolverxGsl: should not be here with Adol-C");
+#else
         /*A will be modified by LU decomposition. Use copy*/
         IssmDouble* Acopy = xNew<IssmDouble>(n*n);
         xMemCpy<IssmDouble>(Acopy,A,n*n);
@@ -69 +74 @@
         gsl_permutation_free(p);
         gsl_vector_free(x);
         *pX=X;
+#endif
 }/*}}}*/