ComputeSparsityPattern.cpp File Reference

Detailed Description

Copyright

Copyright (c) 2012-2024, OpenGeoSys Community (http://www.opengeosys.org) Distributed under a Modified BSD License. See accompanying file LICENSE.txt or http://www.opengeosys.org/project/license

Definition in file ComputeSparsityPattern.cpp.

#include "ComputeSparsityPattern.h"
#include <numeric>
#include "LocalToGlobalIndexMap.h"
#include "MeshLib/NodeAdjacencyTable.h"
#include "MeshLib/NodePartitionedMesh.h"

Include dependency graph for ComputeSparsityPattern.cpp:

Go to the source code of this file.

Namespaces
namespace	NumLib

Functions
GlobalSparsityPattern	computeSparsityPatternPETSc (NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh)
GlobalSparsityPattern	NumLib::computeSparsityPattern (LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh)
	Computes a sparsity pattern for the given inputs.

Function Documentation

computeSparsityPatternPETSc()

GlobalSparsityPattern computeSparsityPatternPETSc	(	NumLib::LocalToGlobalIndexMap const &	dof_table,
		MeshLib::Mesh const &	mesh )

Definition at line 21 of file ComputeSparsityPattern.cpp.

{
    assert(dynamic_cast<MeshLib::NodePartitionedMesh const*>(&mesh));
    auto const& npmesh =
        *static_cast<MeshLib::NodePartitionedMesh const*>(&mesh);
 
    auto const max_nonzeroes = dof_table.getNumberOfGlobalComponents() *
                               npmesh.getMaximumNConnectedNodesToNode();
 
    // The sparsity pattern is misused here in the sense that it will only
    // contain a single value.
    return GlobalSparsityPattern(1, max_nonzeroes);
}

References NumLib::LocalToGlobalIndexMap::getNumberOfGlobalComponents().

Referenced by NumLib::computeSparsityPattern().