diff --git a/ProcessLib/CalculateSurfaceFlux/CalculateSurfaceFlux.cpp b/ProcessLib/CalculateSurfaceFlux/CalculateSurfaceFlux.cpp index 51557f5d4ed3e3964324e46c5048414f74025c8b..2582fafdc23efd95c3e86eaaa804a400e0cd5c44 100644 --- a/ProcessLib/CalculateSurfaceFlux/CalculateSurfaceFlux.cpp +++ b/ProcessLib/CalculateSurfaceFlux/CalculateSurfaceFlux.cpp @@ -41,13 +41,18 @@ CalculateSurfaceFlux::CalculateSurfaceFlux(MeshLib::Mesh& boundary_mesh, new NumLib::LocalToGlobalIndexMap(std::move(all_mesh_subsets), NumLib::ComponentOrder::BY_LOCATION)); - boost::optional<MeshLib::PropertyVector<std::size_t> const&> - bulk_element_ids(boundary_mesh.getProperties() - .template getPropertyVector<std::size_t>( - "OriginalSubsurfaceElementIDs")); - boost::optional<MeshLib::PropertyVector<std::size_t> const&> bulk_face_ids( + auto const bulk_element_ids = boundary_mesh.getProperties().template getPropertyVector<std::size_t>( - "OriginalFaceIDs")); + "OriginalSubsurfaceElementIDs"); + if (!bulk_element_ids) + OGS_FATAL( + "OriginalSubsurfaceElementIDs boundary mesh property not found."); + auto const bulk_face_ids = + boundary_mesh.getProperties().template getPropertyVector<std::size_t>( + "OriginalFaceIDs"); + if (!bulk_face_ids) + OGS_FATAL("OriginalFaceIDs boundary mesh property not found."); + const unsigned integration_order = 2; ProcessLib::createLocalAssemblers<CalculateSurfaceFluxLocalAssembler>( boundary_mesh.getDimension() + 1, // or bulk_mesh.getDimension()? diff --git a/ProcessLib/CalculateSurfaceFlux/MapBulkElementPoint.cpp b/ProcessLib/CalculateSurfaceFlux/MapBulkElementPoint.cpp index 88ee57a121ce0509a596672d88b98beb4b719ecf..e179cddcec2b7f9292e450e3aa7ff69a26dc2c8d 100644 --- a/ProcessLib/CalculateSurfaceFlux/MapBulkElementPoint.cpp +++ b/ProcessLib/CalculateSurfaceFlux/MapBulkElementPoint.cpp @@ -13,7 +13,7 @@ namespace ProcessLib { -MathLib::Point3d getBulkElementPoint(MeshLib::Quad const& quad, +MathLib::Point3d getBulkElementPoint(MeshLib::Quad const& /*quad*/, std::size_t const face_id, MathLib::WeightedPoint1D const& wp) { @@ -31,7 +31,7 @@ MathLib::Point3d getBulkElementPoint(MeshLib::Quad const& quad, } } -MathLib::Point3d getBulkElementPoint(MeshLib::Hex const& hex, +MathLib::Point3d getBulkElementPoint(MeshLib::Hex const& /*hex*/, std::size_t const face_id, MathLib::WeightedPoint2D const& wp) { @@ -89,10 +89,10 @@ MathLib::Point3d getBulkElementPoint(MeshLib::Mesh const& mesh, } // TODO disable the 3d elements in the local assembler creator -MathLib::Point3d getBulkElementPoint(MeshLib::Mesh const& mesh, - std::size_t bulk_element_id, - std::size_t bulk_face_id, - MathLib::WeightedPoint3D const& wp) +MathLib::Point3d getBulkElementPoint(MeshLib::Mesh const& /*mesh*/, + std::size_t /*bulk_element_id*/, + std::size_t /*bulk_face_id*/, + MathLib::WeightedPoint3D const& /*wp*/) { return MathLib::ORIGIN; } diff --git a/ProcessLib/Process.h b/ProcessLib/Process.h index 7139941fb81778eb368534bfbecfc7583004fb33..bdb72e08cfbd0b9dfae2e9e7b146360da7ab02a1 100644 --- a/ProcessLib/Process.h +++ b/ProcessLib/Process.h @@ -102,7 +102,7 @@ public: return _secondary_variables; } - // used as call back for CalculateSurfaceFlux process + // Used as a call back for CalculateSurfaceFlux process. virtual std::vector<double> getFlux(std::size_t /*element_id*/, MathLib::Point3d const& /*p*/, GlobalVector const& /*x*/) const