/[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 1711 by artak, Tue Aug 19 03:37:25 2008 UTC revision 1778 by gross, Tue Sep 9 07:46:02 2008 UTC
# Line 83  err_t Paso_Solver_TFQMR( Line 83  err_t Paso_Solver_TFQMR(
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;
85    dim_t n = Paso_SystemMatrix_getTotalNumRows(A);    dim_t n = Paso_SystemMatrix_getTotalNumRows(A);
86    double  *u1=NULL, *u2=NULL, *y1=NULL, *y2=NULL, *d=NULL, *w=NULL, *v=NULL, *v_old=NULL, *tmp=NULL ;    double  *u1=NULL, *u2=NULL, *y1=NULL, *y2=NULL, *d=NULL, *w=NULL, *v=NULL, *v_old=NULL;
87    
88    double eta,theta,tau,rho,beta,alpha,sigma,rhon,c;    double eta,theta,tau,rho,beta,alpha,sigma,rhon,c;
89    
# Line 105  err_t Paso_Solver_TFQMR( Line 105  err_t Paso_Solver_TFQMR(
105    v=TMPMEMALLOC(n,double);    v=TMPMEMALLOC(n,double);
106    v_old=TMPMEMALLOC(n,double);    v_old=TMPMEMALLOC(n,double);
107        
     
   tmp=TMPMEMALLOC(n,double);  
     
108    
109   if (u1 ==NULL || u2== NULL || y1 == NULL || y2== NULL || d==NULL || w==NULL || v==NULL || v_old==NULL) {   if (u1 ==NULL || u2== NULL || y1 == NULL || y2== NULL || d==NULL || w==NULL || v==NULL || v_old==NULL) {
110       status=SOLVER_MEMORY_ERROR;       status=SOLVER_MEMORY_ERROR;

Legend:
Removed from v.1711  
changed lines
  Added in v.1778

  ViewVC Help
Powered by ViewVC 1.1.26