TwoDMeshTopology.h

Go to the documentation of this file.

// -*- 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 std;
using namespace libdap;

namespace libdap {
    class BaseType;
    class Array;
    class MeshDataVariable;
    class Grid;
}

namespace ugrid {

class TwoDMeshTopology
{

private:
        BaseType *d_myMeshVar;

        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;

        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(string meshVarName, DDS &dds);
        TwoDMeshTopology();
        ~TwoDMeshTopology();
        void init(string meshVarName, libdap::DDS *dds);
        // void addDataVariable(MeshDataVariable *mdt);
        string name(){ return getMeshVariable()->name();}
    libdap::BaseType *getMeshVariable(){ return d_myMeshVar; }

    void buildBasicGfTopology();
    void applyRestrictOperator(locationType loc, string filterExpression);
    // void buildRestrictedGfTopology(locationType loc, string filterExpression);
    // void buildGridFieldsTopology();

    int getInputGridSize(locationType location);
    int getResultGridSize(locationType location);
    // void convertResultGridFieldToDapObjects(vector<libdap::BaseType *> *results);
    void convertResultGridFieldStructureToDapObjects(vector<libdap::BaseType *> *results);
    // void convertResultRangeVarsToDapObjects(vector<BaseType *> *results);
    // void restrictRange(vector<libdap::BaseType *> *results);
    // libdap::Array *restrictDapArray(libdap::Array*, libdap::Array::Dim_iter locationCoordinateDim, locationType gridLocation);
    void setLocationCoordinateDimension(MeshDataVariable *mdv);

    void addIndexVariable(locationType location);
    void getResultIndex(locationType location, void *target);

    // void getResultGFAttributeValues(libdap::Array *templateArray, locationType rank, void *target);
    void getResultGFAttributeValues(string attrName, libdap::Type type, locationType rank, void *target);
};

} // namespace ugrid

#endif // _TwoDMeshTopology_h