diff --git a/src/serac/physics/tests/CMakeLists.txt b/src/serac/physics/tests/CMakeLists.txt index a1b24f86a..580ae53d6 100644 --- a/src/serac/physics/tests/CMakeLists.txt +++ b/src/serac/physics/tests/CMakeLists.txt @@ -44,6 +44,7 @@ set(physics_parallel_test_sources blt_list_append(TO physics_parallel_test_sources ELEMENTS contact_patch.cpp + contact_patch_tied.cpp contact_beam.cpp contact_solid_adjoint.cpp IF TRIBOL_FOUND) diff --git a/src/serac/physics/tests/contact_patch_tied.cpp b/src/serac/physics/tests/contact_patch_tied.cpp index 0eb767cf8..fbee17b8d 100644 --- a/src/serac/physics/tests/contact_patch_tied.cpp +++ b/src/serac/physics/tests/contact_patch_tied.cpp @@ -55,7 +55,7 @@ TEST_P(ContactPatchTied, patch) // }; // #elif defined(MFEM_USE_STRUMPACK) #ifdef MFEM_USE_STRUMPACK - LinearSolverOptions linear_options{.linear_solver = LinearSolver::Strumpack, .print_level = 1}; + LinearSolverOptions linear_options{.linear_solver = LinearSolver::Strumpack, .print_level = 0}; #else LinearSolverOptions linear_options{}; SLIC_INFO_ROOT("Contact requires MFEM built with strumpack."); @@ -64,8 +64,8 @@ TEST_P(ContactPatchTied, patch) NonlinearSolverOptions nonlinear_options{.nonlin_solver = NonlinearSolver::Newton, .relative_tol = 1.0e-10, - .absolute_tol = 1.0e-10, - .max_iterations = 20, + .absolute_tol = 5.0e-10, + .max_iterations = 25, .print_level = 1}; ContactOptions contact_options{.method = ContactMethod::SingleMortar, diff --git a/src/serac/physics/tests/contact_solid_adjoint.cpp b/src/serac/physics/tests/contact_solid_adjoint.cpp index 1dc3d9bab..61d6b1bbc 100644 --- a/src/serac/physics/tests/contact_solid_adjoint.cpp +++ b/src/serac/physics/tests/contact_solid_adjoint.cpp @@ -64,7 +64,10 @@ std::unique_ptr createContactSolver(const NonlinearSolverOptions& no auto solid = std::make_unique(nonlinear_opts, solid_mechanics::direct_linear_options, dyn_opts, physics_prefix + std::to_string(iter++), mesh_tag, std::vector{}, 0, 0.0, false, true); - solid->setMaterial(mat); + + + Domain whole_mesh = EntireDomain(StateManager::mesh(mesh_tag)); + solid->setMaterial(mat, whole_mesh); solid->setDisplacementBCs({2}, [](const mfem::Vector& /*X*/, double /*t*/, mfem::Vector& disp) { disp = 0.0; }); solid->setDisplacementBCs({4}, [](const mfem::Vector& /*X*/, double /*t*/, mfem::Vector& disp) {