00001 #ifndef _MY_MixedHybridCompressibleNewton_ 00002 #define _MY_MixedHybridCompressibleNewton_ 00003 00004 00017 class MixedHybridCompressibleNewton : public MixedHybridBase, public DynamicBase 00018 { 00019 private: 00020 OrthoMesh &m_mesh; 00022 Function1D &m_density; 00023 const VecDouble &m_vK; 00024 double m_dt; 00026 SparsityPattern m_sparsePattern; 00027 SparseMatrix<double> m_A; 00028 VecDouble m_B; 00029 VecDouble m_Pn; 00030 VecDouble m_Pnn; 00031 LinearSolver &m_solver; 00033 00034 00035 00036 int _debugLevel; 00037 00038 00039 HybridFaceBC m_bc; 00040 00041 00042 protected: 00043 00044 public: 00045 MixedHybridCompressibleNewton(OrthoMesh &mesh,Function3D &fPrescribedVelocities, Function3D &fPrescribedPression, const VecDouble &K, Function1D &density,LinearSolver &solver,int debugLevel); 00046 ~MixedHybridCompressibleNewton(); 00047 virtual void iterate(TransportBase &trans); 00048 00049 }; 00050 00051 #endif