160 |
h[INDEX2(k,k-1,l)]=0; |
h[INDEX2(k,k-1,l)]=0; |
161 |
ApplyGivensRotations(2,&(g[k-1]),&(c[k-1]),&(s[k-1])); |
ApplyGivensRotations(2,&(g[k-1]),&(c[k-1]),&(s[k-1])); |
162 |
} |
} |
163 |
norm_of_residual=abs(g[k]); |
norm_of_residual=fabs(g[k]); |
164 |
maxIterFlag = (k>=iter_max); |
maxIterFlag = (k>=iter_max); |
165 |
convergeFlag = (abs(g[k]) <= abs_tol); |
convergeFlag = (fabs(g[k]) <= abs_tol); |
166 |
printf("FGMRES step %d: error %e (tol=%d)\n",k,abs(g[k]),abs_tol); |
printf("FGMRES step %d: error %e (tol=%d)\n",k,fabs(g[k]),abs_tol); |
167 |
} |
} |
168 |
} |
} |
169 |
/* |
/* |