anonymous_namespace{ParallelVectorMatrixAssembler.cpp} Namespace Reference

Functions
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_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data, std::vector< GlobalIndexType > &indices, ProcessLib::AbstractJacobianAssembler &jacobian_assembler, ProcessLib::Assembly::MultiMatrixElementCache &cache)
int	getNumberOfThreads ()

Function Documentation

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_M_data,
		std::vector< double > &	local_K_data,
		std::vector< double > &	local_b_data,
		std::vector< double > &	local_Jac_data,
		std::vector< GlobalIndexType > &	indices,
		ProcessLib::AbstractJacobianAssembler &	jacobian_assembler,
		ProcessLib::Assembly::MultiMatrixElementCache &	cache )

Definition at line 25 of file ParallelVectorMatrixAssembler.cpp.

{
    indices = NumLib::getIndices(mesh_item_id, dof_table);
 
    local_M_data.clear();
    local_K_data.clear();
    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_M_data,
        local_K_data, 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_M_data, local_K_data, local_b_data, local_Jac_data,
              indices);
}

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

getNumberOfThreads()

int anonymous_namespace{ParallelVectorMatrixAssembler.cpp}::getNumberOfThreads

(

)

Definition at line 61 of file ParallelVectorMatrixAssembler.cpp.

{
    char const* const num_threads_env = std::getenv("OGS_ASM_THREADS");
 
    if (!num_threads_env)
    {
        return 1;
    }
 
    if (std::strlen(num_threads_env) == 0)
    {
        OGS_FATAL("The environment variable OGS_ASM_THREADS is set but empty.");
    }
 
    std::string num_threads_str{num_threads_env};
    BaseLib::trim(num_threads_str);
 
    std::istringstream num_threads_iss{num_threads_str};
    int num_threads = -1;
 
    num_threads_iss >> num_threads;
 
    if (!num_threads_iss)
    {
        OGS_FATAL("Error parsing OGS_ASM_THREADS (= \"{}\").", num_threads_env);
    }
 
    if (!num_threads_iss.eof())
    {
        OGS_FATAL(
            "Error parsing OGS_ASM_THREADS (= \"{}\"): not read entirely, the "
            "remainder is \"{}\"",
            num_threads_env,
            num_threads_iss.str().substr(num_threads_iss.tellg()));
    }
 
    if (num_threads < 1)
    {
        OGS_FATAL(
            "You asked (via OGS_ASM_THREADS) to assemble with {} < 1 thread.",
            num_threads);
    }
 
    return num_threads;
}

References OGS_FATAL, and BaseLib::trim().