Functions
void	assembleOneElement (const std::size_t mesh_item_id, ProcessLib::LocalAssemblerInterface &local_assembler, const NumLib::LocalToGlobalIndexMap &dof_table, const double t, const double dt, const GlobalVector &x, const GlobalVector &x_prev, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, ProcessLib::Assembly::MultiMatrixElementCache< 2 > &cache)
void	assembleForStaggeredSchemeOneElement (const std::size_t mesh_item_id, ProcessLib::LocalAssemblerInterface &local_assembler, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, const double t, const double dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, ProcessLib::Assembly::MultiMatrixElementCache< 2 > &cache)
void	assembleWithJacobianOneElement (const std::size_t mesh_item_id, ProcessLib::LocalAssemblerInterface &local_assembler, const NumLib::LocalToGlobalIndexMap &dof_table, const double t, const double dt, const GlobalVector &x, const GlobalVector &x_prev, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data, ProcessLib::AbstractJacobianAssembler &jacobian_assembler, ProcessLib::Assembly::MultiMatrixElementCache< 1 > &cache)
void	assembleWithJacobianForStaggeredSchemeOneElement (const std::size_t mesh_item_id, ProcessLib::LocalAssemblerInterface &local_assembler, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, const double t, const double dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data, ProcessLib::AbstractJacobianAssembler &jacobian_assembler, ProcessLib::Assembly::MultiMatrixElementCache< 1 > &cache)
void	runAssemblyImpl (BaseLib::PolymorphicRandomAccessContainerView< ProcessLib::LocalAssemblerInterface > const &local_assemblers, ThreadException &exception, auto local_matrix_output, auto assemble, std::ptrdiff_t const num_active_local_assemblers, auto get_element_id_callback)
	Runs the passed `assemble` functor for each active local assembler.
void	runAssembly (BaseLib::PolymorphicRandomAccessContainerView< ProcessLib::LocalAssemblerInterface > const &local_assemblers, std::vector< std::size_t > const *const active_elements, ThreadException &exception, auto local_matrix_output, auto assemble)

Function Documentation

◆ assembleForStaggeredSchemeOneElement()

void anonymous_namespace{ParallelVectorMatrixAssembler.cpp}::assembleForStaggeredSchemeOneElement	(	const std::size_t	mesh_item_id,
		ProcessLib::LocalAssemblerInterface &	local_assembler,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_tables,
		const double	t,
		const double	dt,
		std::vector< GlobalVector * > const &	x,
		std::vector< GlobalVector * > const &	x_prev,
		int const	process_id,
		std::vector< double > &	local_M_data,
		std::vector< double > &	local_K_data,
		std::vector< double > &	local_b_data,
		ProcessLib::Assembly::MultiMatrixElementCache< 2 > &	cache )

Definition at line 48 of file ParallelVectorMatrixAssembler.cpp.

{
    std::vector<std::vector<GlobalIndexType>> indices_of_processes;
    indices_of_processes.reserve(dof_tables.size());
    transform(cbegin(dof_tables), cend(dof_tables),
              back_inserter(indices_of_processes), [&](auto const* dof_table)
              { return NumLib::getIndices(mesh_item_id, *dof_table); });
 
    auto local_coupled_xs =
        ProcessLib::getCoupledLocalSolutions(x, indices_of_processes);
    auto const local_x = MathLib::toVector(local_coupled_xs);
 
    auto local_coupled_x_prevs =
        ProcessLib::getCoupledLocalSolutions(x_prev, indices_of_processes);
    auto const local_x_prev = MathLib::toVector(local_coupled_x_prevs);
 
    local_M_data.clear();
    local_K_data.clear();
    local_b_data.clear();
 
    local_assembler.assembleForStaggeredScheme(t, dt, local_x, local_x_prev,
                                               process_id, local_M_data,
                                               local_K_data, local_b_data);
 
    auto const& indices = indices_of_processes[process_id];
    cache.add(local_M_data, local_K_data, local_b_data, indices);
}

References ProcessLib::LocalAssemblerInterface::assembleForStaggeredScheme(), ProcessLib::getCoupledLocalSolutions(), NumLib::getIndices(), and MathLib::toVector().

◆ assembleOneElement()

void anonymous_namespace{ParallelVectorMatrixAssembler.cpp}::assembleOneElement	(	const std::size_t	mesh_item_id,
		ProcessLib::LocalAssemblerInterface &	local_assembler,
		const NumLib::LocalToGlobalIndexMap &	dof_table,
		const double	t,
		const double	dt,
		const GlobalVector &	x,
		const GlobalVector &	x_prev,
		std::vector< double > &	local_M_data,
		std::vector< double > &	local_K_data,
		std::vector< double > &	local_b_data,
		ProcessLib::Assembly::MultiMatrixElementCache< 2 > &	cache )

Definition at line 23 of file ParallelVectorMatrixAssembler.cpp.

{
    std::vector<GlobalIndexType> indices =
        NumLib::getIndices(mesh_item_id, dof_table);
 
    local_M_data.clear();
    local_K_data.clear();
    local_b_data.clear();
 
    auto const local_x = x.get(indices);
    auto const local_x_prev = x_prev.get(indices);
    local_assembler.assemble(t, dt, local_x, local_x_prev, local_M_data,
                             local_K_data, local_b_data);
 
    cache.add(local_M_data, local_K_data, local_b_data, std::move(indices));
}

References ProcessLib::LocalAssemblerInterface::assemble(), MathLib::EigenVector::get(), and NumLib::getIndices().

◆ assembleWithJacobianForStaggeredSchemeOneElement()

void anonymous_namespace{ParallelVectorMatrixAssembler.cpp}::assembleWithJacobianForStaggeredSchemeOneElement	(	const std::size_t	mesh_item_id,
		ProcessLib::LocalAssemblerInterface &	local_assembler,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_tables,
		const double	t,
		const double	dt,
		std::vector< GlobalVector * > const &	x,
		std::vector< GlobalVector * > const &	x_prev,
		int const	process_id,
		std::vector< double > &	local_b_data,
		std::vector< double > &	local_Jac_data,
		ProcessLib::AbstractJacobianAssembler &	jacobian_assembler,
		ProcessLib::Assembly::MultiMatrixElementCache< 1 > &	cache )

Definition at line 116 of file ParallelVectorMatrixAssembler.cpp.

{
    std::vector<std::vector<GlobalIndexType>> indices_of_processes;
    indices_of_processes.reserve(dof_tables.size());
    transform(cbegin(dof_tables), cend(dof_tables),
              back_inserter(indices_of_processes),
              [&](auto const* dof_table)
              { return NumLib::getIndices(mesh_item_id, *dof_table); });
 
    auto local_coupled_xs =
        ProcessLib::getCoupledLocalSolutions(x, indices_of_processes);
    auto const local_x = MathLib::toVector(local_coupled_xs);
 
    auto local_coupled_x_prevs =
        ProcessLib::getCoupledLocalSolutions(x_prev, indices_of_processes);
    auto const local_x_prev = MathLib::toVector(local_coupled_x_prevs);
 
    std::vector<GlobalIndexType> const& indices =
        indices_of_processes[process_id];
 
    local_b_data.clear();
    local_Jac_data.clear();
 
    jacobian_assembler.assembleWithJacobianForStaggeredScheme(
        local_assembler, t, dt, local_x, local_x_prev, process_id, local_b_data,
        local_Jac_data);
 
    if (local_Jac_data.empty())
    {
        OGS_FATAL(
            "No Jacobian has been assembled! This might be due to "
            "programming errors in the local assembler of the "
            "current process.");
    }
 
    cache.add(local_b_data, local_Jac_data, indices);
}

References ProcessLib::AbstractJacobianAssembler::assembleWithJacobianForStaggeredScheme(), ProcessLib::getCoupledLocalSolutions(), NumLib::getIndices(), OGS_FATAL, and MathLib::toVector().

◆ assembleWithJacobianOneElement()

void anonymous_namespace{ParallelVectorMatrixAssembler.cpp}::assembleWithJacobianOneElement	(	const std::size_t	mesh_item_id,
		ProcessLib::LocalAssemblerInterface &	local_assembler,
		const NumLib::LocalToGlobalIndexMap &	dof_table,
		const double	t,
		const double	dt,
		const GlobalVector &	x,
		const GlobalVector &	x_prev,
		std::vector< double > &	local_b_data,
		std::vector< double > &	local_Jac_data,
		ProcessLib::AbstractJacobianAssembler &	jacobian_assembler,
		ProcessLib::Assembly::MultiMatrixElementCache< 1 > &	cache )

Definition at line 84 of file ParallelVectorMatrixAssembler.cpp.

{
    std::vector<GlobalIndexType> const& indices =
        NumLib::getIndices(mesh_item_id, dof_table);
 
    local_b_data.clear();
    local_Jac_data.clear();
 
    auto const local_x = x.get(indices);
    auto const local_x_prev = x_prev.get(indices);
    jacobian_assembler.assembleWithJacobian(local_assembler, t, dt, local_x,
                                            local_x_prev, local_b_data,
                                            local_Jac_data);
 
    if (local_Jac_data.empty())
    {
        OGS_FATAL(
            "No Jacobian has been assembled! This might be due to "
            "programming errors in the local assembler of the "
            "current process.");
    }
 
    cache.add(local_b_data, local_Jac_data, indices);
}

References ProcessLib::AbstractJacobianAssembler::assembleWithJacobian(), MathLib::EigenVector::get(), NumLib::getIndices(), and OGS_FATAL.

◆ runAssembly()

void anonymous_namespace{ParallelVectorMatrixAssembler.cpp}::runAssembly	(	BaseLib::PolymorphicRandomAccessContainerView< ProcessLib::LocalAssemblerInterface > const &	local_assemblers,
		std::vector< std::size_t > const *const	active_elements,
		ThreadException &	exception,
		auto	local_matrix_output,
		auto	assemble )

Runs the passed assemble functor for each active local assembler.

Forwards all calls to runAssemblyImpl().

Definition at line 200 of file ParallelVectorMatrixAssembler.cpp.

{
    [[likely]] if (!active_elements)
    {
        std::ptrdiff_t const n_elements =
            static_cast<std::ptrdiff_t>(local_assemblers.size());
 
        runAssemblyImpl(local_assemblers, exception, local_matrix_output,
                        assemble, n_elements, [](auto i) { return i; });
    }
    else
    {
        std::ptrdiff_t const n_elements =
            static_cast<std::ptrdiff_t>(active_elements->size());
 
        runAssemblyImpl(local_assemblers, exception, local_matrix_output,
                        assemble, n_elements, [active_elements](auto i)
                        { return (*active_elements)[i]; });
    }
}

References runAssemblyImpl().

◆ runAssemblyImpl()

void anonymous_namespace{ParallelVectorMatrixAssembler.cpp}::runAssemblyImpl	(	BaseLib::PolymorphicRandomAccessContainerView< ProcessLib::LocalAssemblerInterface > const &	local_assemblers,
		ThreadException &	exception,
		auto	local_matrix_output,
		auto	assemble,
		std::ptrdiff_t const	num_active_local_assemblers,
		auto	get_element_id_callback )

Runs the passed assemble functor for each active local assembler.

Definition at line 163 of file ParallelVectorMatrixAssembler.cpp.

{
#pragma omp for nowait
    for (std::ptrdiff_t i = 0; i < num_active_local_assemblers; ++i)
    {
        [[unlikely]] if (exception)
        {
            continue;
        }
 
        auto const element_id = get_element_id_callback(i);
        auto& loc_asm = local_assemblers[element_id];
 
        try
        {
            assemble(element_id, loc_asm);
        }
        catch (...)
        {
            exception.capture();
            continue;
        }
 
        local_matrix_output(element_id);
    }
}

References ThreadException::capture().

Referenced by runAssembly().