MAST::Numerics::libMeshWrapper Namespace Reference

Functions
template<typename P1 , int P2, typename P3 >
void	init_sparse_matrix (const libMesh::DofMap &dof_map, Eigen::SparseMatrix< P1, P2, P3 > &m)
template<typename IntType >
void	unconstrained_dofs (const libMesh::DofMap &dof_map, std::vector< IntType > &dofs)
	copies the unconstrained dofs for the dof_map in the vector dofs. More...

Function Documentation

init_sparse_matrix()

template<typename P1 , int P2, typename P3 >

void MAST::Numerics::libMeshWrapper::init_sparse_matrix	(	const libMesh::DofMap &	dof_map,
		Eigen::SparseMatrix< P1, P2, P3 > &	m
	)

Definition at line 36 of file sparse_matrix_initialization.hpp.

                                                           {

        Assert1(dof_map.comm().size() == 1,
                dof_map.comm().size(),
                "Eigen matrix can only be used for MPI communicator with rank 1.");
        
        m.resize(dof_map.n_dofs(), dof_map.n_dofs());
        
        m.reserve(dof_map.get_n_nz());
}

unconstrained_dofs()

template<typename IntType >

void MAST::Numerics::libMeshWrapper::unconstrained_dofs	(	const libMesh::DofMap &	dof_map,
		std::vector< IntType > &	dofs
	)

copies the unconstrained dofs for the dof_map in the vector dofs.

This is a collective operation and each rank will store only its own set of unconstrained dofs.

Definition at line 39 of file unconstrained_dofs.hpp.

                                                     {
    
    std::set<IntType>
    local_non_condensed_dofs_set;
    
    for (IntType i=dof_map.first_dof(); i<dof_map.end_dof(); i++) {
        
        if (!dof_map.is_constrained_dof(i))
            local_non_condensed_dofs_set.insert(i);
    }

    typename std::set<IntType>::const_iterator
    it   = local_non_condensed_dofs_set.begin(),
    end  = local_non_condensed_dofs_set.end();
    
    dofs.clear();
    dofs.reserve(local_non_condensed_dofs_set.size());
    
    for ( ; it != end; it++)
        dofs.push_back(*it);
}