[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;

阅读更多
个人分类: G+smo
上一篇[G+smo] gsDofMapper 的使用
下一篇[G+smo]细分网格时控制点的排序
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭