[G+smo]输出控制点

103 阅读 0 评论 68 点赞

我是靠谱客的博主无情大雁，最近开发中收集的这篇文章主要介绍[G+smo]输出控制点，觉得挺不错的，现在分享给大家，希望可以做个参考。

概述

在一个gsAssemblerBase<T>的派生类的成员函数中：

注意：

for (unsigned int i=0; i<m_patches.nPatches(); ++i)

{
gsGeometry<> & geo = m_patches.patch(i);
gsDebug << "control points " << geo.coefs() << endl;

}

inline void evaluate(gsBasis<T> const & basis,
gsGeometryEvaluator<T> & geoEval,
gsMatrix<T> const      & quNodes)
{
            gsDebug << "gauss points " << quNodes << endl;
basis.active_into(quNodes.col(0), m_actives);

numActive = dof_per_node * m_actives.rows();

            gsMatrix<T> allCp = geoEval.geometry().coefs();
            control_point_coor.setZero(m_actives.rows(), allCp.cols());
            for (int i=0; i<m_actives.rows(); ++i)
            {
                control_point_coor.row(i) = allCp.row(m_actives(i, 0));
            }
            gsDebug << control_point_coor << endl;

//Calculate the 0,1,2 derivatives!
basis.evalAllDers_into(quNodes, 2, m_basisData);

geoEval.evaluateAt(quNodes);

localMat.setZero(numActive, numActive);
            localMat_localCoor = localMat;
}