qf11to25.cpp

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

QuadraturePoint(0.0208943071440/2,0.0379866714177,0.4157413128558),
QuadraturePoint(0.0208943251956/2,0.4157414028965,0.0379867061535),
QuadraturePoint(0.0214864573885/2,0.6507106491463,0.0420141226713),
QuadraturePoint(0.0214864586007/2,0.0420141133438,0.6507105645084),
QuadraturePoint(0.0222218133036/2,0.0425548444254,0.2920626023484),
QuadraturePoint(0.0222218160203/2,0.2920627107240,0.0425548546753),
QuadraturePoint(0.0223345305455/2,0.5389729538180,0.4193031469005),
QuadraturePoint(0.0223345378739/2,0.4193031828489,0.5389729093610),
QuadraturePoint(0.0224758924946/2,0.6549472009700,0.3007352636162),
QuadraturePoint(0.0224758980440/2,0.3007352790917,0.6549471812731),
QuadraturePoint(0.0229701395845/2,0.3752400771585,0.3453980130752),
QuadraturePoint(0.0229703394438/2,0.3453980282786,0.3752400695673),
QuadraturePoint(0.0232798376102/2,0.0994532168761,0.1598308695187),
QuadraturePoint(0.0232798427506/2,0.1598309359585,0.0994531960132),
QuadraturePoint(0.0269483199647/2,0.1797326661667,0.7124585430924),
QuadraturePoint(0.0269483307107/2,0.7124584461943,0.1797327722240),
QuadraturePoint(0.0280438758010/2,0.1066065678636,0.7001701784175),
QuadraturePoint(0.0280438764607/2,0.7001701904096,0.1066065855677),
QuadraturePoint(0.0287526270172/2,0.0993303629801,0.6065647984796),
QuadraturePoint(0.0287526387271/2,0.6065648052521,0.0993303896769),
QuadraturePoint(0.0298980829063/2,0.1023223542704,0.2533381579528),
QuadraturePoint(0.0298980922759/2,0.2533382324938,0.1023223826189),
QuadraturePoint(0.0309004358516/2,0.6166226715217,0.2769502060575),
QuadraturePoint(0.0309004385956/2,0.2769500693109,0.6166227900624),
QuadraturePoint(0.0314031017088/2,0.0904184571873,0.4981522637001),
QuadraturePoint(0.0314031073955/2,0.4981522767248,0.0904185045149),
QuadraturePoint(0.0319191553024/2,0.0928231860168,0.3738418516908),
QuadraturePoint(0.0319191668378/2,0.3738418699229,0.0928232584790),
QuadraturePoint(0.0321429924062/2,0.2521678840407,0.2521680925697),
QuadraturePoint(0.0330395601388/2,0.5087500218708,0.3905580544330),
QuadraturePoint(0.0330395631829/2,0.3905579116731,0.5087501437661),
QuadraturePoint(0.0356169095589/2,0.1706141469096,0.5266738039554),
QuadraturePoint(0.0356169276054/2,0.5266737761312,0.1706142257537),
QuadraturePoint(0.0365741189998/2,0.3487581527629,0.2588055084886),
QuadraturePoint(0.0365741515204/2,0.2588053596017,0.3487583491703),
QuadraturePoint(0.0365977646990/2,0.1696614558053,0.3013522183964),
QuadraturePoint(0.0365978053889/2,0.3013521806875,0.1696615963219),
QuadraturePoint(0.0369945680114/2,0.2580202409759,0.4584741774478),
QuadraturePoint(0.0369945775059/2,0.4584740860198,0.2580203819011),
QuadraturePoint(0.0374053623787/2,0.1848898683498,0.1848898704551),
QuadraturePoint(0.0375550258317/2,0.6130740338465,0.1921611994069),
QuadraturePoint(0.0375550312530/2,0.1921611750994,0.6130740398389),
QuadraturePoint(0.0388887693486/2,0.4180541160599,0.1650613336416),
QuadraturePoint(0.0388887708342/2,0.1650612642036,0.4180541199244),
QuadraturePoint(0.0392705643548/2,0.5159205739625,0.2982719005229),
QuadraturePoint(0.0392705802517/2,0.2982718935750,0.5159205534362),
QuadraturePoint(0.0398766879831/2,0.4098894602340,0.4098894317792),
};
const QuadratureFormular QuadratureFormular_T_23(23 ,105,P_QuadratureFormular_T_23);


static  QuadraturePoint P_QuadratureFormular_T_25[] = {
QuadraturePoint(0.0014873417859/2,0.0082881595033,0.9848202768869),
QuadraturePoint(0.0014889035262/2,0.4618422030241,0.5381577969759),
QuadraturePoint(0.0015005944380/2,0.0071066441239,0.0080842361390),
QuadraturePoint(0.0015059208313/2,0.9847613141699,0.0070015755134),
QuadraturePoint(0.0015318868715/2,0.5374447869049,0.4625552130951),
QuadraturePoint(0.0023032634487/2,0.0000000000000,0.4887676880140),
QuadraturePoint(0.0023649067042/2,0.4914131929361,0.0000000000000),
QuadraturePoint(0.0028751143611/2,0.0070345937020,0.9574158053697),
QuadraturePoint(0.0029862488735/2,0.9564734714228,0.0364655449485),
QuadraturePoint(0.0030384162737/2,0.0370198792045,0.0070908577166),
QuadraturePoint(0.0032092459688/2,0.1024124542747,0.8936125594937),
QuadraturePoint(0.0037029598435/2,0.5928065811509,0.0049451705600),
QuadraturePoint(0.0037407186035/2,0.0050948422371,0.0996676659189),
QuadraturePoint(0.0038452543223/2,0.0081562023689,0.0415561148784),
QuadraturePoint(0.0038670778668/2,0.0424936107568,0.9494865260352),
QuadraturePoint(0.0039192555178/2,0.9495543500844,0.0081794507292),
QuadraturePoint(0.0039573282688/2,0.8932787471239,0.0053224326262),
QuadraturePoint(0.0044032251724/2,0.0069317612927,0.9065401020433),
QuadraturePoint(0.0045907108173/2,0.9035839030665,0.0894771171077),
QuadraturePoint(0.0047023669435/2,0.0905665738209,0.0070525342005),
QuadraturePoint(0.0050014843818/2,0.0083929332787,0.6663179931111),
QuadraturePoint(0.0052387830156/2,0.6261245686071,0.0092197583153),
QuadraturePoint(0.0054422104092/2,0.0062801592979,0.8335207460527),
QuadraturePoint(0.0056931248912/2,0.8272539257367,0.1665134939330),
QuadraturePoint(0.0059107422989/2,0.0062005875353,0.7424693255229),
QuadraturePoint(0.0059687967687/2,0.1676900311185,0.0065717743528),
QuadraturePoint(0.0067262190287/2,0.7199353069567,0.0064354534962),
QuadraturePoint(0.0068307848624/2,0.2749740090237,0.7185296120719),
QuadraturePoint(0.0069531259112/2,0.0079257582005,0.1766411374714),
QuadraturePoint(0.0072460270642/2,0.0069981220752,0.2704767254004),
QuadraturePoint(0.0072728189613/2,0.8125248773263,0.0082299533210),
QuadraturePoint(0.0073008930847/2,0.0073536969970,0.5934167875453),
QuadraturePoint(0.0073604666776/2,0.7283665935411,0.2648817553752),
QuadraturePoint(0.0074119923255/2,0.1800642304565,0.8115848976682),
QuadraturePoint(0.0074892214336/2,0.2658102467762,0.0068553525429),
QuadraturePoint(0.0078604067260/2,0.0070892364520,0.3757632659744),
QuadraturePoint(0.0078621726423/2,0.3774054302043,0.6148573533757),
QuadraturePoint(0.0080506361066/2,0.0369649608668,0.9210792302893),
QuadraturePoint(0.0081442860473/2,0.9203194109805,0.0426025082114),
QuadraturePoint(0.0081478804152/2,0.0425477806431,0.0372689941794),
QuadraturePoint(0.0092444146612/2,0.6191278394983,0.3724055713809),
QuadraturePoint(0.0094674635165/2,0.3762697209178,0.0081436422011),
QuadraturePoint(0.0097132210137/2,0.0956111149690,0.8771098372601),
QuadraturePoint(0.0099753581151/2,0.0302473410377,0.0943858903393),
QuadraturePoint(0.0103367803673/2,0.8739905691754,0.0313198990883),
QuadraturePoint(0.0112263277166/2,0.8604133734958,0.1049019782046),
QuadraturePoint(0.0114309118745/2,0.0347307852352,0.8609856462886),
QuadraturePoint(0.0115550567487/2,0.1043606608343,0.0357152881004),
QuadraturePoint(0.0135575856957/2,0.7797622824754,0.1872318199265),
QuadraturePoint(0.0135984962900/2,0.0185865164256,0.4834397678794),
QuadraturePoint(0.0137754813837/2,0.0324585286618,0.7783474916042),
QuadraturePoint(0.0137961015942/2,0.8371293901157,0.0804060570156),
QuadraturePoint(0.0138408839904/2,0.0836602075315,0.8421414817051),
QuadraturePoint(0.0140634019977/2,0.0784070242501,0.0849927089145),
QuadraturePoint(0.0140991451009/2,0.4929238648458,0.4892855914710),
QuadraturePoint(0.0142004111991/2,0.1870637584073,0.0345210858281),
QuadraturePoint(0.0144518424517/2,0.4892636967025,0.0190774755077),
QuadraturePoint(0.0150245979639/2,0.0401982618372,0.1691143187109),
QuadraturePoint(0.0152817804122/2,0.7894259278865,0.0412206731484),
QuadraturePoint(0.0155550724169/2,0.1686260456429,0.7894860640585),
QuadraturePoint(0.0164570886000/2,0.3750901913174,0.5895318272013),
QuadraturePoint(0.0165275759573/2,0.0356362876880,0.3681256217699),
QuadraturePoint(0.0166847554451/2,0.5887548164804,0.0359968962541),
QuadraturePoint(0.0167409312985/2,0.0373308082182,0.6790704673533),
QuadraturePoint(0.0168674663361/2,0.2820769993374,0.0373639992361),
QuadraturePoint(0.0168882230165/2,0.6819277603320,0.2803330345725),
QuadraturePoint(0.0172087112691/2,0.0374938324382,0.2634016180014),
QuadraturePoint(0.0174681068264/2,0.6984079204127,0.0364154673322),
QuadraturePoint(0.0176663899614/2,0.2654390894079,0.6980717436193),
QuadraturePoint(0.0182967621475/2,0.1429848440800,0.7612254618453),
QuadraturePoint(0.0183576852459/2,0.7623554007647,0.0943741220275),
QuadraturePoint(0.0186392569521/2,0.0934222022749,0.1479799836832),
QuadraturePoint(0.0189781060590/2,0.5759004479923,0.3821329641698),
QuadraturePoint(0.0191847922578/2,0.3822427332525,0.0426716362301),
QuadraturePoint(0.0194080442044/2,0.0411414081675,0.5718082874432),
QuadraturePoint(0.0194720072193/2,0.0802462538379,0.7702204382042),
QuadraturePoint(0.0200855080495/2,0.7625229819410,0.1559420577362),
QuadraturePoint(0.0201673909332/2,0.1524941445131,0.0842965421322),
QuadraturePoint(0.0221742162761/2,0.0622159195833,0.4538181318873),
QuadraturePoint(0.0229702440508/2,0.1109539036076,0.4586014071171),
QuadraturePoint(0.0233465117399/2,0.4575627212057,0.4795313560210),
QuadraturePoint(0.0234883135338/2,0.4322865136374,0.1230591237472),
QuadraturePoint(0.0240682099018/2,0.5865002850241,0.0834119779793),
QuadraturePoint(0.0240910792953/2,0.0869359250818,0.6755677013351),
QuadraturePoint(0.0245677049481/2,0.0929594906936,0.2326500892727),
QuadraturePoint(0.0246536315719/2,0.6661932141454,0.2448294007406),
QuadraturePoint(0.0246756530052/2,0.4780306362227,0.0661749044835),
QuadraturePoint(0.0249704602710/2,0.4372215294577,0.4442145585244),
QuadraturePoint(0.0250026544082/2,0.6779224504669,0.0929096534577),
QuadraturePoint(0.0250490869426/2,0.2423431255660,0.0889793655129),
QuadraturePoint(0.0250936250125/2,0.2288925420305,0.6780053081672),
QuadraturePoint(0.0251482076226/2,0.3315065049959,0.5847381559741),
QuadraturePoint(0.0255010290447/2,0.3424200526607,0.5139245722736),
QuadraturePoint(0.0256544511979/2,0.0862630046475,0.3340976249234),
QuadraturePoint(0.0257974750630/2,0.5113188946635,0.1380154720554),
QuadraturePoint(0.0270007753993/2,0.1538977841001,0.6788062619562),
QuadraturePoint(0.0274431536844/2,0.6779951348472,0.1663358925269),
QuadraturePoint(0.0277072401488/2,0.1664600469411,0.1582214504849),
QuadraturePoint(0.0278284415364/2,0.0950910318888,0.5666590332543),
QuadraturePoint(0.0287207381105/2,0.3436048136712,0.0978960873457),
QuadraturePoint(0.0288826834956/2,0.5560417025366,0.3468917820947),
QuadraturePoint(0.0293302729759/2,0.1452404029513,0.3599534491052),
QuadraturePoint(0.0318902879557/2,0.1619685156238,0.5810131373330),
QuadraturePoint(0.0319083660286/2,0.5800164844262,0.2560674640672),
QuadraturePoint(0.0320938960329/2,0.2450201223288,0.5881469552102),
QuadraturePoint(0.0321618608780/2,0.2557621891794,0.1652244065047),
QuadraturePoint(0.0322424127534/2,0.2205239985511,0.3496507466106),
QuadraturePoint(0.0327072446421/2,0.4940183111285,0.2549448448453),
QuadraturePoint(0.0329946316695/2,0.2531570689798,0.2543369115017),
QuadraturePoint(0.0331828096025/2,0.5846891116357,0.1666603916479),
QuadraturePoint(0.0334857162651/2,0.1660333602278,0.2523240191705),
QuadraturePoint(0.0335468472792/2,0.2505426292461,0.4959007627528),
QuadraturePoint(0.0337049042988/2,0.3519336802182,0.1805380367800),
QuadraturePoint(0.0340361462767/2,0.3502668835419,0.4358582329881),
QuadraturePoint(0.0342465235323/2,0.4400892485512,0.2120576104941),
QuadraturePoint(0.0345528817251/2,0.4680855471546,0.3552681570774),
QuadraturePoint(0.0356782875703/2,0.1770237763947,0.4670352922266),
QuadraturePoint(0.0364656225016/2,0.3900920779501,0.3323152819300),
QuadraturePoint(0.0365172708706/2,0.2805847774120,0.3898041176680),
QuadraturePoint(0.0371924811018/2,0.3361523347440,0.2778500044356),
};
const QuadratureFormular QuadratureFormular_T_25(25 ,120,P_QuadratureFormular_T_25);

//  Hack to do something at initialisation time
//   to add the name myfunction to the freefem++ table 
class Init { public:
  Init();
};
Init init;
Init::Init(){
  /*
grep QuadratureFormular QF.cpp|grep ^const|awk -F"[_(]" '{print "Global.New(@qf" $3 "pT@,CConstant<const QuadratureFormular *>(&QuadratureFormular_T_"$3"));"}'|sed -e 's/@/"/g'
   */
	     Global.New("qf11pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_11));
	     Global.New("qf13pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_13));
	     Global.New("qf14pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_14));
	     Global.New("qf16pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_16));
	     Global.New("qf18pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_18));
	     Global.New("qf20pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_20));
	     Global.New("qf21pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_21));
	     Global.New("qf23pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_23));
	     Global.New("qf25pT",CConstant<const QuadratureFormular *>(&QuadratureFormular_T_25));
	     Global.New("qf11pE",CConstant<const QuadratureFormular1d *>(new QuadratureFormular1d(-1+2*11,11,GaussLegendre(11),true)));
	     Global.New("qf12pE",CConstant<const QuadratureFormular1d *>(new QuadratureFormular1d(-1+2*12,12,GaussLegendre(12),true)));
	     Global.New("qf13pE",CConstant<const QuadratureFormular1d *>(new QuadratureFormular1d(-1+2*13,13,GaussLegendre(13),true)));
}