/[escript]/trunk/paso/src/MKL.c
ViewVC logotype

Diff of /trunk/paso/src/MKL.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 633 by dhawcroft, Thu Mar 23 05:37:00 2006 UTC revision 727 by gross, Fri May 12 06:31:06 2006 UTC
# Line 134  void Paso_MKL(Paso_SystemMatrix* A, Line 134  void Paso_MKL(Paso_SystemMatrix* A,
134                  &n, A->val, A->pattern->ptr, A->pattern->index, &idum, &nrhs,                  &n, A->val, A->pattern->ptr, A->pattern->index, &idum, &nrhs,
135                  iparm, &msglvl, in, out, &error);                  iparm, &msglvl, in, out, &error);
136             if (error != MKL_ERROR_NO) {             if (error != MKL_ERROR_NO) {
137               Paso_setError(ZERO_DIVISION_ERROR,"factorization in paradiso library failed.");               Paso_setError(ZERO_DIVISION_ERROR,"factorization in paradiso library failed. Most likely the matrix is singular.");
138               Paso_MKL_free(A);               Paso_MKL_free(A);
139             }             }
140             if (options->verbose) printf("timing MKL: LDU factorization: %.4e sec.\n",Paso_timer()-time0);             if (options->verbose) printf("timing MKL: LDU factorization: %.4e sec.\n",Paso_timer()-time0);
# Line 149  void Paso_MKL(Paso_SystemMatrix* A, Line 149  void Paso_MKL(Paso_SystemMatrix* A,
149                   iparm, &msglvl, in, out, &error);                   iparm, &msglvl, in, out, &error);
150          if (options->verbose) printf("timing MKL: solve: %.4e sec\n",Paso_timer()-time0);          if (options->verbose) printf("timing MKL: solve: %.4e sec\n",Paso_timer()-time0);
151          if (error != MKL_ERROR_NO) {          if (error != MKL_ERROR_NO) {
152                Paso_setError(VALUE_ERROR,"forward/backward substition in paradiso library failed.");                Paso_setError(VALUE_ERROR,"forward/backward substition in paradiso library failed. Most likely the matrix is singular.");
153          }          }
154       }       }
155       Performance_stopMonitor(pp,PERFORMANCE_ALL);       Performance_stopMonitor(pp,PERFORMANCE_ALL);

Legend:
Removed from v.633  
changed lines
  Added in v.727

  ViewVC Help
Powered by ViewVC 1.1.26