00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023 #if !defined(SCHURLAPOP_HPP)
00024 #define SCHURLAPOP_HPP
00025
00026 #include "linop.hpp"
00027 #include "elem2d.hpp"
00028 #include "massop.hpp"
00029 #include "stokesop.hpp"
00030
00031 class SchurLapOp: public LinOp
00032 {
00033 public:
00034 SchurLapOp();
00035 SchurLapOp(Elem2D *ve, Elem2D *pe);
00036
00037 ~SchurLapOp();
00038
00039 GVector operator*(GVector) ;
00040 GVector *GetMassMatrix();
00041 void OpVec_prod(GVector &x, GVector &ret);
00042 void SetConst(GDOUBLE c);
00043 void SetElem(Elem2D *ve, Elem2D *pe);
00044 void SetVBdyData(GINT idir, GVector *vbdy);
00045 void DoCentralOp(GBOOL bop);
00046 GVector *GetiMass();
00047 Elem2D *GetElem();
00048 #if 0
00049 Elem2D *GetPElem();
00050 #endif
00051
00052
00053
00054 private:
00055 GBOOL QuadOp(GVector *x, GVector *newx);
00056 GBOOL TriangleOp(GVector *x, GVector *newx);
00057 void DeleteDynamic();
00058 void Reset();
00059 void Init(Elem2D *ve, Elem2D *pe);
00060
00061 GBOOL bCentralOp_;
00062 GDOUBLE mConst;
00063 Elem2D *velem;
00064 Elem2D *pelem;
00065 GVector *vtmpp0;
00066 GVector *vtmpv0;
00067 GVector *vtmpv1;
00068 GVector *v1bdy;
00069 GVector *v2bdy;
00070 GVector *M;
00071 GVector *iM;
00072 StokesOp *pD1;
00073 StokesOp *pD2;
00074 StokesOp *pDT1;
00075 StokesOp *pDT2;
00076
00077 };
00078
00079 #endif