simplexD1linked.h

Go to the documentation of this file.

#ifndef SIMPLEXD1LINKED_H
#define SIMPLEXD1LINKED_H


#include <cassert>
#include <sstream>
using namespace std;

#include <typedefs.h>

template< typename INDX=uint>
class simplexD1linked
{
public:

  INDX pi[2];
  INDX ni[2];

  simplexD1linked() {};
  simplexD1linked
  (
    INDX const a, INDX const b, INDX const na=0, INDX const nb=0
  );

  void construct 
  (
    INDX const a, INDX const b, INDX const na=0, INDX const nb=0
  );

  void togglelinedirection()
  {
    INDX p0 = pi[0];
    INDX n0 = ni[0];
    pi[0] = pi[1];
    ni[0] = ni[1];
    pi[1] = p0;
    ni[1] = n0;
  }

  // Query Functions

  INDX const piInverse( INDX const gpt ) const
  {
    assert(gpt!=0);

    if (pi[0]==gpt) return 0;
    if (pi[1]==gpt) return 1;

    assert(false);
    return 3; // Return crap.
  }

  boolc isonboundary() const
  {
    if (ni[0]==0) return true;
    if (ni[1]==0) return true;
    return false;
  }
  
  void setnull()
    { pi[0] = pi[1] = ni[0] = ni[1] = 0; }

  boolc isnull() const
  {
    if (pi[0]!=0) return false;
    if (pi[1]!=0) return false;
    if (ni[0]!=0) return false;
    if (ni[1]!=0) return false;
    return true;
  }

  operator stringc () const;

  void serializeInverse(stringc & s)
  {
    stringstream w(s);
    w >> pi[0];
    w >> pi[1];
    w >> ni[0];
    w >> ni[1];
  }
};


//---------------------------------------------------------
// Implementation

template< typename INDX >
simplexD1linked<INDX>::operator stringc () const
{
  stringstream ss;
  ss << pi[0] << " " << pi[1] << " " << ni[0] << " " << ni[1];
  return ss.str();
}

template< typename INDX >
simplexD1linked<INDX>::simplexD1linked
(
  INDX const a, INDX const b, INDX const na, INDX const nb
)
{
  pi[0] = a;
  pi[1] = b;
  ni[0] = na;
  ni[1] = nb;
}

template< typename INDX >
void simplexD1linked<INDX>::construct
(
  INDX const a, INDX const b, INDX const na, INDX const nb
)
{
  pi[0] = a;
  pi[1] = b;
  ni[0] = na;
  ni[1] = nb;
}




#endif