298   Vec                         X, F, Y, G, W, GradF, YNtmp;  in SNESSolve_NEWTONTRDC()  local
322   YNtmp  = snes->work[4];  /* Newton solution */  in SNESSolve_NEWTONTRDC()
327   PetscCall(VecGetBlockSize(YNtmp, &bs));  in SNESSolve_NEWTONTRDC()
375     PetscCall(KSPSolve(snes->ksp, F, YNtmp)); /* Quasi Newton Solution */  in SNESSolve_NEWTONTRDC()
384       PetscCall(VecStrideNormAll(YNtmp, NORM_INFINITY, inorms));  in SNESSolve_NEWTONTRDC()
390         PetscCall(VecStrideScale(YNtmp, j, 1.0 / inorms[j]));  in SNESSolve_NEWTONTRDC()
405     PetscCall(VecNorm(YNtmp, NORM_2, &ynnorm)); /* ynnorm <- || Y_newton || */  in SNESSolve_NEWTONTRDC()
411         PetscCall(VecCopy(YNtmp, Y));  in SNESSolve_NEWTONTRDC()
429           PetscCall(VecAXPY(YNtmp, -1.0, YCtmp)); /* YCtmp = A, YNtmp = B */  in SNESSolve_NEWTONTRDC()
430           PetscCall(VecNorm(YNtmp, NORM_2, &c0)); /* this could be improved */  in SNESSolve_NEWTONTRDC()
432           PetscCall(VecDotRealPart(YCtmp, YNtmp, &c1));  in SNESSolve_NEWTONTRDC()
440           PetscCall(VecWAXPY(W, tau, YNtmp, YCtmp));  in SNESSolve_NEWTONTRDC()
447         PetscCall(VecScale(YNtmp, auk));  in SNESSolve_NEWTONTRDC()
448         PetscCall(VecCopy(YNtmp, Y)); /* this could be improved (many VecCopy, VecNorm)*/  in SNESSolve_NEWTONTRDC()