/[escript]/trunk-mpi-branch/paso/src/Solver.c
ViewVC logotype

Diff of /trunk-mpi-branch/paso/src/Solver.c

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

revision 150 by jgs, Thu Sep 15 03:44:45 2005 UTC revision 153 by jgs, Tue Oct 25 01:51:20 2005 UTC
# Line 187  void Paso_Solver(Paso_SystemMatrix* A,do Line 187  void Paso_Solver(Paso_SystemMatrix* A,do
187                          finalizeIteration = FALSE;                          finalizeIteration = FALSE;
188                    } else if (errorCode==SOLVER_MAXITER_REACHED) {                    } else if (errorCode==SOLVER_MAXITER_REACHED) {
189                           Paso_setError(WARNING,"maximum number of iteration step reached.\nReturned solution does not fulfil stopping criterion.");                           Paso_setError(WARNING,"maximum number of iteration step reached.\nReturned solution does not fulfil stopping criterion.");
190                             if (options->verbose) printf("Maximum number of iterations reached.!\n");
191                    } else if (errorCode == SOLVER_INPUT_ERROR ) {                    } else if (errorCode == SOLVER_INPUT_ERROR ) {
192                           Paso_setError(SYSTEM_ERROR,"illegal dimension in iterative solver.");                           Paso_setError(SYSTEM_ERROR,"illegal dimension in iterative solver.");
193                             if (options->verbose) printf("Internal error!\n");
194                    } else if ( errorCode == SOLVER_BREAKDOWN ) {                    } else if ( errorCode == SOLVER_BREAKDOWN ) {
195                    if (cntIter <= 1) {                    if (cntIter <= 1) {
196                             Paso_setError(ZERO_DIVISION_ERROR, "fatal break down in iterative solver.");                             Paso_setError(ZERO_DIVISION_ERROR, "fatal break down in iterative solver.");
197                               if (options->verbose) printf("Uncurable break down!\n");
198                        } else {                        } else {
199                        if (options->verbose) printf("Breakdown at iter %d (residual = %e). Restarting ...\n", cntIter+totIter, tol);                        if (options->verbose) printf("Breakdown at iter %d (residual = %e). Restarting ...\n", cntIter+totIter, tol);
200                            finalizeIteration = FALSE;                            finalizeIteration = FALSE;
201                        }                        }
202                    } else if (errorCode == SOLVER_MEMORY_ERROR) {                    } else if (errorCode == SOLVER_MEMORY_ERROR) {
203                         Paso_setError(MEMORY_ERROR,"memory allocation failed.");                     Paso_setError(MEMORY_ERROR,"memory allocation failed.");
204                           if (options->verbose) printf("Memory allocation failed!\n");
205                    } else if (errorCode !=SOLVER_NO_ERROR ) {                    } else if (errorCode !=SOLVER_NO_ERROR ) {
206                         Paso_setError(SYSTEM_ERROR,"unidentified error in iterative solver.");                     Paso_setError(SYSTEM_ERROR,"unidentified error in iterative solver.");
207                           if (options->verbose) printf("Unidentified error!\n");
208                    }                    }
209               } else {               } else {
210                  if (options->verbose) printf(". convergence! \n");                  if (options->verbose) printf(". convergence! \n");
# Line 215  void Paso_Solver(Paso_SystemMatrix* A,do Line 220  void Paso_Solver(Paso_SystemMatrix* A,do
220          }          }
221     }     }
222  }  }
   
 /*  
  * $Log$  
  * Revision 1.2  2005/09/15 03:44:40  jgs  
  * Merge of development branch dev-02 back to main trunk on 2005-09-15  
  *  
  * Revision 1.1.2.2  2005/09/07 00:59:09  gross  
  * some inconsistent renaming fixed to make the linking work.  
  *  
  * Revision 1.1.2.1  2005/09/05 06:29:49  gross  
  * These files have been extracted from finley to define a stand alone libray for iterative  
  * linear solvers on the ALTIX. main entry through Paso_solve. this version compiles but  
  * has not been tested yet.  
  *  
  *  
  */  

Legend:
Removed from v.150  
changed lines
  Added in v.153

  ViewVC Help
Powered by ViewVC 1.1.26