00001 #ifndef MULTIPHASEC_H
00002 #define MULTIPHASEC_H
00003
00004 #include "genfile.h"
00005
00006 class multiphc
00007 {
00008 public:
00009 multiphc();
00010 ~multiphc();
00011
00012 void matcond_u (matrix &d,long ri,long ci,long ipp);
00013 void matcap_u (matrix &d,long ri,long ci,long ipp);
00014 void matcond1d_u (matrix &d,long ri,long ci,long ipp);
00015 void matcond2d_u (matrix &d,long ri,long ci,long ipp);
00016 void matcond3d_u (matrix &d,long ri,long ci,long ipp);
00017 void matcap1d_u (matrix &d,long ri,long ci,long ipp);
00018 void matcap2d_u (matrix &d,long ri,long ci,long ipp);
00019 void matcap3d_u (matrix &d,long ri,long ci,long ipp);
00020
00021 void matcond_l (matrix &d,long ri,long ci,long ipp);
00022 void matcap_l (matrix &d,long ri,long ci,long ipp);
00023 void matcond1d_l (matrix &d,long ri,long ci,long ipp);
00024 void matcond2d_l (matrix &d,long ri,long ci,long ipp);
00025 void matcond3d_l (matrix &d,long ri,long ci,long ipp);
00026 void matcap1d_l (matrix &d,long ri,long ci,long ipp);
00027 void matcap2d_l (matrix &d,long ri,long ci,long ipp);
00028 void matcap3d_l (matrix &d,long ri,long ci,long ipp);
00029
00030 void gaspress_check(double pc,double &pg,double t,long ipp);
00031 void cappress_check(double &pc,double pg,double t,long ipp);
00032 void cappress_stop(double &pc,double pg,double t,long ipp);
00033
00034 void values_correction (vector &nv);
00035 double get_kcu(double pc,double pg,double t,long ipp);
00036 double get_capcu(double pc,double pg,double t,long ipp);
00037 double get_kgu(double pc,double pg,double t,long ipp);
00038 double get_capgu(double pc,double pg,double t,long ipp);
00039 double get_ktu(double pc,double pg,double t,long ipp);
00040 double get_captu(double pc,double pg,double t,long ipp);
00041
00042 double get_kug(double pc,double pg,double t,long ipp);
00043 double get_capug(double pc,double pg,double t,long ipp);
00044 double get_kuc(double pc,double pg,double t,long ipp);
00045 double get_capuc(double pc,double pg,double t,long ipp);
00046 double get_kut(double pc,double pg,double t,long ipp);
00047 double get_caput(double pc,double pg,double t,long ipp);
00048
00049 void rhs_u1 (matrix &d,long ri,long ci,long ipp);
00050 void rhs1d_u1 (matrix &d,long ri,long ci,long ipp);
00051 void rhs2d_u1 (matrix &d,long ri,long ci,long ipp);
00052 void rhs3d_u1 (matrix &d,long ri,long ci,long ipp);
00053
00054 double get_fuc1(double pc,double pg,double t,long ipp);
00055 double get_fug1(double pc,double pg,double t,long ipp);
00056 double get_fut1(double pc,double pg,double t,long ipp);
00057
00058 void rhs_u2 (matrix &d,long ri,long ci,long ipp);
00059 void rhs1d_u2 (matrix &d,long ri,long ci,long ipp);
00060 void rhs2d_u2 (matrix &d,long ri,long ci,long ipp);
00061 void rhs3d_u2 (matrix &d,long ri,long ci,long ipp);
00062
00063 double get_fu2(double pc,double pg,double t,long ipp);
00064
00065 private:
00066
00067 double scale_pc,scale_pg,scale_t;
00068
00069 };
00070
00071 #endif