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

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

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

revision 1703 by artak, Thu Aug 14 05:34:25 2008 UTC revision 1706 by artak, Thu Aug 14 06:09:50 2008 UTC
# Line 75  err_t Paso_Solver_TFQMR( Line 75  err_t Paso_Solver_TFQMR(
75      Paso_Performance* pp) {      Paso_Performance* pp) {
76    
77    /* Local variables */    /* Local variables */
78      
79      int m=1;  
80      int j=0;
81      
82    dim_t num_iter=0,maxit;    dim_t num_iter=0,maxit;
83    bool_t breakFlag=FALSE, maxIterFlag=FALSE, convergeFlag=FALSE;    bool_t breakFlag=FALSE, maxIterFlag=FALSE, convergeFlag=FALSE;
84    err_t status = SOLVER_NO_ERROR;    err_t status = SOLVER_NO_ERROR;
# Line 84  err_t Paso_Solver_TFQMR( Line 88  err_t Paso_Solver_TFQMR(
88    double eta,theta,tau,rho,beta,alpha,sigma,rhon,c;    double eta,theta,tau,rho,beta,alpha,sigma,rhon,c;
89    
90    double norm_of_residual;    double norm_of_residual;
91      
92  /*                                                                 */  /*                                                                 */
93  /*-----------------------------------------------------------------*/  /*-----------------------------------------------------------------*/
94  /*                                                                 */  /*                                                                 */
# Line 150  err_t Paso_Solver_TFQMR( Line 154  err_t Paso_Solver_TFQMR(
154    tau = Paso_l2(n,r,A->mpi_info);    tau = Paso_l2(n,r,A->mpi_info);
155        
156    rho = tau * tau;    rho = tau * tau;
157    int m=1;          
   int j=0;  
       
158    norm_of_residual=tau*sqrt ( m + 1 );    norm_of_residual=tau*sqrt ( m + 1 );
159        
160    while (!(convergeFlag || maxIterFlag || breakFlag || (status !=SOLVER_NO_ERROR) ))    while (!(convergeFlag || maxIterFlag || breakFlag || (status !=SOLVER_NO_ERROR) ))

Legend:
Removed from v.1703  
changed lines
  Added in v.1706

  ViewVC Help
Powered by ViewVC 1.1.26