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 23 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 44 of file DOFTableUtil.cpp.

{
    double res =
        norm(x, global_component, dof_table, mesh,
             [](double res, double value) { return res + std::abs(value); });
 
#ifdef USE_PETSC
    double global_result = 0.0;
    MPI_Allreduce(&res, &global_result, 1, MPI_DOUBLE, MPI_SUM,
                  PETSC_COMM_WORLD);
    res = global_result;
#endif
    return res;
}

References 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 60 of file DOFTableUtil.cpp.

{
    double res =
        norm(x, global_component, dof_table, mesh,
             [](double res, double value) { return res + value * value; });
 
#ifdef USE_PETSC
    double global_result = 0.0;
    MPI_Allreduce(&res, &global_result, 1, MPI_DOUBLE, MPI_SUM,
                  PETSC_COMM_WORLD);
    res = global_result;
#endif
    return std::sqrt(res);
}

References 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 76 of file DOFTableUtil.cpp.

{
    double res = norm(x, global_component, dof_table, mesh,
                      [](double res, double value)
                      { return std::max(res, std::abs(value)); });
 
#ifdef USE_PETSC
    double global_result = 0.0;
    MPI_Allreduce(&res, &global_result, 1, MPI_DOUBLE, MPI_MAX,
                  PETSC_COMM_WORLD);
    res = global_result;
#endif
    return res;
}

References NumLib::norm(), and normInfinity().

Referenced by normInfinity().