NumLib::anonymous_namespace{DOFTableUtil.cpp} Namespace Reference

Functions
template<typename CalculateNorm >
double	norm (GlobalVector const &x, unsigned const global_component, LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, CalculateNorm calculate_norm)
double	norm1 (GlobalVector const &x, unsigned const global_component, LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh)
double	norm2 (GlobalVector const &x, unsigned const global_component, LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh)
double	normInfinity (GlobalVector const &x, unsigned const global_component, LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh)

Function Documentation

norm()

template<typename CalculateNorm >

double NumLib::anonymous_namespace{DOFTableUtil.cpp}::norm	(	GlobalVector const &	x,
		unsigned const	global_component,
		LocalToGlobalIndexMap const &	dof_table,
		MeshLib::Mesh const &	mesh,
		CalculateNorm	calculate_norm )

Definition at line 24 of file DOFTableUtil.cpp.

{
#ifdef USE_PETSC
    x.setLocalAccessibleVector();
#endif
 
    double res = 0.0;
    auto const& ms = dof_table.getMeshSubset(global_component);
 
    assert(ms.getMeshID() == mesh.getID());
    for (MeshLib::Node const* node : ms.getNodes())
    {
        auto const value = getNonGhostNodalValue(
            x, mesh, dof_table, node->getID(), global_component);
        res = calculate_norm(res, value);
    }
    return res;
}

References MeshLib::Mesh::getID(), NumLib::LocalToGlobalIndexMap::getMeshSubset(), and NumLib::getNonGhostNodalValue().

norm1()

double NumLib::anonymous_namespace{DOFTableUtil.cpp}::norm1	(	GlobalVector const &	x,
		unsigned const	global_component,
		LocalToGlobalIndexMap const &	dof_table,
		MeshLib::Mesh const &	mesh )

Definition at line 45 of file DOFTableUtil.cpp.

{
    double const res =
        norm(x, global_component, dof_table, mesh,
             [](double res, double value) { return res + std::abs(value); });
 
#ifdef USE_PETSC
    return BaseLib::MPI::allreduce(res, MPI_SUM, BaseLib::MPI::Mpi{});
#endif
    return res;
}

References BaseLib::MPI::allreduce(), and NumLib::norm().

norm2()

double NumLib::anonymous_namespace{DOFTableUtil.cpp}::norm2	(	GlobalVector const &	x,
		unsigned const	global_component,
		LocalToGlobalIndexMap const &	dof_table,
		MeshLib::Mesh const &	mesh )

Definition at line 58 of file DOFTableUtil.cpp.

{
    double const res =
        norm(x, global_component, dof_table, mesh,
             [](double res, double value) { return res + value * value; });
 
#ifdef USE_PETSC
    return std::sqrt(
        BaseLib::MPI::allreduce(res, MPI_SUM, BaseLib::MPI::Mpi{}));
#endif
    return std::sqrt(res);
}

References BaseLib::MPI::allreduce(), and NumLib::norm().

normInfinity()

double NumLib::anonymous_namespace{DOFTableUtil.cpp}::normInfinity	(	GlobalVector const &	x,
		unsigned const	global_component,
		LocalToGlobalIndexMap const &	dof_table,
		MeshLib::Mesh const &	mesh )

Definition at line 72 of file DOFTableUtil.cpp.

{
    double const res =
        norm(x, global_component, dof_table, mesh, [](double res, double value)
             { return std::max(res, std::abs(value)); });
 
#ifdef USE_PETSC
    return BaseLib::MPI::allreduce(res, MPI_MAX, BaseLib::MPI::Mpi{});
#endif
    return res;
}

References BaseLib::MPI::allreduce(), NumLib::norm(), and normInfinity().

Referenced by normInfinity().