| Rev | Line | |
|---|
| [12679] | 1 | Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/modules/Solverx/SolverxGsl.cpp
|
|---|
| 2 | ===================================================================
|
|---|
| 3 | --- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/modules/Solverx/SolverxGsl.cpp (revision 12577)
|
|---|
| 4 | +++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/modules/Solverx/SolverxGsl.cpp (revision 12578)
|
|---|
| 5 | @@ -42,7 +42,12 @@
|
|---|
| 6 | gsl_vector_view b;
|
|---|
| 7 | gsl_vector *x = NULL;
|
|---|
| 8 | gsl_permutation *p = NULL;
|
|---|
| 9 | -
|
|---|
| 10 | +#ifdef _HAVE_ADOLC_
|
|---|
| 11 | + // if we use Adol-C then the IssmDouble will be an adouble
|
|---|
| 12 | + // and the calls to gsl_... will not work
|
|---|
| 13 | + // and we should call a suitable wrapped solve instead
|
|---|
| 14 | + _error2_("SolverxGsl: should not be here with Adol-C");
|
|---|
| 15 | +#else
|
|---|
| 16 | /*A will be modified by LU decomposition. Use copy*/
|
|---|
| 17 | IssmDouble* Acopy = xNew<IssmDouble>(n*n);
|
|---|
| 18 | xMemCpy<IssmDouble>(Acopy,A,n*n);
|
|---|
| 19 | @@ -69,4 +74,5 @@
|
|---|
| 20 | gsl_permutation_free(p);
|
|---|
| 21 | gsl_vector_free(x);
|
|---|
| 22 | *pX=X;
|
|---|
| 23 | +#endif
|
|---|
| 24 | }/*}}}*/
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.