370 |
memcpy(P_PRES_dot_AP,&loc_dots[1],sizeof(double)*order); |
memcpy(P_PRES_dot_AP,&loc_dots[1],sizeof(double)*order); |
371 |
#endif |
#endif |
372 |
R_PRES_dot_AP0=R_PRES_dot_P_PRES[0]; |
R_PRES_dot_AP0=R_PRES_dot_P_PRES[0]; |
|
|
|
373 |
/*** if sum_BREAKF is equal to zero a breakdown occurs. |
/*** if sum_BREAKF is equal to zero a breakdown occurs. |
374 |
*** iteration procedure can be continued but R_PRES is not the |
*** iteration procedure can be continued but R_PRES is not the |
375 |
*** residual of X_PRES approximation. |
*** residual of X_PRES approximation. |