26   DIFFPAR_MORE *neP;  in SNESDiffParameterCreate_More()  local
33   PetscCall(PetscNew(&neP));  in SNESDiffParameterCreate_More()
35   neP->function_count = 0;  in SNESDiffParameterCreate_More()
36   neP->fnoise_min     = 1.0e-20;  in SNESDiffParameterCreate_More()
37   neP->hopt_min       = 1.0e-8;  in SNESDiffParameterCreate_More()
38   neP->h_first_try    = 1.0e-3;  in SNESDiffParameterCreate_More()
39   neP->fnoise_resets  = 0;  in SNESDiffParameterCreate_More()
40   neP->hopt_resets    = 0;  in SNESDiffParameterCreate_More()
43   PetscCall(VecDuplicateVecs(x, 3, &neP->workv));  in SNESDiffParameterCreate_More()
44   w = neP->workv[0];  in SNESDiffParameterCreate_More()
54   if (flg) neP->fp = fopen(noise_file, "w");  in SNESDiffParameterCreate_More()
55   else neP->fp = fopen("noise.out", "w");  in SNESDiffParameterCreate_More()
56   PetscCheck(neP->fp, PETSC_COMM_SELF, PETSC_ERR_FILE_OPEN, "Cannot open file");  in SNESDiffParameterCreate_More()
59   *outneP = neP;  in SNESDiffParameterCreate_More()
65   DIFFPAR_MORE *neP = (DIFFPAR_MORE *)nePv;  in SNESDiffParameterDestroy_More()  local
70   PetscCall(VecDestroyVecs(3, &neP->workv));  in SNESDiffParameterDestroy_More()
71   err = fclose(neP->fp);  in SNESDiffParameterDestroy_More()
73   PetscCall(PetscFree(neP));  in SNESDiffParameterDestroy_More()
79   DIFFPAR_MORE *neP = (DIFFPAR_MORE *)nePv;  in SNESDiffParameterCompute_More()  local
97   w    = neP->workv[0];  in SNESDiffParameterCompute_More()
98   xp   = neP->workv[1];  in SNESDiffParameterCompute_More()
99   fvec = neP->workv[2];  in SNESDiffParameterCompute_More()
100   fp   = neP->fp;  in SNESDiffParameterCompute_More()
105   h        = neP->h_first_try;  in SNESDiffParameterCompute_More()
108   fcount   = neP->function_count;  in SNESDiffParameterCompute_More()
114     neP->h_first_try = h;  in SNESDiffParameterCompute_More()
122       neP->function_count++;  in SNESDiffParameterCompute_More()
177   if (*fnoise < neP->fnoise_min) {  in SNESDiffParameterCompute_More()
178 …p, "Resetting fnoise: fnoise1 = %g, fnoise_min = %g\n", (double)*fnoise, (double)neP->fnoise_min));  in SNESDiffParameterCompute_More()
179     *fnoise = neP->fnoise_min;  in SNESDiffParameterCompute_More()
180     neP->fnoise_resets++;  in SNESDiffParameterCompute_More()
182   if (*hopt < neP->hopt_min) {  in SNESDiffParameterCompute_More()
183 …tf(comm, fp, "Resetting hopt: hopt1 = %g, hopt_min = %g\n", (double)*hopt, (double)neP->hopt_min));  in SNESDiffParameterCompute_More()
184     *hopt = neP->hopt_min;  in SNESDiffParameterCompute_More()
185     neP->hopt_resets++;  in SNESDiffParameterCompute_More()
200   fcount = neP->function_count - fcount;  in SNESDiffParameterCompute_More()
201 …t_cum = %" PetscInt_FMT ", rerrf=%g, sqrt(noise)=%g, h_more=%g\n", fcount, neP->function_count, (d…  in SNESDiffParameterCompute_More()