|
@@ -19,7 +19,7 @@
|
|
|
template <typename MT, typename DerivedV, typename DerivedF, typename DerivedE, typename DerivedEMAP>
|
|
template <typename MT, typename DerivedV, typename DerivedF, typename DerivedE, typename DerivedEMAP>
|
|
|
void igl::crouzeix_raviart_massmatrix(
|
|
void igl::crouzeix_raviart_massmatrix(
|
|
|
const Eigen::MatrixBase<DerivedV> & V,
|
|
const Eigen::MatrixBase<DerivedV> & V,
|
|
|
- const Eigen::MatrixBase<DerivedF> & F,
|
|
|
|
|
|
|
+ const Eigen::MatrixBase<DerivedF> & F,
|
|
|
Eigen::SparseMatrix<MT> & M,
|
|
Eigen::SparseMatrix<MT> & M,
|
|
|
Eigen::PlainObjectBase<DerivedE> & E,
|
|
Eigen::PlainObjectBase<DerivedE> & E,
|
|
|
Eigen::PlainObjectBase<DerivedEMAP> & EMAP)
|
|
Eigen::PlainObjectBase<DerivedEMAP> & EMAP)
|
|
@@ -27,7 +27,7 @@ void igl::crouzeix_raviart_massmatrix(
|
|
|
// All occurrences of directed "facets"
|
|
// All occurrences of directed "facets"
|
|
|
Eigen::Matrix<typename DerivedF::Scalar, Eigen::Dynamic, Eigen::Dynamic> allE;
|
|
Eigen::Matrix<typename DerivedF::Scalar, Eigen::Dynamic, Eigen::Dynamic> allE;
|
|
|
oriented_facets(F,allE);
|
|
oriented_facets(F,allE);
|
|
|
- Eigen::VectorXi _1;
|
|
|
|
|
|
|
+ Eigen::Matrix<typename DerivedF::Scalar, Eigen::Dynamic, 1> _1;
|
|
|
unique_simplices(allE,E,_1,EMAP);
|
|
unique_simplices(allE,E,_1,EMAP);
|
|
|
return crouzeix_raviart_massmatrix(V,F,E,EMAP,M);
|
|
return crouzeix_raviart_massmatrix(V,F,E,EMAP,M);
|
|
|
}
|
|
}
|