ParallelSolver.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$
 * -----------------------------------------------------------------------------
 *  Description:
 */


#ifndef __PARALLSOLVER_HPP
#define __PARALLSOLVER_HPP

#include "Router.hpp"
#include "NumericObject.hpp"
#include "MemoryManager.hpp"
#include "Convertor.hpp"
#include <Diana.hpp>

class CParallelSolver : virtual public IParallelSolver
{
private:
   int m_nMyRank; 
   int m_nProcQuantity; 
   CRouter m_Router; 
   CMemoryManager m_MemoryManager; 

   CapeString m_szConnectionsFileName;
   std::vector<CNumericObject*> m_Objects; 
   CapeArrayBoolean m_ObjectsConvergences; 

   CapeArrayLong m_ModelDisplacement;

   CapeLong m_lVerboseLevel;
   CapeBoolean m_bNeedToConsistentInitialize;
   
   CapeLong m_lErrorStatus; 
private:
   void InitialLoadBalancing()
        throw (Common::Error::ECapeUnknown, Common::Error::ECapeInvalidArgument);
 
   void TransmitData();


   void SaveState()
      throw (Common::Error::ECapeUnknown);
   void LoadState()
      throw (Common::Error::ECapeUnknown);

   CapeBoolean CheckConvergence();
   
   void ReduceError();

   CapeLong MoveNumericObject( CapeLong _nObjectIndex, 
      CapeLong _nSourcePE, CapeLong _nDestinationPE );

   CapeDouble CalculateNewWindowSize( CapeLong _lIterationNumber, CapeDouble _dOldWindowSize,
      CapeBoolean _bWasConverged );

   void ShowStat( CapeBoolean _bShowTotal );
   

public:
   CParallelSolver( const Common::Types::CapeString& _szConnectionsFileName );
   ~CParallelSolver();

   CapeLong FindInitialConditions( CapeDouble _dStartTime, CapeDouble _dEndTime,
      CapeDouble _dInitialStepSize = 0 )
      throw (Common::Error::ECapeUnknown, Common::Error::ECapeInvalidArgument);

   void Run( CapeDouble _dInitialWindowSize, CapeDouble _dEndTime, 
      CapeDouble _dInitialStepSize = 0 )
      throw (Common::Error::ECapeUnknown, Common::Error::ECapeInvalidArgument);

   Numeric::Solvers::Solver::ICapeNumericSolver* 
      AddModel( Numeric::Solvers::Model::ICapeNumericModel* _Model, 
      const Common::Types::CapeString& ModelName )
      throw (Common::Error::ECapeUnknown, Common::Error::ECapeInvalidArgument);


   void ShowInfo();

   Common::Types::CapeArrayString GetAllModelNames();
   Common::Types::CapeArrayString GetAllESONames();

};

#endif