surface.hpp serif/solver/preconditioner/recipes/surface.hpp serif/solver/preconditioner/computations/inverse.hpp serif/solver/preconditioner/backend/mfem/coupling.hpp serif/solver/preconditioner/backend/mfem/surface/fitted_h1.hpp src/include/serif/solver/preconditioner/graph.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/fields.hpp serif/discretization/blocks/system_layout.hpp serif/discretization/domain/concepts.hpp serif/discretization/domain/physical_domains.hpp serif/discretization/domain/types.hpp serif/discretization/functional_space/space.hpp serif/geometry/deformation/failure.hpp serif/geometry/deformation/models.hpp serif/geometry/deformation/reference_surface.hpp serif/geometry/deformation/spaces.hpp serif/geometry/deformation/surface_embedding.hpp serif/geometry/deformation/workspace.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/operators/backend/mfem/linearized_block.hpp serif/operators/linearized_system_view.hpp serif/solver/preconditioner/backend/mfem/coupling.hpp serif/solver/preconditioner/backend/mfem/inverse.hpp serif/solver/preconditioner/backend/mfem/inverse/amg.hpp serif/solver/preconditioner/backend/mfem/inverse/chebyshev.hpp serif/solver/preconditioner/backend/mfem/inverse/dense.hpp serif/solver/preconditioner/backend/mfem/inverse/diagonal.hpp serif/solver/preconditioner/backend/mfem/surface/calibrated_mass.hpp serif/solver/preconditioner/backend/mfem/surface/fitted_h1.hpp serif/solver/preconditioner/backend/mfem/surface/reference_mass.hpp serif/solver/preconditioner/computations/inverse.hpp serif/solver/preconditioner/computations/physical_blocks.hpp src/include/serif/solver/preconditioner/computations/surface.hpp serif/solver/preconditioner/contract.hpp serif/solver/preconditioner/failure.hpp serif/solver/preconditioner/keys.hpp serif/solver/preconditioner/recipes/block.hpp serif/solver/preconditioner/recipes/surface.hpp serif/solver/preconditioner/strategies/coupling/order.hpp serif/solver/preconditioner/strategies/inverse/adaptive_amg.hpp serif/solver/preconditioner/strategies/inverse/chebyshev.hpp serif/solver/preconditioner/strategies/inverse/dense.hpp serif/solver/preconditioner/strategies/inverse/diagonal.hpp serif/solver/preconditioner/strategies/inverse/diagonal_regularization.hpp serif/solver/preconditioner/strategies/inverse/fixed_amg.hpp serif/solver/preconditioner/strategies/surface/calibrated_mass.hpp serif/solver/preconditioner/strategies/surface/fitted_h1.hpp serif/solver/preconditioner/strategies/surface/reference_mass.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/computation.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/resource.hpp serif/utils/dependencies/snapshot.hpp serif/utils/error/codes.hpp serif/utils/error/errors.hpp serif/utils/misc/concepts/concat.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/std_helper/cleaning.hpp serif/utils/misc/std_helper/variadic_cleaning.hpp serif/utils/misc/terminal/colors.hpp serif/utils/misc/traits/type_uniqueness.hpp serif/utils/types/strings.hpp serif/utils/types/type_list.hpp algorithm array compare concepts cstddef cstdint cstdio cstdlib exception expected format functional iostream magic_enum/magic_enum_all.hpp memory mfem/linalg/densemat.hpp mfem/linalg/operator.hpp mfem/linalg/sparsemat.hpp nameof.hpp numeric optional span string string_view tuple type_traits unistd.h unordered_map utility vector version src/include/serif/solver/preconditioner/computations/surface.hpp src/include/serif/solver/preconditioner/graph.hpp src/include/serif/solver/preconditioner/preconditioner.hpp src/include/serif/solver/solver.hpp serif serif::solver serif::solver::preconditioner serif::solver::preconditioner::computations #pragmaonce #include"serif/solver/preconditioner/recipes/surface.hpp" #include"serif/solver/preconditioner/computations/inverse.hpp" #include"serif/solver/preconditioner/backend/mfem/coupling.hpp" #include"serif/solver/preconditioner/backend/mfem/surface/fitted_h1.hpp" namespaceserif::solver::preconditioner::computations { template<typenameField,typenameRow,typenameStrategy,typenameBackend> requires(!std::same_as<Strategy,strategies::surface::ReferenceSurfaceMass>) [[nodiscard]]autosurface_target(constSurfacePreconditioner<Field,Row,Strategy,Backend>&choice) { usingnamespaceutils::dependencies; usingRecipe=SurfacePreconditioner<Field,Row,Strategy,Backend>; usingTargets=backend::mfem::SurfaceTargets; returnmake_computation<keys::SurfaceTarget<typenameRecipe::Rows,typenameRecipe::Values>,InputList<keys::Block<typenameRecipe::Rows,typenameRecipe::Values>>>( [target=choice.surface.target](constoperators::backend::mfem::LinearizedBlock&physical)->std::expected<Targets,Failure> { if(target!=strategies::surface::CalibrationTarget::surface_jacobian) { returnstd::unexpected(Failure{.code=utils::errors::SERiFErrorCode::PRECONDITIONER_INCOMPATIBLE_STRATEGY,.message="AstandalonesurfacefitrequiresthesurfaceJacobiantarget.Coupleittoamaterialrecipetofitmaterialrelaxation."}); } returnTargets{.surface_jacobian=physical.action()}; }); } template<typenameField,typenameRow,typenameStrategy,typenameBackend,typenameMaterial> requires(!std::same_as<Strategy,strategies::surface::ReferenceSurfaceMass>) [[nodiscard]]autosurface_target(constSurfacePreconditioner<Field,Row,Strategy,Backend>&choice,constMaterial&) { usingnamespaceutils::dependencies; usingRecipe=SurfacePreconditioner<Field,Row,Strategy,Backend>; usingBlock=operators::backend::mfem::LinearizedBlock; usingTargets=backend::mfem::SurfaceTargets; returnmake_computation<keys::SurfaceTarget<typenameRecipe::Rows,typenameRecipe::Values>,InputList<keys::Block<typenameRecipe::Rows,typenameRecipe::Values>,typenameMaterial::Result,keys::Block<typename Material::Rows, typename Recipe::Values>,keys::Block<typename Recipe::Rows, typename Material::Values>>>( [target=choice.surface.target](constBlock&physical,constbackend::mfem::Inverse&material,constBlock&material_from_surface,constBlock&surface_from_material)->std::expected<Targets,Failure> { if(target==strategies::surface::CalibrationTarget::surface_jacobian) { returnTargets{.surface_jacobian=physical.action()}; } if(target!=strategies::surface::CalibrationTarget::material_schur) { returnstd::unexpected(Failure{.code=utils::errors::SERiFErrorCode::PRECONDITIONER_INVALID_OPTIONS,.message="ChoosethesurfaceJacobianormaterialSchurtargetforsurfacecalibration."}); } constautoschur=backend::mfem::try_schur_operator(physical.action(),material,{.first_from_second=material_from_surface.action(),.second_from_first=surface_from_material.action()}); if(!schur) { returnstd::unexpected(schur.error()); } returnTargets{.material_schur=*schur}; }); } template<typenameField,typenameRow,typenameStrategy,typenameBackend> requires(!std::same_as<Strategy,strategies::surface::FittedSurfaceH1>) [[nodiscard]]autosurface_mass(constSurfacePreconditioner<Field,Row,Strategy,Backend>&choice) { usingnamespaceutils::dependencies; usingRecipe=SurfacePreconditioner<Field,Row,Strategy,Backend>; usingResult=keys::SurfaceMassInverse<typenameRecipe::Rows,typenameRecipe::Values>; ifconstexpr(std::same_as<Strategy,strategies::surface::ReferenceSurfaceMass>) { returnmake_computation<Result,InputList<keys::SurfaceMetric<typenameRecipe::Values>>>( [options=choice.surface](conststd::vector<double>&metric) { returnbackend::mfem::try_reference_surface_mass_inverse(metric,options); }); } else { returnmake_computation<Result,InputList<keys::SurfaceMetric<typenameRecipe::Values>,keys::SurfaceTarget<typenameRecipe::Rows,typenameRecipe::Values>>>( [options=choice.surface](conststd::vector<double>&metric,constbackend::mfem::SurfaceTargets&targets) { returnbackend::mfem::try_calibrated_surface_mass_inverse(metric,targets,options); }); } } template<typenameField,typenameRow,typenameBackend> [[nodiscard]]autosurface_h1_fit(constSurfacePreconditioner<Field,Row,strategies::surface::FittedSurfaceH1,Backend>&choice) { usingnamespaceutils::dependencies; usingRecipe=SurfacePreconditioner<Field,Row,strategies::surface::FittedSurfaceH1,Backend>; returnmake_computation<keys::SurfaceH1Fit<typenameRecipe::Rows,typenameRecipe::Values>,InputList<keys::SurfaceMatrices<typenameRecipe::Values>,keys::SurfaceTarget<typenameRecipe::Rows,typenameRecipe::Values>>>( [options=choice.surface](conststd::shared_ptr<constbackend::mfem::SurfaceH1Matrices>&matrices,constbackend::mfem::SurfaceTargets&targets)->std::expected<backend::mfem::FittedSurfaceOperator,Failure> { if(!matrices) { returnstd::unexpected(Failure{.code=utils::errors::SERiFErrorCode::PRECONDITIONER_INVALID_OPERATOR,.message="SurfaceH1preparationrequiresacceptedreferencemassandstiffnessmatrices."}); } returnbackend::mfem::try_fitted_surface_h1(*matrices,targets,options); }); } template<typenameField,typenameRow,typenameBackend> [[nodiscard]]autosurface_positive_inverse(constSurfacePreconditioner<Field,Row,strategies::surface::FittedSurfaceH1,Backend>&choice) { usingnamespaceutils::dependencies; usingRecipe=SurfacePreconditioner<Field,Row,strategies::surface::FittedSurfaceH1,Backend>; returnmake_computation<keys::SurfacePositiveInverse<typenameRecipe::Rows,typenameRecipe::Values>,InputList<keys::SurfaceH1Fit<typenameRecipe::Rows,typenameRecipe::Values>>>( [strategy=choice.inverse](constbackend::mfem::FittedSurfaceOperator&fitted) { returninvert_matrix(fitted.positive_operator,strategy); }); } template<typenameField,typenameRow,typenameStrategy,typenameBackend> [[nodiscard]]autosurface_inverse(constSurfacePreconditioner<Field,Row,Strategy,Backend>&) { usingnamespaceutils::dependencies; usingRecipe=SurfacePreconditioner<Field,Row,Strategy,Backend>; ifconstexpr(std::same_as<Strategy,strategies::surface::FittedSurfaceH1>) { returnmake_computation<typenameRecipe::Result,InputList<keys::SurfaceH1Fit<typenameRecipe::Rows,typenameRecipe::Values>,keys::SurfacePositiveInverse<typenameRecipe::Rows,typenameRecipe::Values>>>( [](constbackend::mfem::FittedSurfaceOperator&fitted,constbackend::mfem::Inverse&positive) { returnbackend::mfem::try_signed_surface_inverse(positive,fitted.fit.sign); }); } else { returnmake_computation<typenameRecipe::Result,InputList<keys::SurfaceMassInverse<typenameRecipe::Rows,typenameRecipe::Values>>>( [](constbackend::mfem::SurfaceMassInverse&mass) { returnmass.inverse; }); } } }//namespaceserif::solver::preconditioner::computations