Search in sources :

Example 6 with VisPolygon

use of org.vcell.vis.vismesh.thrift.VisPolygon in project vcell by virtualcell.

the class MovingBoundaryMeshMapping method fromReader.

public VisMesh fromReader(MovingBoundaryReader reader, DomainType domainType, int timeIndex) {
    MeshInfo meshInfo = reader.getMeshInfo();
    ISize size = new ISize(meshInfo.xinfo.number(), meshInfo.yinfo.number(), 1);
    Plane plane = reader.getPlane(timeIndex);
    int numX = size.getX();
    int numY = size.getY();
    int dimension = 2;
    Vect3D origin = new Vect3D(meshInfo.xinfo.start, meshInfo.yinfo.start, 0.0);
    Vect3D extent = new Vect3D(meshInfo.xinfo.end - meshInfo.xinfo.start, meshInfo.yinfo.end - meshInfo.yinfo.start, 1.0);
    // invoke VisMesh() constructor
    VisMesh visMesh = new VisMesh(dimension, origin, extent);
    int currPointIndex = 0;
    HashMap<String, Integer> pointDict = new HashMap<String, Integer>();
    PointIndex pointIndex = reader.getPointIndex();
    if (domainType == DomainType.INSIDE || domainType == DomainType.OUTSIDE) {
        int volumeIndex = 0;
        for (int j = 0; j < numY; j++) {
            for (int i = 0; i < numX; i++) {
                Element mbElement = plane.get(i, j);
                if ((domainType == DomainType.INSIDE && (mbElement.position == Position.INSIDE || mbElement.position == Position.BOUNDARY)) || (domainType == DomainType.OUTSIDE && mbElement.position == Position.OUTSIDE)) {
                    ArrayList<Integer> elementVisIndices = new ArrayList<Integer>();
                    int[] elementPointIndices = mbElement.boundary();
                    // first and last index is same
                    int numUniqueElementIndices = elementPointIndices.length - 1;
                    for (int p = 0; p < numUniqueElementIndices; p++) {
                        Vect3Didx point = pointIndex.lookup(elementPointIndices[p]);
                        VisPoint visPoint = new VisPoint(point.x, point.y, point.z);
                        String visPointKey = toStringKey(visPoint);
                        Integer visPointIndex = pointDict.get(visPointKey);
                        if (visPointIndex == null) {
                            visPointIndex = currPointIndex++;
                            pointDict.put(visPointKey, visPointIndex);
                            visMesh.addToPoints(visPoint);
                        }
                        elementVisIndices.add(visPointIndex);
                    }
                    VisPolygon polygon = new VisPolygon(elementVisIndices);
                    polygon.setMovingBoundaryVolumeIndex(new MovingBoundaryVolumeIndex(volumeIndex));
                    visMesh.addToPolygons(polygon);
                    volumeIndex++;
                }
            // end if
            }
        // end i
        }
    // end j
    } else if (domainType == DomainType.MEMBRANE) {
        throw new RuntimeException("DomainType MEMBRANE not yet implemented for moving boundary vtk processing");
    }
    return visMesh;
}
Also used : Plane(cbit.vcell.solvers.mb.MovingBoundaryTypes.Plane) HashMap(java.util.HashMap) ISize(org.vcell.util.ISize) Element(cbit.vcell.solvers.mb.MovingBoundaryTypes.Element) ArrayList(java.util.ArrayList) VisPoint(org.vcell.vis.vismesh.thrift.VisPoint) Vect3D(org.vcell.vis.vismesh.thrift.Vect3D) VisMesh(org.vcell.vis.vismesh.thrift.VisMesh) PointIndex(cbit.vcell.solvers.mb.PointIndex) VisPolygon(org.vcell.vis.vismesh.thrift.VisPolygon) MeshInfo(cbit.vcell.solvers.mb.MovingBoundaryTypes.MeshInfo) Vect3Didx(cbit.vcell.solvers.mb.Vect3Didx) MovingBoundaryVolumeIndex(org.vcell.vis.vismesh.thrift.MovingBoundaryVolumeIndex) VisPoint(org.vcell.vis.vismesh.thrift.VisPoint)

Aggregations

VisPoint (org.vcell.vis.vismesh.thrift.VisPoint)6 VisPolygon (org.vcell.vis.vismesh.thrift.VisPolygon)6 HashMap (java.util.HashMap)4 ISize (org.vcell.util.ISize)4 VisMesh (org.vcell.vis.vismesh.thrift.VisMesh)4 Vect3D (org.vcell.vis.vismesh.thrift.Vect3D)3 ArrayList (java.util.ArrayList)2 Box3D (org.vcell.vis.core.Box3D)2 FiniteVolumeIndex (org.vcell.vis.vismesh.thrift.FiniteVolumeIndex)2 VisLine (org.vcell.vis.vismesh.thrift.VisLine)2 Element (cbit.vcell.solvers.mb.MovingBoundaryTypes.Element)1 MeshInfo (cbit.vcell.solvers.mb.MovingBoundaryTypes.MeshInfo)1 Plane (cbit.vcell.solvers.mb.MovingBoundaryTypes.Plane)1 PointIndex (cbit.vcell.solvers.mb.PointIndex)1 Vect3Didx (cbit.vcell.solvers.mb.Vect3Didx)1 ChomboBoundaries (org.vcell.vis.chombo.ChomboBoundaries)1 ChomboBox (org.vcell.vis.chombo.ChomboBox)1 ChomboLevel (org.vcell.vis.chombo.ChomboLevel)1 Covering (org.vcell.vis.chombo.ChomboLevel.Covering)1 ChomboLevelData (org.vcell.vis.chombo.ChomboLevelData)1