~simula-acdc/dolfin/dolfin-spacetime

#ifndef __COMBUSTION_KRONECKER_H

#define __COMBUSTION_KRONECKER_H

 

// Include dolfin libraries

#include <dolfin/tensorproduct/TensorProductFunctionSpace.h>

 

// Include the generated code (includes a lot of other stuff)

#include "Combustion.h"

 

namespace CombustionKronecker

{

 

class KroneckerForm_0 : public dolfin::TensorProductForm

{

  public:

 

    // Constructor

    KroneckerForm_0(boost::shared_ptr<const dolfin::TensorProductFunctionSpace> V1, boost::shared_ptr<const dolfin::TensorProductFunctionSpace> V0)

     : dolfin::TensorProductForm(2, 0, 2, 4)

    {

      // Store function spaces

      _tensor_spaces[0] = V0;

      _tensor_spaces[1] = V1;

 

      // Create forms in flattened structure (a0, a1, b0, b1, c0, c1) etc.

      _forms[0] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_0(V1->factor_space(0),V0->factor_space(0)));

      _forms[2] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_2(V1->factor_space(0),V0->factor_space(0)));

      _forms[4] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_4(V1->factor_space(0),V0->factor_space(0)));

      _forms[6] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_6(V1->factor_space(0),V0->factor_space(0)));

      _forms[1] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_1(V1->factor_space(1),V0->factor_space(1)));

      _forms[3] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_3(V1->factor_space(1),V0->factor_space(1)));

      _forms[5] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_5(V1->factor_space(1),V0->factor_space(1)));

      _forms[7] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_7(V1->factor_space(1),V0->factor_space(1)));

 

    }

 

    // Destructor

    ~KroneckerForm_0()

    {}

 

};

 

class KroneckerForm_1 : public dolfin::TensorProductForm

{

  public:

 

    // Constructor

    KroneckerForm_1(boost::shared_ptr<const dolfin::TensorProductFunctionSpace> V0)

     : dolfin::TensorProductForm(1, 0, 2, 1)

    {

      // Store function spaces

      _tensor_spaces[0] = V0;

 

      // Create forms in flattened structure (a0, a1, b0, b1, c0, c1) etc.

      _forms[0] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_8(V0->factor_space(0)));

      _forms[1] = boost::shared_ptr<dolfin::Form>(new Combustion::Form_9(V0->factor_space(1)));

 

    }

 

    // Destructor

    ~KroneckerForm_1()

    {}

 

};

 

typedef KroneckerForm_1 LinearForm;

typedef KroneckerForm_0 BilinearForm;

 

}

#endif