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

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

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

revision 3005 by gross, Thu Apr 22 05:59:31 2010 UTC revision 3006 by gross, Thu Apr 22 08:31:28 2010 UTC
# Line 35  err_t Paso_FunctionDerivative(double* J0 Line 35  err_t Paso_FunctionDerivative(double* J0
35     norm_w=Paso_lsup(n,w,F->mpi_info);     norm_w=Paso_lsup(n,w,F->mpi_info);
36     tt=EPSILON*norm_x0;     tt=EPSILON*norm_x0;
37     ttt=sqrt(EPSILON)*norm_w;     ttt=sqrt(EPSILON)*norm_w;
38     #pragma omp parallel private(local_s, local_norm_w)     #pragma omp parallel private(local_s)
39     {     {
40        local_s=s;        local_s=s;
41        #pragma omp for private(i, aw)        #pragma omp for private(i, aw)
# Line 61  err_t Paso_FunctionDerivative(double* J0 Line 61  err_t Paso_FunctionDerivative(double* J0
61         norm_w=v[1];         norm_w=v[1];
62     }     }
63     #endif     #endif
64  printf("s ::  = %e, %e \n",s, norm_x0/norm_w);  /* printf("s ::  = %e, %e \n",s, norm_x0/norm_w); */
65     if (norm_w>0) {     if (norm_w>0) {
66          s=s*epsnew;          s=s*epsnew;
67  printf("s = %e\n",s);  /* printf("s = %e\n",s); */
68          Paso_LinearCombination(n,setoff,1.,x0,s,w);          Paso_LinearCombination(n,setoff,1.,x0,s,w);
69          err=Paso_FunctionCall(F,J0w,setoff,pp);          err=Paso_FunctionCall(F,J0w,setoff,pp);
70          if (err==SOLVER_NO_ERROR) {          if (err==SOLVER_NO_ERROR) {

Legend:
Removed from v.3005  
changed lines
  Added in v.3006

  ViewVC Help
Powered by ViewVC 1.1.26