Changeset 23198
- Timestamp:
- 08/28/18 14:14:41 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/solutionsequences/solutionsequence_nonlinear.cpp
r23197 r23198 20 20 Vector<IssmDouble>* df = NULL; 21 21 Vector<IssmDouble>* ys = NULL; 22 23 Vec res;24 double resNorm;25 26 22 27 23 Loads* savedloads=NULL; … … 75 71 femmodel->profiler->Stop(SOLVER); 76 72 77 78 // VecDuplicate(uf->pvector->vector,&res);79 // VecCopy(uf->pvector->vector,res);80 // VecAssemblyBegin(res);81 // VecAssemblyEnd(res);82 // VecScale(pf->pvector->vector,-1.);83 // MatMultAdd(Kff->pmatrix->matrix,uf->pvector->vector,pf->pvector->vector,res);84 // VecScale(pf->pvector->vector,-1.);85 // VecView(pf->pvector->vector,PETSC_VIEWER_STDOUT_WORLD);86 // VecNorm(res,NORM_INFINITY,&resNorm);87 // PetscPrintf(PETSC_COMM_WORLD,"Count = %d, Res. norm: %g\n", count,resNorm);88 89 90 91 92 73 Mergesolutionfromftogx(&ug, uf,ys,femmodel->nodes,femmodel->parameters);delete ys; 93 74
Note:
See TracChangeset
for help on using the changeset viewer.