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

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

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

revision 3119 by jfenwick, Thu Jan 28 02:03:15 2010 UTC revision 3120 by gross, Mon Aug 30 10:48:00 2010 UTC
# Line 23  Line 23 
23  Paso_Function * Paso_Function_LinearSystem_alloc(Paso_SystemMatrix* A, double* b, Paso_Options* options)  Paso_Function * Paso_Function_LinearSystem_alloc(Paso_SystemMatrix* A, double* b, Paso_Options* options)
24  {  {
25      Paso_Function * out=NULL;      Paso_Function * out=NULL;
26      Paso_Solver_setPreconditioner(A,options);      Paso_SystemMatrix_setPreconditioner(A,options);
27      if (! Paso_noError()) return NULL;      if (! Paso_noError()) return NULL;
28      out=MEMALLOC(1,Paso_Function);      out=MEMALLOC(1,Paso_Function);
29      if (! Paso_checkPtr(out)) {      if (! Paso_checkPtr(out)) {
# Line 60  err_t Paso_Function_LinearSystem_call(Pa Line 60  err_t Paso_Function_LinearSystem_call(Pa
60      Paso_SystemMatrix* A=(Paso_SystemMatrix*)(F->more);      Paso_SystemMatrix* A=(Paso_SystemMatrix*)(F->more);
61      Paso_Copy(F->n,F->tmp,F->b); /* tmp=b */      Paso_Copy(F->n,F->tmp,F->b); /* tmp=b */
62      Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(PASO_ONE, A, arg,-PASO_ONE, F->tmp); /* tmp=(A*arg-tmp) */      Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(PASO_ONE, A, arg,-PASO_ONE, F->tmp); /* tmp=(A*arg-tmp) */
63      Paso_Solver_solvePreconditioner(A,value,F->tmp);  /* value=P*tmp */      Paso_SystemMatrix_solvePreconditioner(A,value,F->tmp);  /* value=P*tmp */
64      return NO_ERROR;      return NO_ERROR;
65  }  }

Legend:
Removed from v.3119  
changed lines
  Added in v.3120

  ViewVC Help
Powered by ViewVC 1.1.26