gsC1SurfEdge.hpp

 
#pragma once
 
#include <gsUnstructuredSplines/src/gsG1AuxiliaryPatch.h>
 
namespace gismo
{
 
    template<short_t d,class T>
    void gsC1SurfEdge<d,T>::computeAuxTopology(){
        //gsMultiPatch<T> auxTop;
        for(size_t i = 0; i <  auxGeom.size(); i++){
            if(auxGeom[i].getPatch().orientation() == -1)
            {
                auxGeom[i].swapAxis();
//                gsInfo << "Changed axis on patch: " << auxGeom[i].getGlobalPatchIndex() << "\n";
            }
            //auxTop.addPatch(auxGeom[i].getPatch());
        }
        //auxTop.computeTopology();
        //return auxTop;
    }
 
 
    template<short_t d,class T>
    void gsC1SurfEdge<d,T>::reparametrizeInterface(const boundaryInterface & item){
        computeAuxTopology();
 
//        if(item.second().side().index() == 1 && item.first().side().index() == 3)
//            return repTop;
 
        // Right patch along the interface. Patch 0 -> v coordinate. Edge west along interface
        switch (item.first().side().index())
        {
            case 1:
                break;
            case 4: auxGeom[0].rotateParamClock();
                break;
            case 3: auxGeom[0].rotateParamAntiClock();
                break;
            case 2: auxGeom[0].rotateParamAntiClockTwice();
                break;
            default:
                break;
        }
 
        // Left patch along the interface. Patch 1 -> u coordinate. Edge south along interface
        switch (item.second().side().index())
        {
            case 3:
                break;
            case 4: auxGeom[1].rotateParamAntiClockTwice();
                break;
            case 2: auxGeom[1].rotateParamAntiClock();
                break;
            case 1: auxGeom[1].rotateParamClock();
                break;
            default:
                break;
        }
 
        computeAuxTopology();
    }
 
 
    template<short_t d,class T>
    void gsC1SurfEdge<d,T>::reparametrizeBoundary(const index_t bInd){
        computeAuxTopology();
 
        if(auxGeom[0].getOrient())
        {
            switch (bInd)
            {
                case 3:
                    break;
                case 2:
                    auxGeom[0].rotateParamClock();
                    break;
                case 4:
                    auxGeom[0].rotateParamAntiClockTwice();
                    break;
                case 1:
                    auxGeom[0].rotateParamAntiClock();
                    break;
            }
        }
        else {
            switch (bInd) {
                case 1:
                    break;
                case 4:
                    auxGeom[0].rotateParamClock();
                    break;
                case 2:
                    auxGeom[0].rotateParamAntiClockTwice();
                    break;
                case 3:
                    auxGeom[0].rotateParamAntiClock();
                    break;
            }
        }
 
        computeAuxTopology();
    }
 
} // namespace gismo