qf11to25.cpp

FreeFem++可以生成高质量的有限元网格。可以用于流体力学

QuadraturePoint(0.0173070172095/2,0.0918740717058,0.0395513001973),
QuadraturePoint(0.0174524546493/2,0.8604888296191,0.0966224057079),
QuadraturePoint(0.0177217222159/2,0.0439842178673,0.8561886349107),
QuadraturePoint(0.0282824024023/2,0.2011017606735,0.7449115835626),
QuadraturePoint(0.0284996712488/2,0.7449993726263,0.0536865638166),
QuadraturePoint(0.0285005646539/2,0.0532186641310,0.1963754275935),
QuadraturePoint(0.0300647223478/2,0.7453984647401,0.1982065805550),
QuadraturePoint(0.0302031277082/2,0.1957289932876,0.0555713833156),
QuadraturePoint(0.0303987136077/2,0.1092532057988,0.6100036182413),
QuadraturePoint(0.0305668796074/2,0.0567625702001,0.7409121894959),
QuadraturePoint(0.0306067413002/2,0.0483837933475,0.6075135660978),
QuadraturePoint(0.0309330068201/2,0.1080612809760,0.1122081510437),
QuadraturePoint(0.0309773820835/2,0.6185605900991,0.2698753703035),
QuadraturePoint(0.0313146250545/2,0.7721296013497,0.1114117395333),
QuadraturePoint(0.0313573493392/2,0.6115734801133,0.3389367677931),
QuadraturePoint(0.0314320469287/2,0.3381326103376,0.0494693938787),
QuadraturePoint(0.0315182143894/2,0.1173084128254,0.7696451309795),
QuadraturePoint(0.0324248137985/2,0.2674551260596,0.1115718808154),
QuadraturePoint(0.0347512152386/2,0.6542100160026,0.1906548314700),
QuadraturePoint(0.0350393454927/2,0.0538297481158,0.3358616826849),
QuadraturePoint(0.0350717420310/2,0.1848840324117,0.1551831523851),
QuadraturePoint(0.0352129215334/2,0.3376267104744,0.6081402596294),
QuadraturePoint(0.0352615504981/2,0.6067102034499,0.0542632795598),
QuadraturePoint(0.0366403220343/2,0.4612614085496,0.0688176670722),
QuadraturePoint(0.0367733107670/2,0.1525465365671,0.6510240845749),
QuadraturePoint(0.0371675662937/2,0.0700582543543,0.4661904392742),
QuadraturePoint(0.0373371571606/2,0.4704201379032,0.4634826455353),
QuadraturePoint(0.0403973346588/2,0.1216461693746,0.2381494875516),
QuadraturePoint(0.0413580040638/2,0.6371404052702,0.1238399384513),
QuadraturePoint(0.0421957791870/2,0.2379904515119,0.6370216452326),
QuadraturePoint(0.0495451004037/2,0.1483929857177,0.4894188577780),
QuadraturePoint(0.0500419261141/2,0.3598069571550,0.1452880866253),
QuadraturePoint(0.0505794587115/2,0.4941441055095,0.3610216383818),
QuadraturePoint(0.0520037210188/2,0.1440630687981,0.3513508341887),
QuadraturePoint(0.0521533567886/2,0.5019764440004,0.1435491663293),
QuadraturePoint(0.0524899152358/2,0.3555423834298,0.5016491599502),
QuadraturePoint(0.0599159762516/2,0.2443439540771,0.2406052129104),
QuadraturePoint(0.0599609997426/2,0.2437064989342,0.5109017277055),
QuadraturePoint(0.0599915272129/2,0.5122200807321,0.2452737973543),
QuadraturePoint(0.0634133183449/2,0.2526038315178,0.3700319555094),
QuadraturePoint(0.0635311861108/2,0.3759895652851,0.2505406611631),
QuadraturePoint(0.0637206605672/2,0.3729077987144,0.3753750277549),
};
const QuadratureFormular QuadratureFormular_T_20(20 ,78,P_QuadratureFormular_T_20);


static  QuadraturePoint P_QuadratureFormular_T_21[] = {
QuadraturePoint(0.0006704436439/2,0.0035524391922,0.0035524391922),
QuadraturePoint(0.0006704436439/2,0.0035524391922,0.9928951216156),
QuadraturePoint(0.0006704436439/2,0.9928951216156,0.0035524391922),
QuadraturePoint(0.0045472608074/2,0.9553548273730,0.0087898929093),
QuadraturePoint(0.0045472608074/2,0.0358552797177,0.0087898929093),
QuadraturePoint(0.0045472608074/2,0.9553548273730,0.0358552797177),
QuadraturePoint(0.0045472608074/2,0.0087898929093,0.0358552797177),
QuadraturePoint(0.0045472608074/2,0.0087898929093,0.9553548273730),
QuadraturePoint(0.0045472608074/2,0.0358552797177,0.9553548273730),
QuadraturePoint(0.0052077585320/2,0.8865264879047,0.1082329745017),
QuadraturePoint(0.0052077585320/2,0.8865264879047,0.0052405375935),
QuadraturePoint(0.0052077585320/2,0.0052405375935,0.1082329745017),
QuadraturePoint(0.0052077585320/2,0.0052405375935,0.8865264879047),
QuadraturePoint(0.0052077585320/2,0.1082329745017,0.8865264879047),
QuadraturePoint(0.0052077585320/2,0.1082329745017,0.0052405375935),
QuadraturePoint(0.0065435432887/2,0.0466397432150,0.9067205135700),
QuadraturePoint(0.0065435432887/2,0.0466397432150,0.0466397432150),
QuadraturePoint(0.0065435432887/2,0.9067205135700,0.0466397432150),
QuadraturePoint(0.0092737841533/2,0.2075720456946,0.0082759241284),
QuadraturePoint(0.0092737841533/2,0.2075720456946,0.7841520301770),
QuadraturePoint(0.0092737841533/2,0.7841520301770,0.2075720456946),
QuadraturePoint(0.0092737841533/2,0.0082759241284,0.7841520301770),
QuadraturePoint(0.0092737841533/2,0.0082759241284,0.2075720456946),
QuadraturePoint(0.0092737841533/2,0.7841520301770,0.0082759241284),
QuadraturePoint(0.0095937782623/2,0.0858119489725,0.0314836947701),
QuadraturePoint(0.0095937782623/2,0.8827043562574,0.0314836947701),
QuadraturePoint(0.0095937782623/2,0.0314836947701,0.0858119489725),
QuadraturePoint(0.0095937782623/2,0.0858119489725,0.8827043562574),
QuadraturePoint(0.0095937782623/2,0.8827043562574,0.0858119489725),
QuadraturePoint(0.0095937782623/2,0.0314836947701,0.8827043562574),
QuadraturePoint(0.0114247809167/2,0.6688778233826,0.0095150760625),
QuadraturePoint(0.0114247809167/2,0.0095150760625,0.3216071005550),
QuadraturePoint(0.0114247809167/2,0.0095150760625,0.6688778233826),
QuadraturePoint(0.0114247809167/2,0.6688778233826,0.3216071005550),
QuadraturePoint(0.0114247809167/2,0.3216071005550,0.6688778233826),
QuadraturePoint(0.0114247809167/2,0.3216071005550,0.0095150760625),
QuadraturePoint(0.0117216964174/2,0.4379999543113,0.0099859785681),
QuadraturePoint(0.0117216964174/2,0.0099859785681,0.5520140671206),
QuadraturePoint(0.0117216964174/2,0.4379999543113,0.5520140671206),
QuadraturePoint(0.0117216964174/2,0.0099859785681,0.4379999543113),
QuadraturePoint(0.0117216964174/2,0.5520140671206,0.4379999543113),
QuadraturePoint(0.0117216964174/2,0.5520140671206,0.0099859785681),
QuadraturePoint(0.0188197155232/2,0.7974931072148,0.0405093994119),
QuadraturePoint(0.0188197155232/2,0.0405093994119,0.1619974933734),
QuadraturePoint(0.0188197155232/2,0.0405093994119,0.7974931072148),
QuadraturePoint(0.0188197155232/2,0.1619974933734,0.7974931072148),
QuadraturePoint(0.0188197155232/2,0.7974931072148,0.1619974933734),
QuadraturePoint(0.0188197155232/2,0.1619974933734,0.0405093994119),
QuadraturePoint(0.0235260980271/2,0.3864215551955,0.3864215551955),
QuadraturePoint(0.0235260980271/2,0.3864215551955,0.2271568896090),
QuadraturePoint(0.0235260980271/2,0.2271568896090,0.3864215551955),
QuadraturePoint(0.0235571466151/2,0.8090129379329,0.0954935310336),
QuadraturePoint(0.0235571466151/2,0.0954935310336,0.8090129379329),
QuadraturePoint(0.0235571466151/2,0.0954935310336,0.0954935310336),
QuadraturePoint(0.0268246207430/2,0.2745425238718,0.0479840480721),
QuadraturePoint(0.0268246207430/2,0.0479840480721,0.6774734280561),
QuadraturePoint(0.0268246207430/2,0.6774734280561,0.0479840480721),
QuadraturePoint(0.0268246207430/2,0.6774734280561,0.2745425238718),
QuadraturePoint(0.0268246207430/2,0.2745425238718,0.6774734280561),
QuadraturePoint(0.0268246207430/2,0.0479840480721,0.2745425238718),
QuadraturePoint(0.0314289776779/2,0.4053472446667,0.5429849622344),
QuadraturePoint(0.0314289776779/2,0.0516677930989,0.4053472446667),
QuadraturePoint(0.0314289776779/2,0.4053472446667,0.0516677930989),
QuadraturePoint(0.0314289776779/2,0.5429849622344,0.0516677930989),
QuadraturePoint(0.0314289776779/2,0.0516677930989,0.5429849622344),
QuadraturePoint(0.0314289776779/2,0.5429849622344,0.4053472446667),
QuadraturePoint(0.0337196192159/2,0.1877738615539,0.1068148267588),
QuadraturePoint(0.0337196192159/2,0.7054113116872,0.1877738615539),
QuadraturePoint(0.0337196192159/2,0.7054113116872,0.1068148267588),
QuadraturePoint(0.0337196192159/2,0.1068148267588,0.7054113116872),
QuadraturePoint(0.0337196192159/2,0.1877738615539,0.7054113116872),
QuadraturePoint(0.0337196192159/2,0.1068148267588,0.1877738615539),
QuadraturePoint(0.0427745294213/2,0.1195059712009,0.3057122990643),
QuadraturePoint(0.0427745294213/2,0.1195059712009,0.5747817297348),
QuadraturePoint(0.0427745294213/2,0.5747817297348,0.1195059712009),
QuadraturePoint(0.0427745294213/2,0.5747817297348,0.3057122990643),
QuadraturePoint(0.0427745294213/2,0.3057122990643,0.5747817297348),
QuadraturePoint(0.0427745294213/2,0.3057122990643,0.1195059712009),
QuadraturePoint(0.0441138932737/2,0.5981245743363,0.2009377128319),
QuadraturePoint(0.0441138932737/2,0.2009377128319,0.5981245743363),
QuadraturePoint(0.0441138932737/2,0.2009377128319,0.2009377128319),
QuadraturePoint(0.0461469594684/2,0.2160775200005,0.3121360256673),
QuadraturePoint(0.0461469594684/2,0.3121360256673,0.2160775200005),
QuadraturePoint(0.0461469594684/2,0.2160775200005,0.4717864543321),
QuadraturePoint(0.0461469594684/2,0.3121360256673,0.4717864543321),
QuadraturePoint(0.0461469594684/2,0.4717864543321,0.3121360256673),
QuadraturePoint(0.0461469594684/2,0.4717864543321,0.2160775200005),
QuadraturePoint(0.0469152468624/2,0.4376579903849,0.4376579903849),
QuadraturePoint(0.0469152468624/2,0.4376579903849,0.1246840192303),
QuadraturePoint(0.0469152468624/2,0.1246840192303,0.4376579903849),
QuadraturePoint(0.0551199980347/2,0.3333333333333,0.3333333333333),
};
const QuadratureFormular QuadratureFormular_T_21(21 ,91,P_QuadratureFormular_T_21);


static  QuadraturePoint P_QuadratureFormular_T_23[] = {
QuadraturePoint(0.0006438298261/2,0.0087809303836,0.9903676436772),
QuadraturePoint(0.0006438413076/2,0.9903675314220,0.0087809216232),
QuadraturePoint(0.0010134735710/2,0.0027029276450,0.0335914404439),
QuadraturePoint(0.0010134752576/2,0.0335909214524,0.0027028946710),
QuadraturePoint(0.0019679929935/2,0.0091675068606,0.0091676353051),
QuadraturePoint(0.0033467313784/2,0.9675568182558,0.0084737176656),
QuadraturePoint(0.0033467339208/2,0.0084737200688,0.9675569435345),
QuadraturePoint(0.0042873323375/2,0.0078781948792,0.0676784943862),
QuadraturePoint(0.0042873459885/2,0.0676785477700,0.0078781659291),
QuadraturePoint(0.0043003801372/2,0.9470266955047,0.0442974541187),
QuadraturePoint(0.0043003849098/2,0.0442974755680,0.9470266676487),
QuadraturePoint(0.0056934629205/2,0.9144243214882,0.0081735455132),
QuadraturePoint(0.0056934640134/2,0.0081735424459,0.9144244234031),
QuadraturePoint(0.0061643868015/2,0.2497452292741,0.3833232434720),
QuadraturePoint(0.0061644756418/2,0.3833232646055,0.2497451268005),
QuadraturePoint(0.0062014513591/2,0.8876850353557,0.1035328809446),
QuadraturePoint(0.0062014531952/2,0.1035329228297,0.8876849931840),
QuadraturePoint(0.0069636330294/2,0.0077255923618,0.1403190991974),
QuadraturePoint(0.0069636331842/2,0.1403192425107,0.0077255934624),
QuadraturePoint(0.0075066257720/2,0.8104591009652,0.1809642523926),
QuadraturePoint(0.0075066264565/2,0.1809643003717,0.8104590515334),
QuadraturePoint(0.0079074768339/2,0.8330767948684,0.0083010939677),
QuadraturePoint(0.0079074772485/2,0.0083010907126,0.8330768545392),
QuadraturePoint(0.0080353344623/2,0.0348407706147,0.0348406969482),
QuadraturePoint(0.0087963441074/2,0.2740287679608,0.7173981847948),
QuadraturePoint(0.0087963448112/2,0.7173982224778,0.2740287304386),
QuadraturePoint(0.0091304195716/2,0.2394976858234,0.0081859182262),
QuadraturePoint(0.0091304213611/2,0.0081859185845,0.2394975566677),
QuadraturePoint(0.0092821748751/2,0.0068836152075,0.4843740892687),
QuadraturePoint(0.0092821815662/2,0.4843741485699,0.0068836232949),
QuadraturePoint(0.0094499806178/2,0.4960767772741,0.4960767529507),
QuadraturePoint(0.0094627468484/2,0.6112936776245,0.3804323691239),
QuadraturePoint(0.0094627485294/2,0.3804323980345,0.6112936466533),
QuadraturePoint(0.0095555772285/2,0.7303890713524,0.0083987179701),
QuadraturePoint(0.0095555792843/2,0.0083987168639,0.7303890895407),
QuadraturePoint(0.0096138842488/2,0.6128525675612,0.0075475979695),
QuadraturePoint(0.0096138846826/2,0.0075475961037,0.6128525484582),
QuadraturePoint(0.0099991524212/2,0.0079525316513,0.3559773826721),
QuadraturePoint(0.0099991551850/2,0.3559774870460,0.0079525358502),
QuadraturePoint(0.0100301319277/2,0.9110236977966,0.0437233665345),
QuadraturePoint(0.0100301346636/2,0.0437233605166,0.9110236807446),
QuadraturePoint(0.0124936676185/2,0.0388480061835,0.0967030908282),
QuadraturePoint(0.0124936726125/2,0.0967032117936,0.0388479942386),
QuadraturePoint(0.0140197309137/2,0.0873226911312,0.0873226620391),
QuadraturePoint(0.0143336216896/2,0.0421445202084,0.8485617789108),
QuadraturePoint(0.0143336272125/2,0.8485617974961,0.0421445420915),
QuadraturePoint(0.0153604142740/2,0.8477921333864,0.1067435942472),
QuadraturePoint(0.0153604183425/2,0.1067435889398,0.8477921328146),
QuadraturePoint(0.0184523825614/2,0.1833966521991,0.0416340521608),
QuadraturePoint(0.0184523863146/2,0.0416340541167,0.1833965196930),
QuadraturePoint(0.0195833983573/2,0.7611632251560,0.1941599202852),
QuadraturePoint(0.0195834019994/2,0.1941599254144,0.7611632153938),
QuadraturePoint(0.0197632751342/2,0.7579378747173,0.0439826608586),
QuadraturePoint(0.0197632766677/2,0.0439826512395,0.7579378242308),
QuadraturePoint(0.0198806391019/2,0.0369760535918,0.5363186076436),
QuadraturePoint(0.0198806485776/2,0.5363187134342,0.0369760780935),
QuadraturePoint(0.0207181838484/2,0.1001256948921,0.7912267093545),
QuadraturePoint(0.0207181934893/2,0.7912266693524,0.1001257554673),