00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022 #if !defined(UZAWALAPOP_HPP)
00023 #define UZAWALAPOP_HPP
00024
00025 #include "linop.hpp"
00026 #include "elem2d.hpp"
00027 #include "massop.hpp"
00028 #include "stokesop.hpp"
00029 #include "ilaplacianop.hpp"
00030
00031 class UzawaLapOp: public LinOp
00032 {
00033 public:
00034 UzawaLapOp();
00035 UzawaLapOp(Elem2D *ve, Elem2D *pe);
00036
00037 ~UzawaLapOp();
00038
00039 GVector operator*(GVector) ;
00040 void OpVec_prod(GVector &x, GVector &ret);
00041 void SetConst(GDOUBLE c);
00042 void SetElem(Elem2D *ve, Elem2D *pe);
00043 void SetMask(GVector *m);
00044 void SetVBdyData (GINT idir, GVector *vb, GINT *indices);
00045
00046
00047
00048 private:
00049 GBOOL QuadOp(GVector *x, GVector *newx);
00050 GBOOL TriangleOp(GVector *x, GVector *newx);
00051 void DeleteDynamic();
00052 void Init(Elem2D *ve, Elem2D *pe);
00053
00054 GDOUBLE mConst;
00055 GINT *ibdy;
00056 Elem2D *velem;
00057 Elem2D *pelem;
00058 GVector *vtmpp0;
00059 GVector *vtmpv0;
00060 GVector *vtmpv1;
00061 GVector *vbdy[3];
00062 GVector *mask;
00063 iLaplacianOp *iL;
00064 StokesOp *pD1;
00065 StokesOp *pD2;
00066 StokesOp *pDT1;
00067 StokesOp *pDT2;
00068
00069 };
00070
00071 #endif