Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/ParFUM-Tops-Dev/ParFUM_TOPS.h

Go to the documentation of this file.

#ifndef __PARFUM_TOPS___H
#define __PARFUM_TOPS___H

#include <ParFUM.h>
#include <ParFUM_internals.h>


#include "ParFUM_TOPS_Types.h"
#include "ParFUM_TOPS_CUDA.h"


typedef struct{
    FEM_Mesh *mesh;
    void *mAtt;
    AllocTable2d<unsigned char> *ElemData_T;
    AllocTable2d<unsigned char> *GhostElemData_T;
    AllocTable2d<unsigned char> *NodeData_T;
    AllocTable2d<unsigned char> *GhostNodeData_T;
    AllocTable2d<int> *ElemConn_T;
    AllocTable2d<FP_TYPE_LOW> *coord_T;
    AllocTable2d<int> *node_id_T;
    AllocTable2d<int> *elem_id_T;
    AllocTable2d<int> *n2eConn_T;

    // table mapping global ids to local indices
    CkHashtableT<CkHashtableAdaptorT<int>, int>* nodeIDHash;
    CkHashtableT<CkHashtableAdaptorT<int>, int>* elemIDHash;

    unsigned node_attr_size;
    unsigned elem_attr_size;
    unsigned model_attr_size;

    unsigned num_local_elem;
    unsigned num_local_node;

#ifdef CUDA
    TopModelDevice device_model;
#endif

} TopModel;




class TopNodeItr{
public:
        int parfum_index;
        TopModel *model;
};

class TopElemItr{
public:
        int parfum_index;
        TopModel *model;
};

class TopNodeElemItr{
public:
        int current_index;
        int numAdjElem;
        
        TopModel *model;
        
        long node;
};


class TopFacetItr{
public:
        TopModel *model;
        TopElemItr *elemItr;
        int whichFacet; // Will go from 0 to 3
};


int tops_lib_FP_Type_Size();


TopModel* topModel_Create_Init(int _nelnode);

TopModel* topModel_Create_Driver(int elem_attr_sz, int node_attr_sz, int model_attr_sz, void* mAtt);

void topModel_Destroy(TopModel* m);

TopNode topModel_InsertNode(TopModel*, float x, float y, float z);
TopNode topModel_InsertNode(TopModel*, double x, double y, double z);

void topModel_SuggestInitialSize(TopModel* m, unsigned numNodes, unsigned numElements);


void topNode_SetId(TopModel*, TopNode, TopID id);

void topNode_SetAttrib(TopModel*, TopNode, void*);

TopElement topModel_InsertElem(TopModel*, TopElemType, TopNode*);

void topElement_SetId(TopModel*, TopElement, TopID id);

int topElement_GetId (TopModel* m, TopElement e); 

void topElement_SetAttrib(TopModel*, TopElement, void*);

TopNode topModel_GetNodeAtId(TopModel*,TopID);

void* topNode_GetAttrib(TopModel*, TopNode);

void* topElement_GetAttrib(TopModel*, TopElement);

TopNode topElement_GetNode(TopModel*,TopElement,int idx);

TopElement topModel_GetElemAtId(TopModel*,TopID);

int topNode_GetId(TopModel* m, TopNode n);

int topModel_GetNNodes(TopModel *model);

int topElement_GetNNodes(TopModel* model, TopElement elem);
bool topElement_IsCohesive(TopModel* m, TopElement e);

void topNode_GetPosition(TopModel*model, TopNode node,float*x,float*y,float*z);
void topNode_GetPosition(TopModel*model, TopNode node,double*x,double*y,double*z);

TopNodeItr*  topModel_CreateNodeItr(TopModel*);

void topNodeItr_Destroy(TopNodeItr*);

void topNodeItr_Begin(TopNodeItr*);

bool topNodeItr_IsValid(TopNodeItr*);

void topNodeItr_Next(TopNodeItr*);

TopNode topNodeItr_GetCurr(TopNodeItr*);


TopElemItr*  topModel_CreateElemItr(TopModel*);

void topElemItr_Destroy(TopElemItr*);

void topElemItr_Begin(TopElemItr*);

bool topElemItr_IsValid(TopElemItr*);

void topElemItr_Next(TopElemItr*);

TopElement topElemItr_GetCurr(TopElemItr*);

void topModel_TestIterators(TopModel*m);


void top_retrieve_elem_data(TopModel* m);
void top_retrieve_node_data(TopModel* m);
void top_put_node_data(TopModel* m);


//==============================================================
//   New functions that have been implemented

TopNodeElemItr* topModel_CreateNodeElemItr (TopModel* m, TopNode n);
bool topNodeElemItr_IsValid (TopNodeElemItr* neitr);
void topNodeElemItr_Next (TopNodeElemItr* neitr);
TopElement topNodeElemItr_GetCurr (TopNodeElemItr* neitr);
void topNodeElemItr_Destroy (TopNodeElemItr* neitr);


int topFacet_GetNNodes (TopModel* m, TopFacet f);
TopNode topFacet_GetNode (TopModel* m, TopFacet f, int i);
TopElement topFacet_GetElem (TopModel* m, TopFacet f, int i);

bool topElement_IsValid(TopModel* m, TopElement e);

bool topVertex_IsBoundary (TopModel* m, TopVertex v);

TopVertex topNode_GetVertex (TopModel* m, TopNode n);



TopFacetItr* topModel_CreateFacetItr (TopModel* m);
void topFacetItr_Begin(TopFacetItr* itr);
bool topFacetItr_IsValid(TopFacetItr* itr);
void topFacetItr_Next(TopFacetItr* itr);
TopFacet topFacetItr_GetCurr (TopFacetItr* itr);
void topFacetItr_Destroy (TopFacetItr* itr);


//==============================================================
//   New functions to be implemented for the new code to work

// TODO: setup a correct boundary condition after loading the mesh
// TODO: fixup everything to work with tet4s
// TODO: fix everything to work with multiple element types
// TODO: write tests

void topModel_InsertCohesiveAtFacet (TopModel* m, int ElemType, TopFacet f);


#endif

---