detail Namespace Reference

Classes
struct	EigenMatrixType
struct	EigenMatrixType< N, 1 >
struct	EigenMatrixType< 0, M >
struct	EigenMatrixType< 0, 1 >

Functions
void	rotateToLocal (const MeshLib::RotationMatrix &matR2local, std::vector< MathLib::Point3d > &points)
	rotate points to local coordinates More...
MeshLib::RotationMatrix	getRotationMatrixToGlobal (const unsigned element_dimension, const unsigned global_dim, const std::vector< MathLib::Point3d > &points)
template<typename Solutions , typename Vector >
void	applyKnownSolutions (std::vector< Solutions > const *const known_solutions, Vector &x)
	Applies known solutions to the solution vector x. More...

Function Documentation

applyKnownSolutions()

template<typename Solutions , typename Vector >

void detail::applyKnownSolutions	(	std::vector< Solutions > const *const	known_solutions,
		Vector &	x
	)

Applies known solutions to the solution vector x.

Definition at line 22 of file TimeDiscretizedODESystem.cpp.

{
    if (!known_solutions)
    {
        return;
    }
 
    for (auto const& bc : *known_solutions)
    {
        for (std::size_t i = 0; i < bc.ids.size(); ++i)
        {
            // TODO that might have bad performance for some Vector types, e.g.,
            // PETSc.
            MathLib::setVector(x, bc.ids[i], bc.values[i]);
        }
    }
    MathLib::LinAlg::finalizeAssembly(x);
}

References MathLib::LinAlg::finalizeAssembly(), and MathLib::setVector().

Referenced by NumLib::TimeDiscretizedODESystem< ODESystemTag::FirstOrderImplicitQuasilinear, NonlinearSolverTag::Newton >::applyKnownSolutions(), and NumLib::TimeDiscretizedODESystem< ODESystemTag::FirstOrderImplicitQuasilinear, NonlinearSolverTag::Picard >::applyKnownSolutions().

getRotationMatrixToGlobal()

MeshLib::RotationMatrix detail::getRotationMatrixToGlobal	(	const unsigned	element_dimension,
		const unsigned	global_dim,
		const std::vector< MathLib::Point3d > &	points
	)

get a rotation matrix to the global coordinates it computes R in x=R*x' where x is original coordinates and x' is local coordinates

Definition at line 34 of file ElementCoordinatesMappingLocal.cpp.

{
    Eigen::Matrix3d matR;
    // compute R in x=R*x' where x are original coordinates and x' are local
    // coordinates
    if (element_dimension == 1)
    {
        Eigen::Vector3d const xx =
            (Eigen::Map<Eigen::Vector3d const>(points[1].getCoords()) -
             Eigen::Map<Eigen::Vector3d const>(points[0].getCoords()))
                .normalized();
        if (global_dim == 2)
        {
            matR = GeoLib::compute2DRotationMatrixToX(xx);
        }
        else
        {
            matR = GeoLib::compute3DRotationMatrixToX(xx);
        }
        matR.transposeInPlace();
    }
    else if (global_dim == 3 && element_dimension == 2)
    {
        // get plane normal
        auto const [plane_normal, d] = GeoLib::getNewellPlane(points);
        // compute a rotation matrix to XY
        matR = GeoLib::computeRotationMatrixToXY(plane_normal);
        // set a transposed matrix
        matR.transposeInPlace();
    }
    return matR;
}

References GeoLib::compute2DRotationMatrixToX(), GeoLib::compute3DRotationMatrixToX(), GeoLib::computeRotationMatrixToXY(), and GeoLib::getNewellPlane().

Referenced by MeshLib::ElementCoordinatesMappingLocal::ElementCoordinatesMappingLocal().

rotateToLocal()

void detail::rotateToLocal	(	const MeshLib::RotationMatrix &	matR2local,
		std::vector< MathLib::Point3d > &	points
	)

rotate points to local coordinates

Definition at line 24 of file ElementCoordinatesMappingLocal.cpp.

{
    std::transform(points.begin(), points.end(), points.begin(),
                   [&matR2local](auto const& pnt) { return matR2local * pnt; });
}

Referenced by MeshLib::ElementCoordinatesMappingLocal::ElementCoordinatesMappingLocal().