qf11to25.cpp

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

QuadraturePoint(0.0316799866332/2,0.6647637544849,0.0186886898773),
QuadraturePoint(0.0316799866332/2,0.0186886898773,0.6647637544849),
QuadraturePoint(0.0316799866332/2,0.0186886898773,0.3165475556378),
QuadraturePoint(0.0316799866332/2,0.3165475556378,0.6647637544849),
QuadraturePoint(0.0316799866332/2,0.6647637544849,0.3165475556378),
QuadraturePoint(0.0349317947036/2,0.0192662192492,0.4903668903754),
QuadraturePoint(0.0349317947036/2,0.4903668903754,0.0192662192492),
QuadraturePoint(0.0349317947036/2,0.4903668903754,0.4903668903754),
QuadraturePoint(0.0383664533945/2,0.0875134669581,0.8249730660837),
QuadraturePoint(0.0383664533945/2,0.0875134669581,0.0875134669581),
QuadraturePoint(0.0383664533945/2,0.8249730660837,0.0875134669581),
QuadraturePoint(0.0578369491210/2,0.0935526036219,0.2079865423167),
QuadraturePoint(0.0578369491210/2,0.0935526036219,0.6984608540613),
QuadraturePoint(0.0578369491210/2,0.2079865423167,0.0935526036219),
QuadraturePoint(0.0578369491210/2,0.6984608540613,0.0935526036219),
QuadraturePoint(0.0578369491210/2,0.6984608540613,0.2079865423167),
QuadraturePoint(0.0578369491210/2,0.2079865423167,0.6984608540613),
QuadraturePoint(0.0725821687394/2,0.0974892983467,0.5380088595149),
QuadraturePoint(0.0725821687394/2,0.3645018421383,0.0974892983467),
QuadraturePoint(0.0725821687394/2,0.5380088595149,0.0974892983467),
QuadraturePoint(0.0725821687394/2,0.5380088595149,0.3645018421383),
QuadraturePoint(0.0725821687394/2,0.3645018421383,0.5380088595149),
QuadraturePoint(0.0725821687394/2,0.0974892983467,0.3645018421383),
QuadraturePoint(0.0897856524107/2,0.2217145894873,0.5565708210253),
QuadraturePoint(0.0897856524107/2,0.5565708210253,0.2217145894873),
QuadraturePoint(0.0897856524107/2,0.2217145894873,0.2217145894873),
QuadraturePoint(0.1034544533617/2,0.3860471669296,0.2279056661408),
QuadraturePoint(0.1034544533617/2,0.2279056661408,0.3860471669296),
QuadraturePoint(0.1034544533617/2,0.3860471669296,0.3860471669296),
};
const QuadratureFormular QuadratureFormular_T_14(14 ,45,P_QuadratureFormular_T_14);


static  QuadraturePoint P_QuadratureFormular_T_16[] = {
QuadraturePoint(0.0006202599851/2,1.0000000000000,0.0000000000000),
QuadraturePoint(0.0006315174712/2,0.0000000000000,1.0000000000000),
QuadraturePoint(0.0007086601559/2,0.0000000000000,0.0000000000000),
QuadraturePoint(0.0055163716168/2,0.9398863583577,0.0049848744634),
QuadraturePoint(0.0062692407656/2,0.0543806683058,0.9386405618617),
QuadraturePoint(0.0078531408826/2,0.0093940049164,0.0526424462697),
QuadraturePoint(0.0094551483864/2,0.0164345086362,0.9469035517351),
QuadraturePoint(0.0097824511271/2,0.9469487269862,0.0363373677167),
QuadraturePoint(0.0099861643489/2,0.0426604005768,0.0151224541799),
QuadraturePoint(0.0137553818816/2,0.0122269495439,0.8693773510664),
QuadraturePoint(0.0140979178040/2,0.8673696521047,0.1204917285774),
QuadraturePoint(0.0149646864337/2,0.8456744021389,0.0157763967870),
QuadraturePoint(0.0156097503612/2,0.1395759632103,0.8448120870375),
QuadraturePoint(0.0157683693348/2,0.1317821743231,0.0135009605584),
QuadraturePoint(0.0175794546383/2,0.0157955126300,0.1455274938536),
QuadraturePoint(0.0204113840270/2,0.7365462884436,0.0155697540908),
QuadraturePoint(0.0209562878616/2,0.0139688430330,0.7379836894450),
QuadraturePoint(0.0210713412998/2,0.2547895186039,0.7297615689771),
QuadraturePoint(0.0217646760202/2,0.7316386522555,0.2543076683315),
QuadraturePoint(0.0222288408699/2,0.0157253728951,0.2696239795791),
QuadraturePoint(0.0224186693682/2,0.2662302843647,0.0144783956308),
QuadraturePoint(0.0230122616993/2,0.8673504065214,0.0591679410400),
QuadraturePoint(0.0236813902500/2,0.0741493666957,0.8634782575061),
QuadraturePoint(0.0257464643368/2,0.0159285948360,0.4191238955238),
QuadraturePoint(0.0257956801608/2,0.0156061028068,0.5809222921146),
QuadraturePoint(0.0258072327610/2,0.5910094817484,0.0159251452651),
QuadraturePoint(0.0260343232059/2,0.4034771496889,0.5806700368104),
QuadraturePoint(0.0265768141609/2,0.5694745628526,0.4149495146302),
QuadraturePoint(0.0265784761831/2,0.0678493700650,0.0761218678591),
QuadraturePoint(0.0267532329238/2,0.4265968590272,0.0157509692312),
QuadraturePoint(0.0375787806641/2,0.0670982507890,0.7741898312421),
QuadraturePoint(0.0383065894195/2,0.7528310231480,0.0819119495639),
QuadraturePoint(0.0384849695025/2,0.7753727783557,0.1577128457292),
QuadraturePoint(0.0389619825852/2,0.1689073157787,0.7503943099742),
QuadraturePoint(0.0394604111547/2,0.1687335832919,0.0708311507268),
QuadraturePoint(0.0412364778098/2,0.0821244708436,0.1762996626771),
QuadraturePoint(0.0512872438483/2,0.6288705363345,0.0807744953317),
QuadraturePoint(0.0516405641935/2,0.0811413015266,0.3054373589776),
QuadraturePoint(0.0518230042269/2,0.2969112065080,0.6227485988871),
QuadraturePoint(0.0528527988181/2,0.0767542314171,0.6247247149546),
QuadraturePoint(0.0538505573027/2,0.6223022333845,0.3011485821166),
QuadraturePoint(0.0541895329319/2,0.3103786288051,0.0779098365079),
QuadraturePoint(0.0584737146444/2,0.0819218215187,0.4603633038351),
QuadraturePoint(0.0592863168363/2,0.4717022665013,0.0821554006797),
QuadraturePoint(0.0594358276749/2,0.4546603415250,0.4637565033890),
QuadraturePoint(0.0631800255863/2,0.1701091339237,0.6422277808188),
QuadraturePoint(0.0632926845153/2,0.6406004329487,0.1898293537256),
QuadraturePoint(0.0640707361772/2,0.1912267583717,0.1739955685343),
QuadraturePoint(0.0812040595918/2,0.1885315767070,0.4798914070406),
QuadraturePoint(0.0814437513530/2,0.4772929957691,0.3348356598119),
QuadraturePoint(0.0814679201241/2,0.3126974621760,0.4957972197259),
QuadraturePoint(0.0815050548084/2,0.4961225945946,0.1927553668904),
QuadraturePoint(0.0815164664939/2,0.1928805312867,0.3161015807261),
QuadraturePoint(0.0816931059623/2,0.3360041453816,0.1894892801290),
QuadraturePoint(0.0923218334531/2,0.3337280550848,0.3343571021811),
};
const QuadratureFormular QuadratureFormular_T_16(16 ,55,P_QuadratureFormular_T_16);


static QuadraturePoint P_QuadratureFormular_T_18[] = {
QuadraturePoint(0.0025165756986/2,0.0116731059668,0.9812565951289),
QuadraturePoint(0.0025273452007/2,0.9810030858388,0.0071462504863),
QuadraturePoint(0.0033269295333/2,0.0106966317092,0.0115153933376),
QuadraturePoint(0.0081503492125/2,0.9382476983551,0.0495570591341),
QuadraturePoint(0.0086135525742/2,0.0126627518417,0.9370123620615),
QuadraturePoint(0.0087786746179/2,0.0598109409984,0.0121364578922),
QuadraturePoint(0.0097099585562/2,0.0137363297927,0.0612783625597),
QuadraturePoint(0.0102466211915/2,0.9229527959405,0.0141128270602),
QuadraturePoint(0.0108397688341/2,0.0633107354993,0.9220197291727),
QuadraturePoint(0.0129385390176/2,0.0117265100335,0.1500520475229),
QuadraturePoint(0.0136339823583/2,0.1554720587323,0.8325147121589),
QuadraturePoint(0.0138477328147/2,0.8343293888982,0.0125228158759),
QuadraturePoint(0.0139421540105/2,0.8501638031957,0.1371997508736),
QuadraturePoint(0.0144121399968/2,0.0128816350522,0.8477627063479),
QuadraturePoint(0.0153703455534/2,0.1510801608959,0.0136526924039),
QuadraturePoint(0.0162489802253/2,0.0101917879217,0.5770438618345),
QuadraturePoint(0.0169718304280/2,0.2813372399303,0.7066853759623),
QuadraturePoint(0.0170088532421/2,0.7124374628501,0.0124569780990),
QuadraturePoint(0.0170953520675/2,0.2763025250863,0.0121741311386),
QuadraturePoint(0.0173888854559/2,0.0109658368561,0.4194306712466),
QuadraturePoint(0.0174543962439/2,0.4289110517884,0.5599616067469),
QuadraturePoint(0.0178406757287/2,0.4215420555115,0.0116475994785),
QuadraturePoint(0.0178446863879/2,0.5711258590444,0.0118218313989),
QuadraturePoint(0.0179046337552/2,0.5826868270511,0.4057889581177),
QuadraturePoint(0.0181259756201/2,0.0130567806713,0.2725023750868),
QuadraturePoint(0.0184784838882/2,0.0130760400964,0.7224712523233),
QuadraturePoint(0.0185793564371/2,0.7263437062407,0.2602984019251),
QuadraturePoint(0.0203217151777/2,0.0687230068637,0.0631417277210),
QuadraturePoint(0.0213771661809/2,0.8652302101529,0.0720611837338),
QuadraturePoint(0.0231916854098/2,0.0648599071037,0.8590433543910),
QuadraturePoint(0.0274426710859/2,0.1483494943362,0.7888788352240),
QuadraturePoint(0.0290301922340/2,0.0624359898396,0.1493935499354),
QuadraturePoint(0.0294522738505/2,0.7871369011735,0.0656382042757),
QuadraturePoint(0.0299436251629/2,0.0519104921610,0.5255635695605),
QuadraturePoint(0.0307026948119/2,0.1543129927444,0.0716383926917),
QuadraturePoint(0.0325263365863/2,0.2617842745603,0.0621479485288),
QuadraturePoint(0.0327884208506/2,0.7667257872813,0.1658211554831),
QuadraturePoint(0.0331234675192/2,0.2582103676627,0.6800119766139),
QuadraturePoint(0.0346167526875/2,0.0679065925147,0.7571515437782),
QuadraturePoint(0.0347081373976/2,0.5293578274804,0.4121503841107),
QuadraturePoint(0.0347372049404/2,0.0666036150484,0.2612513087886),
QuadraturePoint(0.0348528762454/2,0.0585675461899,0.3902236114535),
QuadraturePoint(0.0348601561186/2,0.0644535360411,0.6373626559761),
QuadraturePoint(0.0355471569975/2,0.6748138429151,0.0637583342061),
QuadraturePoint(0.0360182996383/2,0.3914602310369,0.5503238090563),
QuadraturePoint(0.0362926285843/2,0.6487701492307,0.2836728360263),
QuadraturePoint(0.0381897702083/2,0.3946498220408,0.0605175522554),
QuadraturePoint(0.0392252800118/2,0.5390137151933,0.0611990176936),
QuadraturePoint(0.0482710125888/2,0.1627895082785,0.6861322141035),
QuadraturePoint(0.0489912121566/2,0.6812436322641,0.1567968345899),
QuadraturePoint(0.0497220833872/2,0.1542832878020,0.1667512624020),
QuadraturePoint(0.0507065736986/2,0.2522727750445,0.2504803933395),
QuadraturePoint(0.0509771994043/2,0.2547981532407,0.4994090649043),
QuadraturePoint(0.0521360063667/2,0.1485580549194,0.5756023096087),
QuadraturePoint(0.0523460874925/2,0.2930239606436,0.5656897354162),
QuadraturePoint(0.0524440683552/2,0.2808991272310,0.1437921574248),
QuadraturePoint(0.0527459644823/2,0.4820989592971,0.2518557535865),
QuadraturePoint(0.0529449063728/2,0.5641878245444,0.1462966743153),
QuadraturePoint(0.0542395594501/2,0.1307699644344,0.4489577586117),
QuadraturePoint(0.0543470203419/2,0.1479692221948,0.3001174386829),
QuadraturePoint(0.0547100548639/2,0.5638684222946,0.2813772089298),
QuadraturePoint(0.0557288345913/2,0.4361157428790,0.4252053446420),
QuadraturePoint(0.0577734264233/2,0.3603263935285,0.2599190004889),
QuadraturePoint(0.0585393781623/2,0.4224188334674,0.1453238443303),
QuadraturePoint(0.0609039250680/2,0.3719001833052,0.3780122703567),
QuadraturePoint(0.0637273964449/2,0.2413645006928,0.3847563284940),
};
const QuadratureFormular QuadratureFormular_T_18(18 ,66,P_QuadratureFormular_T_18);


static  QuadraturePoint P_QuadratureFormular_T_20[] = {
QuadraturePoint(0.0021744545399/2,0.0089411337112,0.0086983293702),
QuadraturePoint(0.0028987135265/2,0.9792622629807,0.0102644133744),
QuadraturePoint(0.0030846029337/2,0.0105475382112,0.9785514202515),
QuadraturePoint(0.0034401633104/2,0.0023777061947,0.0636551098604),
QuadraturePoint(0.0041898472012/2,0.0630425115795,0.0041506347509),
QuadraturePoint(0.0044738051498/2,0.9308422496730,0.0048053482263),
QuadraturePoint(0.0047054420814/2,0.0629076555490,0.9316790069481),
QuadraturePoint(0.0048867935750/2,0.9315962246381,0.0626264881801),
QuadraturePoint(0.0051927643369/2,0.0061951689415,0.9293587058564),
QuadraturePoint(0.0074073058981/2,0.0287125819237,0.0310202122997),
QuadraturePoint(0.0079755410301/2,0.9293844478305,0.0342152968219),
QuadraturePoint(0.0083550522910/2,0.0375457566621,0.9257868884669),
QuadraturePoint(0.0096166660864/2,0.0086895739064,0.1584971251510),
QuadraturePoint(0.0096318257850/2,0.1547597053965,0.8363606657688),
QuadraturePoint(0.0098577460758/2,0.8331025294185,0.0089257244824),
QuadraturePoint(0.0102657880301/2,0.8374231073526,0.1529167304078),
QuadraturePoint(0.0103188103111/2,0.1559362505234,0.0094966240058),
QuadraturePoint(0.0106291001630/2,0.0098599642095,0.8342211493596),
QuadraturePoint(0.0106881306895/2,0.4055873733289,0.0074389302008),
QuadraturePoint(0.0106969021010/2,0.5964727898618,0.3956330809311),
QuadraturePoint(0.0109026461714/2,0.0080747800416,0.4031319425903),
QuadraturePoint(0.0109899783575/2,0.0075073977721,0.5851609594681),
QuadraturePoint(0.0113423055229/2,0.3936764519237,0.5974896592899),
QuadraturePoint(0.0120535642930/2,0.5846530726212,0.0087250464968),
QuadraturePoint(0.0139619193821/2,0.4870804112120,0.0202129229912),
QuadraturePoint(0.0141147991536/2,0.2683512811785,0.7202340088668),
QuadraturePoint(0.0141930347046/2,0.7223956288748,0.2662399366456),
QuadraturePoint(0.0144212676268/2,0.2716826742357,0.0112882698808),
QuadraturePoint(0.0144704346855/2,0.0112580842046,0.7169695963325),
QuadraturePoint(0.0144949769872/2,0.0115034734370,0.2740067110166),
QuadraturePoint(0.0145386775694/2,0.7140525900564,0.0113511560497),
QuadraturePoint(0.0145964190926/2,0.4902871053112,0.4936491841468),
QuadraturePoint(0.0147314578466/2,0.0201423425209,0.4832573459601),
QuadraturePoint(0.0167463963304/2,0.0361107464859,0.0935679501582),
QuadraturePoint(0.0168955500458/2,0.8607998819851,0.0397379067075),
QuadraturePoint(0.0169422662884/2,0.1005891526001,0.8586343419352),