quadrature_jacobi_1d.c

一个用来实现偏微分方程中网格的计算库

	    _points[ 3](0)=4.1400214459705974641828724885057e-01L;	_weights[ 3]=6.3146378708891297709533940481199e-02L;
	    _points[ 4](0)=6.0002151327899293019396877288994e-01L;	_weights[ 4]=2.9422211289528614840745351985011e-02L;
	    _points[ 5](0)=7.7351724659143750111199899642603e-01L;	_weights[ 5]=8.1629256323046945072318090687213e-03L;
	    _points[ 6](0)=9.1118316656300276363931736720905e-01L;	_weights[ 6]=8.9268803368920680400902684425561e-04L;

	    return;
	  }
	case FOURTEENTH:
	case FIFTEENTH:
	  {
	    _points.resize (8);
	    _weights.resize(8);
	      
	    _points[ 0](0)=1.6097759551921033532013550096618e-02L;	_weights[ 0]=3.9778957806690559830106455926398e-02L;
	    _points[ 1](0)=8.2900617485651102700366353380388e-02L;	_weights[ 1]=7.6818093267222624804506735068290e-02L;
	    _points[ 2](0)=1.9547516848873991732426660953061e-01L;	_weights[ 2]=8.5284769171938765209441226265186e-02L;
	    _points[ 3](0)=3.4165199147720222027296226216240e-01L;	_weights[ 3]=6.8447183421653272247027200610852e-02L;
	    _points[ 4](0)=5.0559707818448917194006186500608e-01L;	_weights[ 4]=4.0814426388544023186461596152148e-02L;
	    _points[ 5](0)=6.6955227182436145183011451055437e-01L;	_weights[ 5]=1.7246863780234982804243655774862e-02L;
	    _points[ 6](0)=8.1577170358328376075475178697587e-01L;	_weights[ 6]=4.4745217130144079957565210032836e-03L;
	    _points[ 7](0)=9.2850896495990689720101861784923e-01L;	_weights[ 7]=4.6851778403469725578994253231334e-04L;

	    return;
	  }
	case SIXTEENTH:
	case SEVENTEENTH:
	  {
	    _points.resize (9);
	    _weights.resize(9);

	    _points[ 0](0)=1.3165885597114490545190537817972e-02L;	_weights[ 0]=3.2760145111039721461819540084274e-02L;
	    _points[ 1](0)=6.8084529593767587476505856986822e-02L;	_weights[ 1]=6.5489537033396631781939107203912e-02L;
	    _points[ 2](0)=1.6175951676407464206981091232919e-01L;	_weights[ 2]=7.7673569160555504078595170478444e-02L;
	    _points[ 3](0)=2.8589108833922039772798956691264e-01L;	_weights[ 3]=6.9284395689804709215727785348892e-02L;
	    _points[ 4](0)=4.2945364538781279250924800228659e-01L;	_weights[ 4]=4.8540627864506507959528062424865e-02L;
	    _points[ 5](0)=5.7969405635116312626577241331199e-01L;	_weights[ 5]=2.6343280902549835737974696002673e-02L;
	    _points[ 6](0)=7.2326857174033543181796015820021e-01L;	_weights[ 6]=1.0406116579352237052426337235938e-02L;
	    _points[ 7](0)=8.4743684201323732017318042501962e-01L;	_weights[ 7]=2.5743986456087231023317159056037e-03L;
	    _points[ 8](0)=9.4124586421327421141434212713498e-01L;	_weights[ 8]=2.6126234651946294299091864873302e-04L;

	    return;
	  }
	case EIGHTTEENTH:
	case NINTEENTH:
	  {
	    _points.resize (10);
	    _weights.resize(10);

	    _points[ 0](0)=1.0968452456174134250281032990432e-02L;	_weights[ 0]=2.7434088710097496960451648584485e-02L;
	    _points[ 1](0)=5.6898150533657920711434859379121e-02L;	_weights[ 1]=5.6272936402808126282415832152709e-02L;
	    _points[ 2](0)=1.3595023405022895426576342938519e-01L;	_weights[ 2]=7.0069507708666334086372455308183e-02L;
	    _points[ 3](0)=2.4228119613252356111139455536509e-01L;	_weights[ 3]=6.7452219381437824970525794313691e-02L;
	    _points[ 4](0)=3.6800785044933771760476956625017e-01L;	_weights[ 4]=5.2883788766964158209884838942421e-02L;
	    _points[ 5](0)=5.0380712641487390700553064205723e-01L;	_weights[ 5]=3.3854565016814116044101071309373e-02L;
	    _points[ 6](0)=6.3960948865470970932853650616696e-01L;	_weights[ 6]=1.7197575046553008079137502323978e-02L;
	    _points[ 7](0)=7.6534767954810789620953312724460e-01L;	_weights[ 7]=6.4698890685585292297306065609276e-03L;
	    _points[ 8](0)=8.7171007457440848761186134200542e-01L;	_weights[ 8]=1.5455231947365786285853268147186e-03L;
	    _points[ 9](0)=9.5087429264052316644634948461033e-01L;	_weights[ 9]=1.5324003669716084212825702284675e-04L;

	    return;
	  }      
	case TWENTIETH:
	case TWENTYFIRST:
	  {
	    _points.resize (11);
	    _weights.resize(11);

	    _points[ 0](0)=9.2789738313488342029466103013448e-03L;	_weights[ 0]=2.3300850051388377415061404350855e-02L;
	    _points[ 1](0)=4.8249692094286258286483723731132e-02L;	_weights[ 1]=4.8745861030664491864450329028080e-02L;
	    _points[ 2](0)=1.1578862662939522596030591825565e-01L;	_weights[ 2]=6.2979543000367394338333998811129e-02L;
	    _points[ 3](0)=2.0766834159705988244389542530053e-01L;	_weights[ 3]=6.4183559349755944992539371306173e-02L;
	    _points[ 4](0)=3.1811795190623396637744646657914e-01L;	_weights[ 4]=5.4632222321044926070678715192708e-02L;
	    _points[ 5](0)=4.4019839985886246708268132217632e-01L;	_weights[ 5]=3.9290217798435612246535835921329e-02L;
	    _points[ 6](0)=5.6623983915456831356751247193488e-01L;	_weights[ 6]=2.3589663538806928649316661943808e-02L;
	    _points[ 7](0)=6.8832423986295668569553397580809e-01L;	_weights[ 7]=1.1414590918097618249921262806781e-02L;
	    _points[ 8](0)=7.9878435859091464197460872800735e-01L;	_weights[ 8]=4.1400042632190223547676874553245e-03L;
	    _points[ 9](0)=8.9069109935439218915081289990608e-01L;	_weights[ 9]=9.6302057554323908320037373956771e-04L;
	    _points[10](0)=9.5832514378664820192443912466615e-01L;	_weights[10]=9.3800486009778068527692777578977e-05L;

	    return;
	  }
	case TWENTYSECOND:
	case TWENTYTHIRD:
	  {
	    _points.resize (12);
	    _weights.resize(12);

	    _points[ 0](0)=7.9520457026384368762190933400404e-03L;	_weights[ 0]=2.0031112584748752943807724216722e-02L;
	    _points[ 1](0)=4.1427810454294590790770165194795e-02L;	_weights[ 1]=4.2556615367405637772968490900760e-02L;
	    _points[ 2](0)=9.9757625542614157066081703480300e-02L;	_weights[ 2]=5.6584184749474069646478036154467e-02L;
	    _points[ 3](0)=1.7981078905242324954164334430721e-01L;	_weights[ 3]=6.0250607487604518170935343500906e-02L;
	    _points[ 4](0)=2.7727345779931511981753891831052e-01L;	_weights[ 4]=5.4573941164231582825558709780290e-02L;
	    _points[ 5](0)=3.8689200999768979166563470380544e-01L;	_weights[ 5]=4.2764763570188902317265936785879e-02L;
	    _points[ 6](0)=5.0275736044903220327930927843912e-01L;	_weights[ 6]=2.8908051838727360494403492180206e-02L;
	    _points[ 7](0)=6.1862386345845604629354886189114e-01L;	_weights[ 7]=1.6547343644609412101387213540038e-02L;
	    _points[ 8](0)=7.2824645295307268642888033006114e-01L;	_weights[ 8]=7.7163670238864266187877366246530e-03L;
	    _points[ 9](0)=8.2571851421478690176408336295370e-01L;	_weights[ 9]=2.7208421441038393545415361663513e-03L;
	    _points[10](0)=9.0579507354454372680179647039974e-01L;	_weights[10]=6.1995339854451942336673117135933e-04L;
	    _points[11](0)=9.6420653529267155121295530627839e-01L;	_weights[11]=5.9550359808311663832382311702501e-05L;

	    return;
	  }

	case TWENTYFOURTH:
	case TWENTYFIFTH:
	  {
	    _points.resize (13);
	    _weights.resize(13);

	    _points[ 0](0)=6.8908313099587313966815341534070e-03L;	_weights[ 0]=1.7401228610061845488081141733526e-02L;
	    _points[ 1](0)=3.5953362700170095501132811969521e-02L;	_weights[ 1]=3.7426708168844039953666298186474e-02L;
	    _points[ 2](0)=8.6811780157585433284280307690601e-02L;	_weights[ 2]=5.0909830712695019092906847187640e-02L;
	    _points[ 3](0)=1.5709997419650927473052206555148e-01L;	_weights[ 3]=5.6111537130827396185039872203667e-02L;
	    _points[ 4](0)=2.4353289958712844441903410356775e-01L;	_weights[ 4]=5.3340021605871745725401666422880e-02L;
	    _points[ 5](0)=3.4206951264607532347450799656627e-01L;	_weights[ 5]=4.4624600755614954430566337477243e-02L;
	    _points[ 6](0)=4.4810263890570985362873430879398e-01L;	_weights[ 6]=3.2940451678432790083648445397381e-02L;
	    _points[ 7](0)=5.5667462262945501069694664236208e-01L;	_weights[ 7]=2.1254996714756613031113108888658e-02L;
	    _points[ 8](0)=6.6270931356038854232146315212326e-01L;	_weights[ 8]=1.1728259593101880662970381149369e-02L;
	    _points[ 9](0)=7.6124976630827697227665952639343e-01L;	_weights[ 9]=5.3135693388326384480559144318474e-03L;
	    _points[10](0)=8.4769104394723239091031404895142e-01L;	_weights[10]=1.8323911392393785873485644191745e-03L;
	    _points[11](0)=9.1799964553720667953362597028835e-01L;	_weights[11]=4.1072089179088541678428389490865e-04L;
	    _points[12](0)=9.6892889422858896211181181730274e-01L;	_weights[12]=3.9016993264146227750471940564101e-05L;
	
	    return;
	  }

	case TWENTYSIXTH:
	case TWENTYSEVENTH:
	  {
	    _points.resize (14);
	    _weights.resize(14);

	    _points[ 0](0)=6.0288088710666481700887477121600e-03L;	_weights[ 0]=1.5255206119448295373511341972378e-02L;
	    _points[ 1](0)=3.1494259818657590707306712766317e-02L;	_weights[ 1]=3.3139595469960559875930225836517e-02L;
	    _points[ 2](0)=7.6213622474853867586406786340023e-02L;	_weights[ 2]=4.5914772114671178911403799944693e-02L;
	    _points[ 3](0)=1.3836652595947406914513833295825e-01L;	_weights[ 3]=5.2025865410701701832827358849471e-02L;
	    _points[ 4](0)=2.1540975706889557490069553208515e-01L;	_weights[ 4]=5.1389924039471503273808065266413e-02L;
	    _points[ 5](0)=3.0418952814165630963236635490440e-01L;	_weights[ 5]=4.5253432254460881465122373566022e-02L;
	    _points[ 6](0)=4.0107137869715435758118603205199e-01L;	_weights[ 6]=3.5739010292514612984162525810346e-02L;
	    _points[ 7](0)=5.0208913974997841006932539804960e-01L;	_weights[ 7]=2.5216183536839055034306220202219e-02L;
	    _points[ 8](0)=6.0310739672759678368189556654036e-01L;	_weights[ 8]=1.5694332285072957184930434798487e-02L;
	    _points[ 9](0)=6.9999091716990649741073666393849e-01L;	_weights[ 9]=8.4129767371364622096077152778570e-03L;
	    _points[10](0)=7.8877423843456071534030488559891e-01L;	_weights[10]=3.7248756752828588215538124784652e-03L;
	    _points[11](0)=8.6582490073159695367114083518102e-01L;	_weights[11]=1.2617838951574880413283334197260e-03L;
	    _points[12](0)=9.2799573874675232498377812203940e-01L;	_weights[12]=2.7909934094687847008945560658272e-04L;
	    _points[13](0)=9.7277712074118323045296336316727e-01L;	_weights[13]=2.6276161668899854751670304157309e-05L;
	
	    return;
	  }

	case TWENTYEIGHTH:
	case TWENTYNINTH:
	  {
	    _points.resize (15);
	    _weights.resize(15);

	    _points[ 0](0)=5.3190520028437405909978402739774e-03L;	_weights[ 0]=1.3481662725134304090265733542410e-02L;
	    _points[ 1](0)=2.7814561918260051571980709238074e-02L;	_weights[ 1]=2.9527650756798205089388266870014e-02L;
	    _points[ 2](0)=6.7431865838998613384568577571818e-02L;	_weights[ 2]=4.1531307651604265373799456795289e-02L;
	    _points[ 3](0)=1.2274828621549783933958732401303e-01L;	_weights[ 3]=4.8132057668529815888404876719886e-02L;
	    _points[ 4](0)=1.9176570627765362550282684505042e-01L;	_weights[ 4]=4.9041267894409180565196427508457e-02L;
	    _points[ 5](0)=2.7198996592948620305633800793026e-01L;	_weights[ 5]=4.4991555820615933365730152800755e-02L;
	    _points[ 6](0)=3.6052169358705764815778444077840e-01L;	_weights[ 6]=3.7473875792101303095319296790068e-02L;
	    _points[ 7](0)=4.5416123575447709916009624996188e-01L;	_weights[ 7]=2.8326553174241247644068974123766e-02L;
	    _points[ 8](0)=5.4952435559355131365468962196031e-01L;	_weights[ 8]=1.9293416585274520485988288429744e-02L;
	    _points[ 9](0)=6.4316460176878493853949400884187e-01L;	_weights[ 9]=1.1668582484776638563534331197891e-02L;
	    _points[10](0)=7.3169797305252320629933948467521e-01L;	_weights[10]=6.1110501102060490412153142922240e-03L;
	    _points[11](0)=8.1192547727500966668969312120195e-01L;	_weights[11]=2.6556629857743191811282114211333e-03L;
	    _points[12](0)=8.8094952871531360238808055493652e-01L;	_weights[12]=8.8657363523409667398227745346013e-04L;
	    _points[13](0)=9.3628182173551874210251937228715e-01L;	_weights[13]=1.9398726306394554934013378631178e-04L;
	    _points[14](0)=9.7595387433502370956200384127912e-01L;	_weights[14]=1.8128785569508725971591601924058e-05L;
	
	    return;
	  }


	case THIRTIETH:
	case THIRTYFIRST:
	  {
	    _points.resize (16);
	    _weights.resize(16);

	    _points[ 0](0)=4.7276871229345922061506104450645e-03L;	_weights[ 0]=1.1999405624691969586506695404143e-02L;
	    _points[ 1](0)=2.4742967619434984512603456269536e-02L;	_weights[ 1]=2.6460842517339527040437264564557e-02L;
	    _points[ 2](0)=6.0076437716637992275242255328112e-02L;	_weights[ 2]=3.7685824681918572833315081444311e-02L;
	    _points[ 3](0)=1.0960060925877060286789758461950e-01L;	_weights[ 3]=4.4496573472543214976151049587710e-02L;
	    _points[ 4](0)=1.7172475948300210790307417791270e-01L;	_weights[ 4]=4.6506239397362658797012521100333e-02L;
	    _points[ 5](0)=2.4445243174451137730207656160275e-01L;	_weights[ 5]=4.4112657675899371113911374089320e-02L;
	    _points[ 6](0)=3.2544621002412961482683685190181e-01L;	_weights[ 6]=3.8348046951185955673582690910292e-02L;
	    _points[ 7](0)=4.1210296411685959691728421362251e-01L;	_weights[ 7]=3.0612037917837609541948834407964e-02L;
	    _points[ 8](0)=5.0163755615318759370013043235084e-01L;	_weights[ 8]=2.2357185649529930223334734263644e-02L;
	    _points[ 9](0)=5.9117238531650308897076051877634e-01L;	_weights[ 9]=1.4803776657751738340239250997704e-02L;
	    _points[10](0)=6.7782991741456081130251277264148e-01L;	_weights[10]=8.7475454956120699745801506391282e-03L;
	    _points[11](0)=7.5882525794914141335427940535114e-01L;	_weights[11]=4.4948488514798588832066397608164e-03L;
	    _points[12](0)=8.3155588302646195188920635482406e-01L;	_weights[12]=1.9235319260667742273824398535627e-03L;
	    _points[13](0)=8.9368597454880401881486969876816e-01L;	_weights[13]=6.3445999956883686801084896849641e-04L;
	    _points[14](0)=9.4322428571224559352545066239296e-01L;	_weights[14]=1.3757810397555175246959527906570e-04L;
	    _points[15](0)=9.7860643749869701257280091378128e-01L;	_weights[15]=1.2778410569693501244162062287056e-05L;
	
	    return;
	  }


	case THIRTYSECOND:
	case THIRTYTHIRD: