us_astfem_rsa.h

Go to the documentation of this file.

#ifndef US_ASTFEM_RSA_H
#define US_ASTFEM_RSA_H

#include <QtCore>
#include "us_extern.h"
#include "us_model.h"
#include "us_simparms.h"
#include "us_dataIO.h"
#include "us_astfem_math.h"
#include "us_stiffbase.h"

#ifndef DbgLv
#define DbgLv(a) if(dbg_level>=a)qDebug() 
#endif

class US_UTIL_EXTERN US_Astfem_RSA : public QObject
{
   Q_OBJECT

   public:
      US_Astfem_RSA( US_Model&, US_SimulationParameters&, QObject* = 0 );
   
      int  calculate           ( US_DataIO::RawData& );

      void setTimeCorrection   ( bool flag ){ time_correction = flag; }; 

      void setTimeInterpolation( bool flag ){ use_time        = flag; };

      void setStopFlag         ( bool flag ){ stopFlag        = flag; };    

      void set_movie_flag      ( bool flag ){ show_movie      = flag; };

      void set_simout_flag     ( bool flag ){ simout_flag     = flag; };

      void set_debug_flag      ( int  flag ){ dbg_level       = flag; };

   signals:
      void new_scan         ( QVector< double >*, double* );

      void new_time         ( double );

      void current_component( int    );

      void current_speed    ( int    );

      void calc_start       ( int    );

      void calc_progress    ( int    );

      void calc_done        ( void   );

   private:
      bool stopFlag;          
      bool time_correction;   
      
      bool use_time;          
      
      bool show_movie;
      bool simout_flag;
      US_StiffBase stfb0;

      double last_time;      

      double w2t_integral;    
      int    Nx;              
      int    dbg_level;
      
      US_AstfemMath::AstFemParameters af_params;
      US_AstfemMath::MfemData         af_data;
      US_AstfemMath::MfemInitial      af_c0;

      double*                                 xA; //<! x values array
      QVector< double >                       x;  //<! Radii of grid points
      QVector< US_AstfemMath::ReactionGroup > rg;
      US_Model&                               system;
      US_SimulationParameters&                simparams;

      // Functions
      void   update_assocv  ( void );
      void   adjust_limits  ( int speed );
      double stretch        ( double*, int );
      void   initialize_rg  ( void );
      void   initialize_conc( int, US_AstfemMath::MfemInitial&, bool );

      int    calculate_ni   ( double, double, US_AstfemMath::MfemInitial&,
                              US_AstfemMath::MfemData&, bool );
      void   mesh_gen       ( QVector< double >&, int );
      void   mesh_gen_s_pos ( const QVector< double >& );
      void   mesh_gen_s_neg ( const QVector< double >& );
      void   mesh_gen_RefL  ( int, int );
      
      void   ComputeCoefMatrixFixedMesh( double, double, double**, double** );
      void   decompose      ( US_AstfemMath::MfemInitial* );

      void   ComputeCoefMatrixMovingMeshR( double, double, double**, double** );
      void   ComputeCoefMatrixMovingMeshL( double, double, double**, double** );
             
      void   ReactionOneStep_Euler_imp   ( int, double**, double );
             
      void   Reaction_dydt    ( double*, double*  );
      void   Reaction_dfdy    ( double*, double** );
             
      int    calculate_ra2    ( double, double, US_AstfemMath::MfemInitial*,
                                US_AstfemMath::MfemData&, bool );         

      void   GlobalStiff      ( double*, double**, double**,
                                double, double );

      void   load_mfem_data ( US_DataIO::RawData&, US_AstfemMath::MfemData& );         
      void   store_mfem_data( US_DataIO::RawData&, US_AstfemMath::MfemData& );         

#ifdef NEVER
      void GlobalStiff_ellam( QVector <double> *, double **, double **,
                              double, double );
      void adjust_grid( int /*old speed*/, int /*new speed*/,
                        QVector <double> * /*radial grid*/ );
#endif
};


#endif