PDE.hpp

Go to the documentation of this file.

//  This file is part of ff3d - http://www.freefem.org/ff3d
//  Copyright (C) 2001, 2002, 2003 Stéphane Del Pino

//  This program is free software; you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation; either version 2, or (at your option)
//  any later version.

//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.

//  You should have received a copy of the GNU General Public License
//  along with this program; if not, write to the Free Software Foundation,
//  Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  

//  $Id: PDE.hpp,v 1.5 2007/05/20 23:02:47 delpinux Exp $

// This class describes a PDE

#ifndef PDE_HPP
#define PDE_HPP

#include <VectorialPDEOperator.hpp>
#include <ScalarFunctionBuilder.hpp>
#include <ScalarFunctionBase.hpp>

class PDE
{
private:
  ConstReferenceCounting<VectorialPDEOperator>
  __vectorPDEOperator;          
  ConstReferenceCounting<ScalarFunctionBase>
  __secondMember;               
  PDE()
  {
    ;
  }
public:
  ReferenceCounting<PDE>
  operator * (const ConstReferenceCounting<ScalarFunctionBase>& c) const
  {
    ScalarFunctionBuilder functionBuilder;
    functionBuilder.setFunction(__secondMember);
    functionBuilder.setBinaryOperation(BinaryOperation::product,
                                       c);

    PDE* newPDE = new PDE();
    (*newPDE).__secondMember = functionBuilder.getBuiltFunction();

    (*newPDE).__vectorPDEOperator = (*__vectorPDEOperator) * c;
    return newPDE;
  }

  ConstReferenceCounting<ScalarFunctionBase>
  secondMember() const
  {
    return __secondMember;
  }

  ConstReferenceCounting<PDEOperatorSum>
  operator [] (const size_t& i) const
  {
    ASSERT (i<(*__vectorPDEOperator).size());
    return (*__vectorPDEOperator)[i];
  }

  const PDE& operator = (const PDE& pde)
  {
    __vectorPDEOperator=pde.__vectorPDEOperator;
    __secondMember = pde.__secondMember;
    return *this;
  }

  friend std::ostream&
  operator << (std::ostream& os,
               const PDE& pde)
  {
    os << *(pde.__vectorPDEOperator) << " = " << "pde.secondMember" << '\n';
    return os;
  }

  PDE(ConstReferenceCounting<VectorialPDEOperator> A,
      ConstReferenceCounting<ScalarFunctionBase> f)
    : __vectorPDEOperator(A),
      __secondMember(f)
  {
    ;
  }

  PDE(const PDE& pde)
    : __vectorPDEOperator(pde.__vectorPDEOperator),
      __secondMember(pde.__secondMember)
  {
    ;
  }

  ~PDE()
  {
    ;
  }
};

#endif // PDE_HPP

---