MemoryBank.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 __MEMORY_BANK_HPP 
#define __MEMORY_BANK_HPP 

#include "Parallel_defines.hpp"
#include "Router.hpp"
#include "CapeMPI.hpp"
#include <stdio.h>


typedef struct
{
   long Size; 
   double Times[ MAX_TRANSMIT_QUANTITY ];
   long Orders[ MAX_TRANSMIT_QUANTITY ];
   double Steps[ MAX_TRANSMIT_QUANTITY ];
   double Variables[ MAX_TRANSMIT_QUANTITY ];
} TransmitDataStructure; //TDS

typedef struct
{
   CapeBoolean bNeedToWait;
   MPI_Request SizeReq;
   MPI_Request TimesReq;
   MPI_Request VariablesReq;
   MPI_Request OrdersReq;
   MPI_Request StepsReq;
} TransmitRequestsDataStructureRequests; //TRDS

class CMemoryBank
{
private:
   DianaDAESO* m_pOwner; 
   CapeLong m_lOwnerIndex; 
   CapeLong m_lReceiverIndex; 
   CapeArrayLong m_SendVariablesIndexes, m_ReceiveParametersIndexes; 

   CapeLong lSendNewIdx, lSendOldIdx, lReceiveOldIdx, lReceiveNewIdx;
   std::vector<TransmitDataStructure> TDSVector;
   
   CapeLong lSendReqIdx, lReceiveReqIdx;
   std::vector<TransmitRequestsDataStructureRequests> TRDSVector;

   /* Relative and absolute tolerance for convergence test */
   CapeDouble m_dRelTol, m_dAbsTol;
   CapeLong m_lMaxOrder;

   /* This values is used for identifying connections between two objects.
      Each memorybank could have just one connection.
      Identificators should be unique for each connection.*/
   CapeLong m_lConnectionsSendIdentificator, m_lConnectionsReceiveIdentificator; 
   
   CapeDouble m_dLastUpdateTime;

private:
   
   CapeArrayDouble GetValueAtTime( CapeDouble _Time, TransmitDataStructure* _Struct );

   void Initialize();

   void WaitFor( TransmitRequestsDataStructureRequests* _Requests );

public:
   CapeLong m_lIterationNumber; //the value of iteration on this window

   CapeBoolean m_bTempConvergence;

public:
   CMemoryBank();
   CMemoryBank( DianaDAESO* _pOwner, CapeLong _nOwnerIndex, CapeLong _nReceiverIndex, 
      CapeLong _lConnectionsSendIdentificator, CapeLong _lConnectionsReceivedIdentificator );
   ~CMemoryBank();

   void AddSendIndex( CapeLong _Index );


   void AddReceiveIndexes( CapeLong _ParIndex );


   void SetESO ( DianaDAESO* _pOwner );

   void SaveStateNew( double X, int NEQ, int KOLD, double* PHI, double* PSI );

   void SendData( CRouter* _Router );
   
   
   void ReceiveSize( CRouter* _Router );


   void ReceiveData( CRouter* _Router );


   CapeBoolean CheckConvergence();

   
   void SaveInitialState();
   

   void UpdateReceivedValues()
      throw (Common::Error::ECapeUnknown);

   void ResendWaveform();
};


#endif