e2.h

Go to the documentation of this file.

//
// Copyright (c) 2014 Limit Point Systems, Inc. 
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//   http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
//


#ifndef E2_H
#define E2_H

#ifndef SHEAF_DLL_SPEC_H
#include "SheafSystem/sheaf_dll_spec.h"
#endif

#ifndef ED_H
#include "SheafSystem/ed.h"
#endif

namespace fiber_bundle
{

//==============================================================================
// CLASS E2_ROW_DOFS_TYPE
//==============================================================================

template <typename T>
class e2_row_dofs_type
{
public:

  typedef T dof_type;

  typedef e2_row_dofs_type<T> row_dofs_type;

  static const std::string& static_class_name();
  
  static int d();

  T& operator[] (int xindex);

  const T& operator[] (int xindex) const;

  e2_row_dofs_type<T>* clone() const;

  e2_row_dofs_type<T>* clone(bool xnew_state, bool xauto_access) const
  {
    return clone();
  }

  T components[2];

protected:
private:

};


//==============================================================================
// CLASS E2_LITE
//==============================================================================

class SHEAF_DLL_SPEC e2_lite : public ed_lite
{

  //============================================================================
  //============================================================================

public:

  // Typedefs:

  typedef e2_row_dofs_type<double> row_dofs_type;

  e2_lite();

  e2_lite(const e2_lite& xother);

  e2_lite& operator=(const e2_lite& xother);

  virtual ~e2_lite();

  e2_lite(const row_dofs_type& xrow_dofs);

  e2_lite& operator=(const row_dofs_type& xrow_dofs);

  e2_lite(const value_type& x_comp, const value_type& y_comp);

  using ed_lite::put_components;

  virtual void put_components(const value_type& x_comp,
                              const value_type& y_comp);

  operator row_dofs_type& ();

  operator const row_dofs_type& () const;

protected:

  row_dofs_type _row_dofs;

private:



  //============================================================================
  //============================================================================

public:
protected:
private:



  //============================================================================
  //============================================================================

public:
protected:
private:



  //============================================================================
  //============================================================================

public:
protected:
private:



  //============================================================================
  //============================================================================

public:

  typedef e2_lite vector_space_type;
  
  virtual int dd() const;

protected:
private:



  //============================================================================
  //============================================================================

public:

  bool operator==(const e2_lite& xother) const;

  static const tp_lite& static_tp_prototype(int xp);

  virtual const tp_lite& tp_prototype(int xp) const;

  static const atp_lite& static_atp_prototype(int xp);

  virtual const atp_lite& atp_prototype(int xp) const;

  static const stp_lite& static_stp_prototype(int xp);

  virtual const stp_lite& stp_prototype(int xp) const;

protected:
private:



  //============================================================================
  //============================================================================

public:
protected:
private:



  //============================================================================
  //============================================================================

public:

  virtual const std::string& class_name() const;

  static const std::string& static_class_name();
  
  virtual e2_lite* clone() const;

  virtual e2_lite* clone(bool xnew_state, bool xauto_access) const
  {
    return clone();
  }

protected:
private:



  //============================================================================
  //============================================================================

public:

  virtual bool is_ancestor_of(const any_lite& xother) const;

  bool invariant() const;

  virtual void* row_dofs();

  virtual const void* row_dofs() const;

protected:
private:

};


//==============================================================================
// CLASS E2
//==============================================================================

class tp_space;
class tp;
class atp;
class stp;

class SHEAF_DLL_SPEC e2 : public ed
{
  // ===========================================================
  // ===========================================================

public:

  static const poset_path& standard_schema_path();

  static void make_standard_schema(namespace_poset& xns);

  static host_type& new_host(namespace_type& xns, 
                 const poset_path& xhost_path, 
                 const poset_path& xschema_path,
                 const poset_path& xscalar_space_path,
                 bool xauto_access);

  static host_type& standard_host(namespace_type& xns, const std::string& xsuffix, bool xauto_access);

protected:

private:



  //============================================================================
  //============================================================================

public:

  // Typedefs:

  typedef e2_row_dofs_type<double> row_dofs_type;

  typedef e2_lite volatile_type;

  typedef e2 vector_space_type;

  e2();

  // New handle/old state constructors.

  e2(const poset_state_handle* xhost, pod_index_type xhub_id);

  e2(const poset_state_handle* xhost, const scoped_index& xid);

  e2(const poset_state_handle* xhost, const std::string& xname);

  explicit e2(abstract_poset_member* xother);

  e2(poset_state_handle* xhost, bool xauto_access = true);

  e2(poset_state_handle& xhost,
     const row_dofs_type& xrdt,
     bool xauto_access = true);

  e2& operator=(const row_dofs_type& xrdt);

  virtual e2& operator=(const abstract_poset_member& xother);

  e2& operator=(const e2& xother);

  ~e2();

  virtual const volatile_type& lite_prototype() const;

  virtual volatile_type* lite_type() const;


protected:
private:



  //============================================================================
  //============================================================================

public:
protected:
private:



  //============================================================================
  //============================================================================

public:
protected:
private:



  //===========================================================================
  //===========================================================================

public:
protected:
private:



  //===========================================================================
  //===========================================================================

public:
protected:
private:



  //===========================================================================
  //===========================================================================

public:

  using ed::put_components;

  virtual void put_components(dof_type x_comp, dof_type y_comp);

  operator row_dofs_type& ();

  operator const row_dofs_type& () const;

protected:
private:



  //===========================================================================
  //===========================================================================

public:

  static tp* new_tp(tp_space& xvector_space, int xp);

  static atp* new_atp(tp_space& xvector_space, int xp);

  static stp* new_stp(tp_space& xvector_space, int xp);


  virtual tp* new_tp(int xp, bool xauto_access) const;

  virtual atp* new_atp(int xp, bool xauto_access) const;

  virtual stp* new_stp(int xp, bool xauto_access) const;

protected:
private:



  //============================================================================
  //============================================================================

public:

  virtual const std::string& class_name() const;

  static const std::string& static_class_name();
  
  virtual e2* clone() const;

  inline e2* clone(bool xnew_state, bool xauto_access = true) const
  {
    return static_cast<e2*>(ed::clone(xnew_state, xauto_access));
  }

protected:
private:



  //============================================================================
  //============================================================================

public:

  bool is_ancestor_of(const any* xother) const;

  bool invariant() const;

protected:
private:

};

//==============================================================================
// NON-MEMBER FUNCTIONS
//==============================================================================

template <typename T>
std::ostream& operator<<(std::ostream& xos, const e2_row_dofs_type<T>& xrdt);

namespace e2_algebra
{

SHEAF_DLL_SPEC vd_value_type perp_dot(const e2_lite& x0, const e2_lite& x1);

} // namespace e2_algebra

//==============================================================================
// TP FACET
//==============================================================================

// Forward declarations needed for specialization:

class at0_lite;
class at2_e2_lite;
class at3_e2_lite;
class st2_e2_lite;
class st3_e2_lite;
class st4_e2_lite;
class t2_e2_lite;
class t3_e2_lite;
class t4_e2_lite;

class at0;
class at2_e2;
class at3_e2;
class st2_e2;
class st3_e2;
class st4_e2;
class t2_e2;
class t3_e2;
class t4_e2;

template<> class tensor_traits<0, e2_lite>
{
public:
  typedef at0_lite tp_type;

  typedef at0_lite atp_type;

  typedef at0_lite stp_type;
};
  
template<> class tensor_traits<1, e2_lite>
{
public:
  typedef e2_lite tp_type;

  typedef e2_lite atp_type;

  typedef e2_lite stp_type;
};
  
template<> class tensor_traits<2, e2_lite>
{
public:
  typedef t2_e2_lite tp_type;

  typedef at2_e2_lite atp_type;

  typedef st2_e2_lite stp_type;
};
  
template<> class tensor_traits<3, e2_lite>
{
public:
  typedef t3_e2_lite tp_type;

  typedef at3_e2_lite atp_type;

  typedef st3_e2_lite stp_type;
};
  
template<> class tensor_traits<4, e2_lite>
{
public:
  typedef t4_e2_lite tp_type;

  typedef void atp_type;

  typedef st4_e2_lite stp_type;
};

//==============================================================================

template<> class tensor_traits<0, e2>
{
public:
  typedef at0 tp_type;

  typedef at0 atp_type;

  typedef at0 stp_type;
};
  
template<> class tensor_traits<1, e2>
{
public:
  typedef e2 tp_type;

  typedef e2 atp_type;

  typedef e2 stp_type;
};
  
template<> class tensor_traits<2, e2>
{
public:
  typedef t2_e2 tp_type;

  typedef at2_e2 atp_type;

  typedef st2_e2 stp_type;
};
  
template<> class tensor_traits<3, e2>
{
public:
  typedef t3_e2 tp_type;

  typedef at3_e2 atp_type;

  typedef st3_e2 stp_type;
};
  
template<> class tensor_traits<4, e2>
{
public:
  typedef t4_e2 tp_type;

  typedef void atp_type;

  typedef st4_e2 stp_type;
};

} // namespace fiber_bundle

#endif // ifndef E2_H