Function update_directional_field

tutorial/202_Sampling/main.cpp:29–65

void update_directional_field()
{
    using namespace Eigen;
    using namespace std;
    VectorXd rotationAngles;
    Eigen::VectorXi presSingIndices;
    presSingIndices=VectorXi::Zero(field.tb->cycles.rows());
    for (int i=0;i<singVertices.size();i++)
        presSingIndices(singVertices[i])=singIndices[i];
    
    double IPError;
    Eigen::VectorXi currIndices;
    directional::index_prescription(presSingIndices,N,globalRotation, field, rotationAngles, IPError);
    
    if (viewingMode==TRIVIAL_PRINCIPAL_MATCHING)
        directional::principal_matching(field);
    
    if (viewingMode==IMPLICIT_FIELD){
        bc.conservativeResize(b.rows(),3);
        for (int i=0;i<b.size();i++)
            bc.row(i)<<field.extField.block(b(i),0,1,3).normalized();
        
        directional::CartesianField powerField;
        directional::power_field(ftb, b, bc, Eigen::VectorXd::Constant(b.size(),-1), N,powerField);
        directional::power_to_raw(powerField, N, field,true);
        directional::principal_matching(field);
    }
    
    viewer.set_cartesian_field(field);
    
    if (viewingMode==TRIVIAL_ONE_SING)
        viewer.set_singularities(field, singVertices, singIndices);
    
    if ((viewingMode==TRIVIAL_PRINCIPAL_MATCHING)||(viewingMode==IMPLICIT_FIELD))
        viewer.set_singularities(field, field.singLocalCycles, field.singIndices);
    
}


void callbackFunc() {