TwoDMeshTopology.h

// -*- mode: c++; c-basic-offset:4 -*-
 
// This file is part of libdap, A C++ implementation of the OPeNDAP Data
// Access Protocol.
 
// Copyright (c) 2002,2003,2011,2012 OPeNDAP, Inc.
// Authors: Nathan Potter <ndp@opendap.org>
//          James Gallagher <jgallagher@opendap.org>
//          Scott Moe <smeest1@gmail.com>
//          Bill Howe <billhowe@cs.washington.edu>
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
//
// You can contact OPeNDAP, Inc. at PO Box 112, Saunderstown, RI. 02874-0112.
 
#ifndef _TwoDMeshTopology_h
#define _TwoDMeshTopology_h 1
 
#include <gridfields/type.h>
#include <gridfields/gridfield.h>
#include <gridfields/grid.h>
#include <gridfields/cellarray.h>
 
using namespace libdap;
 
namespace libdap {
class BaseType;
class Array;
class MeshDataVariable;
class Grid;
}
 
namespace ugrid {

 
class TwoDMeshTopology {
 
private:
    BaseType *d_meshVar;

    string dimension;

    vector<libdap::Array *> *nodeCoordinateArrays;
    string nodeDimensionName;
    int nodeCount;

    libdap::Array *faceNodeConnectivityArray;
    libdap::Array::Dim_iter fncNodesDim, fncFacesDim;
    string faceDimensionName;
    int faceCount;
 
    vector<MeshDataVariable *> *rangeDataArrays;

    //string edgeNodeConnectivityArrayName;
    //string faceEdgeConnectivityArrayName;
    //Array *faceEdgeConnectivityArray;
    //string faceFaceConnectivityArrayName;
    //Array *faceFaceConnectivityArray;
    // Not used. jhrg 4/15/15 vector<string> *faceCoordinateNames;
    vector<libdap::Array *> *faceCoordinateArrays;

    //vector<string> *edgeCoordinateNames;
    //vector<Array *> *edgeCoordinateArrays;
    GF::Grid *gridTopology;
    GF::GridField *d_inputGridField;
    GF::GridField *resultGridField;
 
    vector<int *> *sharedIntArrays;
    vector<float *> *sharedFloatArrays;
 
    vector<GF::Array *> gfArrays;
 
    GF::Node *fncCellArray;
 
    bool _initialized;
 
    void ingestFaceNodeConnectivityArray(libdap::BaseType *meshTopology, libdap::DDS *dds);
    void ingestNodeCoordinateArrays(libdap::BaseType *meshTopology, libdap::DDS *dds);
    void ingestFaceCoordinateArrays(libdap::BaseType *meshTopology, libdap::DDS *dds);
 
    GF::Node *getFncArrayAsGFCells(libdap::Array *fncVar);
    int getStartIndex(libdap::Array *array);
    GF::CellArray *getFaceNodeConnectivityCells();
 
    libdap::Array *getGFAttributeAsDapArray(libdap::Array *sourceArray, locationType rank,
        GF::GridField *resultGridField);
    libdap::Array *getGridFieldCellArrayAsDapArray(GF::GridField *resultGridField, libdap::Array *sourceFcnArray);
    // libdap::Array *getNewFncDapArray(libdap::Array *templateArray, int N);
 
    void setNodeCoordinateDimension(MeshDataVariable *mdv);
    void setFaceCoordinateDimension(MeshDataVariable *mdv);
 
public:
    TwoDMeshTopology();
    ~TwoDMeshTopology();
 
    void init(string meshVarName, libdap::DDS *dds);
 
    string meshVarName() const
    {
        return getMeshVariable()->name();
    }
 
    libdap::BaseType *getMeshVariable() const
    {
        return d_meshVar;
    }
 
    void buildBasicGfTopology();
    void applyRestrictOperator(locationType loc, string filterExpression);
 
    int getInputGridSize(locationType location);
    int getResultGridSize(locationType location);
 
    void convertResultGridFieldStructureToDapObjects(vector<libdap::BaseType *> *results);
 
    void setLocationCoordinateDimension(MeshDataVariable *mdv);
 
    void addIndexVariable(locationType location);
    void getResultIndex(locationType location, void *target);
 
    void getResultGFAttributeValues(string attrName, libdap::Type type, locationType rank, void *target);
};
 
} // namespace ugrid
 
#endif // _TwoDMeshTopology_h