---

vtkGDBStructReader.h

#ifndef __vtkGDBStructReader_h
#define __vtkGDBStructReader_h

#include "vtkPolyDataSource.h"

extern "C"
  {
  #include "gdb/struct.h"
  }

class VTK_EXPORT vtkGDBStructReader : public vtkPolyDataSource
{
public:
  static vtkGDBStructReader *New();
  vtkTypeMacro(vtkGDBStructReader,vtkPolyDataSource);
  void PrintSelf(ostream& os, vtkIndent indent);
  
  vtkSetMacro(GDBStruct,GDB*);
  vtkGetMacro(GDBStruct,GDB*);

  vtkPolyData* GetVerticesOutput();
  vtkPolyData* GetEdgesOutput();
  vtkPolyData* GetSurfacesOutput();
  vtkPolyData* GetSurfacesPickOutput();

  int* GetSurfaceCellIds(int surfaceId);

  int GetSurfaceId(int pickId);

  static void GetCellCenter(float center[3], vtkCell *cell);

  vtkSetMacro(MergePoints,int);
  vtkGetMacro(MergePoints,int);
  vtkBooleanMacro(MergePoints,int);

  vtkSetMacro(MergeSurfacePoints,int);
  vtkGetMacro(MergeSurfacePoints,int);
  vtkBooleanMacro(MergeSurfacePoints,int);

  B2SPLINE_POINTS GetGrid(int sid);

  int HaveGrid(int sid);

  int InsertNextGrid(int sid, B2SPLINE_POINTS grid);

  void SetGrid(int loc, B2SPLINE_POINTS grid);

  void CleanGridList();

  void UpdateInformation();
  void UpdateData(vtkDataObject *output);

 protected:
  vtkGDBStructReader();
  ~vtkGDBStructReader();
  vtkGDBStructReader(const vtkGDBStructReader&) {};
  void operator=(const vtkGDBStructReader&) {};

  void Execute();

  GDB *GDBStruct; // pointer to GDB structure to read from
  int MergePoints;
  int MergeSurfacePoints;

  B2SPLINE_POINTS *gridList; // must be allocated before use
  vtkIdList *gridId;
  int gridListSize;

  vtkIntArray *Ia;
  vtkIdList *Id;
  vtkTimeStamp ExecuteTime;
};

inline B2SPLINE_POINTS vtkGDBStructReader::GetGrid(int sid)
{
  int gid = this->gridId->IsId( sid );
  if (gid == -1)
    {
    vtkErrorMacro(<<"GetGrid(): sid not in the list");
    gid = 0;
    }
  return this->gridList[gid];
}

inline void vtkGDBStructReader::SetGrid(int loc, B2SPLINE_POINTS grid)
{ this->gridList[loc] = grid; }

inline int vtkGDBStructReader::HaveGrid(int sid)
{ return (this->gridId->IsId(sid) == -1) ? 0 : 1; }

inline int vtkGDBStructReader::GetSurfaceId(int pickId)
{ return this->Id->GetId(pickId); }

inline vtkPolyData* vtkGDBStructReader::GetVerticesOutput()
{ return this->GetOutput(0); }

inline vtkPolyData* vtkGDBStructReader::GetEdgesOutput()
{ return this->GetOutput(1); }

inline vtkPolyData* vtkGDBStructReader::GetSurfacesOutput()
{ return this->GetOutput(2); }

inline vtkPolyData* vtkGDBStructReader::GetSurfacesPickOutput()
{ return this->GetOutput(3); }

#endif  // __vtkGDBStructReader_h

---