HopfPointContinuation.hpp

/* -----------------------------------------------------------------------------
 *  Diana process modelling, simulation and analysis software
 *  Copyright (c) 2005, Michael Krasnyk
 *  All rights reserved.

 *  This file is a part of Diana process modelling, simulation and analysis software
 *
 *  Diana 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 (see accompanying file LICENSE)
 * -----------------------------------------------------------------------------
 *  $Id: HopfPointContinuation.hpp 8160 2009-03-16 18:36:30Z miha $
 * -----------------------------------------------------------------------------
 *  Description:
 */

#ifndef HOPF_POINT_CONTINUATION__HPP
#define HOPF_POINT_CONTINUATION__HPP


#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#include <Diana.hpp>
#include <Fortran.hpp>
#include <DianaContinuation.hpp>
#include <vector>

namespace Diana {

  class HopfPointContinuation : public DianaContinuation {
  public:
    HopfPointContinuation(Numeric::Solvers::Model::ICapeNumericModel* _pModel);
    virtual ~HopfPointContinuation();

    Numeric::Solvers::Solver::SolveReturn Solve()
      throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadInvOrder,
             Common::Error::ECapeSolvingError, Common::Error::ECapeOutOfResources);

    ContiReturn Continuate()
      throw (Common::Error::ECapeUnknown, Common::Error::ECapeSolvingError);
  
    Common::Types::CapeString GetVariableName(Common::Types::CapeLong index)
      throw (Common::Error::ECapeUnknown);

  protected:
    int getModelEquations();

    int getModelVariables();

    void packModelBoundaries(Common::Types::CapeArrayDouble& arrMinY,
                             Common::Types::CapeArrayDouble& arrMaxY);

    void packModelVector(Common::Types::CapeArrayDouble& modelVector);

    void unpackModelVector(const Common::Types::CapeArrayDouble& modelVector);

    void packModelResidual(Common::Types::CapeArrayDouble& modelResidual);

    void packModelJacobian(Diana::DianaSparseArray& Jaco);

    CapeBoolean calcInitialEigenvector();
    CapeBoolean calcEigenvalues();
    CapeBoolean calcEigenvaluesLapack();
    CapeBoolean calcEigenvaluesArpack();

  private:
    Common::Types::CapeLong               nESOSize;
    Common::Types::CapeArrayDouble        arrReEigenVector;
    Common::Types::CapeArrayDouble        arrImEigenVector;
    Common::Types::CapeDouble             dblOmega;

    const Common::Types::CapeLong   nEigMonitorMax;
    Common::Types::CapeLong         nEigMonitor;
    Common::Types::CapeLong         nEigConverged;
    Common::Types::CapeArrayDouble  arrReEigenValues;
    Common::Types::CapeArrayDouble  arrImEigenValues;

  private:
    std::vector<fortrand>           lapA;
    std::vector<fortrand>           lapB;
    std::vector<fortrand>           lapAlphaR;
    std::vector<fortrand>           lapAlphaI;
    std::vector<fortrand>           lapBeta;
    std::vector<fortrand>           lapVl;
    std::vector<fortrand>           lapVr;
    std::vector<fortrand>           lapWork;

  private:
    std::vector<fortrand>           arpResid;
    std::vector<fortrand>           arpV;
    std::vector<fortrand>           arpWorkd;
    std::vector<fortrand>           arpWorkl;
    std::vector<fortrand>           arpTmpX;
    std::vector<fortranl>           arpSelect;
    std::vector<fortrand>           arpWorkev;

    Diana::DianaSparseArray         dsaArpAsM;
    Diana::DianaSparseSolver        dssArpM;
    Diana::DianaSparseSolver        dssArpAsM;
  };

};

#endif // HOPF_POINT_CONTINUATION__HPP