DiscreteExteriorCalculusSolver.h

 
#pragma once

 
#if defined(DiscreteExteriorCalculusSolver_RECURSES)
#error Recursive header files inclusion detected in DiscreteExteriorCalculusSolver.h
#else // defined(DiscreteExteriorCalculusSolver_RECURSES)
#define DiscreteExteriorCalculusSolver_RECURSES
 
#if !defined DiscreteExteriorCalculusSolver_h
#define DiscreteExteriorCalculusSolver_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/base/ConstAlias.h"
#include "DGtal/base/Clone.h"
#include "DGtal/dec/KForm.h"
#include "DGtal/dec/LinearOperator.h"
 
namespace DGtal
{

  // template class DiscreteExteriorCalculusSolver
  template <typename TCalculus, typename TLinearAlgebraSolver, Order order_in, Duality duality_in, Order order_out, Duality duality_out>
  class DiscreteExteriorCalculusSolver
  {
    // ----------------------- Standard services ------------------------------
  public:
    BOOST_STATIC_ASSERT(( order_in  >= 0 ));
    BOOST_STATIC_ASSERT(( order_in  <= TCalculus::dimensionEmbedded ));
    BOOST_STATIC_ASSERT(( order_out >= 0 ));
    BOOST_STATIC_ASSERT(( order_out <= TCalculus::dimensionEmbedded ));
 
    typedef TCalculus Calculus;
    typedef TLinearAlgebraSolver LinearAlgebraSolver;
 
    typedef LinearOperator<Calculus, order_in, duality_in, order_out, duality_out> Operator;
    typedef KForm<Calculus, order_in, duality_in> SolutionKForm;
    typedef KForm<Calculus, order_out, duality_out> InputKForm;

    DiscreteExteriorCalculusSolver();
 
    // ----------------------- Interface --------------------------------------
  public:
    const Calculus* myCalculus;

    LinearAlgebraSolver myLinearAlgebraSolver;

    void selfDisplay(std::ostream& out) const;

    DiscreteExteriorCalculusSolver& compute(const Operator& linear_operator);

    SolutionKForm solve(const InputKForm& input_kform) const;

    bool isValid() const;
 
    // ------------------------- Private Data --------------------------------
  private:
 
    // ------------------------- Hidden services ------------------------------
  protected:
 
    // ------------------------- Internals ------------------------------------
  private:
 
  }; // end of class DiscreteExteriorCalculusSolver
 

  template <typename C, typename S, Order order_in, Duality duality_in, Order order_out, Duality duality_out>
  std::ostream&
  operator<<(std::ostream& out, const DiscreteExteriorCalculusSolver<C, S, order_in, duality_in, order_out, duality_out>& object);
 
} // namespace DGtal
 

// Includes inline functions.
#include "DGtal/dec/DiscreteExteriorCalculusSolver.ih"
 
//                                                                           //
 
#endif // !defined DiscreteExteriorCalculusSolver_h
 
#undef DiscreteExteriorCalculusSolver_RECURSES
#endif // else defined(DiscreteExteriorCalculusSolver_RECURSES)