diff --git a/PRIVATE b/PRIVATE index 5a769ec75..532cfadc1 160000 --- a/PRIVATE +++ b/PRIVATE @@ -1 +1 @@ -Subproject commit 5a769ec759d9dacc1866c35c6663cd0001e198c5 +Subproject commit 532cfadc1fda4f0f7f6a993e2bd029f14a50668f diff --git a/src/mesh/FEM_utilities.f90 b/src/mesh/FEM_utilities.f90 index a3856ccaa..d9e985ebd 100644 --- a/src/mesh/FEM_utilities.f90 +++ b/src/mesh/FEM_utilities.f90 @@ -123,9 +123,7 @@ subroutine FEM_utilities_init call PetscOptionsInsertString(PETSC_NULL_OPTIONS,'-mechanical_snes_type newtonls & &-mechanical_snes_linesearch_type cp -mechanical_snes_ksp_ew & &-mechanical_snes_ksp_ew_rtol0 0.01 -mechanical_snes_ksp_ew_rtolmax 0.01 & - &-mechanical_ksp_type fgmres -mechanical_ksp_max_it 25 & - &-mechanical_pc_type ml -mechanical_mg_levels_ksp_type chebyshev & - &-mechanical_mg_levels_pc_type sor -mechanical_pc_ml_nullspace user',ierr) + &-mechanical_ksp_type fgmres -mechanical_ksp_max_it 25', ierr) CHKERRQ(ierr) call PetscOptionsInsertString(PETSC_NULL_OPTIONS,num_mesh%get_asString('PETSc_options',defaultVal=''),ierr) CHKERRQ(ierr) @@ -151,6 +149,7 @@ subroutine utilities_constitutiveResponse(timeinc,P_av,forwardData) PetscErrorCode :: ierr + print'(/,a)', ' ... evaluating constitutive response ......................................' call homogenization_mechanical_response(timeinc,[1,mesh_maxNips],[1,mesh_NcpElems]) ! calculate P field