test_nuderiv4d.adb running pde44e_nuderiv4d.adb running uxxxx(x,y,z,t) + uyyyy(x,y,z,t) + uzzzz(x,y,z,t) + utttt(x,y,z,t) + u(x,y,z,t) = F(x,y,z,t) F(x,y,z,t) = 2*(exp(x)+exp(y)+exp(z)+exp(t)) uana(x,y)=exp(x)+exp(y)+exp(z)+exp(t) solution vertex, x, y, z, t 0 x= 2.98005491555891E-03, y= 1.49127039027490E-04 z= 2.97383128525201E-02, t= 1.97754215587208E-01 1 x= 8.21907637577074E-04, y= 3.03869902434743E-03 z= 3.13739200594087E-02, t= 2.08630697444505E-01 2 x=-2.77341241550057E-03, y= 1.81039448627537E-03 z= 3.30095272662973E-02, t= 2.19507179301801E-01 3 x=-2.70659836489072E-03, y=-2.18120461122312E-03 z= 3.46451344731860E-02, t= 2.30383661159098E-01 4 x= 1.29651333753636E-03, y=-3.40150378777555E-03 z= 3.62807416800746E-02, t= 2.41260143016394E-01 5 x= 2.53217828031854E-02, y= 1.47031248285148E-03 z= 2.84388448108359E-02, t= 2.52136624873691E-01 6 x= 6.70373299822560E-03, y= 2.55952481939030E-02 z= 2.96656185085583E-02, t= 2.63013106730987E-01 7 x=-2.31919397077539E-02, y= 1.48757285302579E-02 z= 3.08923922062806E-02, t= 2.73889588588284E-01 8 x=-2.21607807732554E-02, y=-1.81533495568647E-02 z= 3.21191659040029E-02, t= 2.84766070445580E-01 9 x= 1.08146675323841E-02, y=-2.77050996403719E-02 z= 3.33459396017253E-02, t= 2.95642552302877E-01 10 x= 4.51689701402223E-02, y= 2.98548821753235E-03 z= 9.84544026894265E-03, t= 3.06519034160173E-01 11 x= 1.15131427590713E-02, y= 4.54378709762631E-02 z= 1.01947946010664E-02, t= 3.17395516017470E-01 12 x=-4.10151626069116E-02, y= 2.58471037238780E-02 z= 1.05441489331902E-02, t= 3.28271997874766E-01 13 x=-3.84908387878038E-02, y=-3.20483324285474E-02 z= 1.08935032653140E-02, t= 3.39148479732063E-01 14 x= 1.91814869647966E-02, y=-4.80020309490219E-02 z= 1.12428575974377E-02, t= 3.50024961589359E-01 15 x= 4.97967614097874E-02, y= 3.69170139124708E-03 z=-2.19499888850803E-02, t= 3.60901443446656E-01 16 x= 1.22349667544167E-02, y= 4.99619867600284E-02 z=-2.26114953994252E-02, t= 3.71777925303952E-01 17 x=-4.50153669652142E-02, y= 2.78673428762470E-02 z=-2.32730019137701E-02, t= 3.82654407161249E-01 18 x=-4.15626478472590E-02, y=-3.51727877556355E-02 z=-2.39345084281149E-02, t= 3.93530889018545E-01 19 x= 2.11773052296168E-02, y=-5.17903096481542E-02 z=-2.45960149424598E-02, t= 4.04407370875841E-01 20 x= 3.15318926141121E-02, y= 2.59142604925538E-03 z=-5.42065117704958E-02, t= 4.15283852733138E-01 21 x= 7.46994634801066E-03, y= 3.15957958486789E-02 z=-5.56262061263897E-02, t= 4.26160334590434E-01 22 x=-2.84490553174695E-02, y= 1.72984842960736E-02 z=-5.70459004822837E-02, t= 4.37036816447731E-01 23 x=-2.58712991812066E-02, y=-2.22514561412578E-02 z=-5.84655948381776E-02, t= 4.47913298305027E-01 24 x= 1.34874672231278E-02, y=-3.22455886076567E-02 z=-5.98852891940715E-02, t= 4.58789780162324E-01 25 x= 4.79134563292288E-02, y= 4.32389188682661E-03 z= 4.39609976777021E-01, t= 1.73359952064038E-01 26 x= 1.09414527813035E-02, y= 4.79907718150468E-02 z= 4.49790418344489E-01, t= 1.77374603585521E-01 27 x=-4.32173624946507E-02, y= 2.58070870612648E-02 z= 4.59970859911957E-01, t= 1.81389255107004E-01 28 x=-3.87377063966394E-02, y=-3.38605260919957E-02 z= 4.70151301479424E-01, t= 1.85403906628487E-01 29 x= 2.06707732387232E-02, y=-4.83295518666820E-02 z= 4.80331743046892E-01, t= 1.89418558149970E-01 30 x= 3.30150959213619E-01, y= 3.24587719434010E-02 z= 3.65281214477906E-01, t= 1.93433209671453E-01 31 x= 7.26288730145978E-02, y= 3.30747529210091E-01 z= 3.72862522702919E-01, t= 1.97447861192936E-01 32 x=-2.98055549548036E-01, y= 1.74763373918250E-01 z= 3.80443830927932E-01, t= 2.01462512714419E-01 33 x=-2.63461638856582E-01, y=-2.34035448626728E-01 z= 3.88025139152945E-01, t= 2.05477164235902E-01 34 x= 1.43924958348298E-01, y=-3.29196482800326E-01 z= 3.95606447377959E-01, t= 2.09491815757385E-01 35 x= 5.30231782003501E-01, y= 5.64160241531277E-02 z= 1.10957099828014E-01, t= 2.13506467278868E-01 36 x= 1.12267862290997E-01, y= 5.31523906578760E-01 z= 1.13043472645293E-01, t= 2.17521118800351E-01 37 x=-4.79506165503554E-01, y= 2.76025104647665E-01 z= 1.15129845462572E-01, t= 2.21535770321834E-01 38 x=-4.18133534931173E-01, y=-3.77459550606994E-01 z= 1.17216218279850E-01, t= 2.25550421843317E-01 39 x= 2.33864843028100E-01, y=-5.23464424812200E-01 z= 1.19302591097129E-01, t= 2.29565073364800E-01 40 x= 5.39769742620943E-01, y= 6.18017079018271E-02 z=-2.44682827295638E-01, t= 2.33579724886283E-01 41 x= 1.09877719505141E-01, y= 5.41600781492923E-01 z=-2.48888313389782E-01, t= 2.37594376407766E-01 42 x=-4.89268712373359E-01, y= 2.76457471192930E-01 z=-2.53093799483925E-01, t= 2.41609027929249E-01 43 x=-4.21000158422726E-01, y=-3.86204547283094E-01 z=-2.57299285578069E-01, t= 2.45623679450732E-01 44 x= 2.41083216778351E-01, y=-5.28233698934358E-01 z=-2.61504771672213E-01, t= 2.49638330972215E-01 45 x= 3.18741081806704E-01, y= 3.90804960130383E-02 z=-5.61749056138165E-01, t= 2.53652982493698E-01 46 x= 6.22993201215676E-02, y= 3.20206371365050E-01 z=-5.70640048393589E-01, t= 2.57667634015181E-01 47 x=-2.89727723582807E-01, y= 1.60664247883917E-01 z=-5.79531040649014E-01, t= 2.61682285536664E-01 48 x=-2.46049348936833E-01, y=-2.29365138410983E-01 z=-5.88422032904438E-01, t= 2.65696937058147E-01 49 x= 1.44252981793420E-01, y=-3.09491334618491E-01 z=-5.97313025159862E-01, t= 2.69711588579630E-01 50 x= 6.84272708814056E-02, y= 8.94599787339966E-03 z= 5.82110058278290E-01, t=-4.67852068257711E-01 51 x= 1.28223832922628E-02, y= 6.88376930191534E-02 z= 5.90647672466372E-01, t=-4.74713898592158E-01 52 x=-6.23952536945842E-02, y= 3.39505799953346E-02 z= 5.99185286654453E-01, t=-4.81575728926604E-01 53 x=-5.23049214433332E-02, y=-4.95461571410902E-02 z= 6.07722900842535E-01, t=-4.88437559261050E-01 54 x= 3.13930016019437E-02, y=-6.59694757516567E-02 z= 6.16260515030617E-01, t=-4.95299389595497E-01 55 x= 4.23110170889851E-01, y= 5.87627036797166E-02 z= 4.61934884095456E-01, t=-5.02161219929943E-01 56 x= 7.58845345153743E-02, y= 4.26307146692049E-01 z= 4.68247037791170E-01, t=-5.09023050264390E-01 57 x=-3.87141988705425E-01, y= 2.06655375749213E-01 z= 4.74559191486885E-01, t=-5.15884880598836E-01 58 x=-3.20379859157494E-01, y=-3.08381859772812E-01 z= 4.80871345182599E-01, t=-5.22746710933283E-01 59 x= 1.96836046287033E-01, y=-4.05245113874359E-01 z= 4.87183498878313E-01, t=-5.29608541267729E-01 60 x= 6.52255893549098E-01, y= 9.59118162899533E-02 z= 1.31006538842357E-01, t=-5.36470371602175E-01 61 x= 1.11751931431735E-01, y= 6.58284176427705E-01 z= 1.32682203874061E-01, t=-5.43332201936622E-01 62 x=-5.99002765496187E-01, y= 3.13614694879439E-01 z= 1.34357868905766E-01, t=-5.50194032271068E-01 63 x=-4.89395724374783E-01, y=-4.78669477348401E-01 z= 1.36033533937471E-01, t=-5.57055862605515E-01 64 x= 3.07754935136838E-01, y=-6.20915382324727E-01 z= 1.37709198969175E-01, t=-5.63917692939961E-01 65 x= 6.41713836429240E-01, y= 9.96126498748738E-02 z=-2.99526973628109E-01, t=-5.70779523274408E-01 66 x= 1.04808243827600E-01, y= 6.48795197737191E-01 z=-3.03127844349321E-01, t=-5.77641353608854E-01 67 x=-5.91598498373138E-01, y= 3.03733007811092E-01 z=-3.06728715070533E-01, t=-5.84503183943300E-01 68 x=-4.77218592882907E-01, y=-4.74288290681378E-01 z=-3.10329585791745E-01, t=-5.91365014277747E-01 69 x=-4.14459244196422E-39, y= 2.83787229166536E-308 z= 2.83335668396214E-308, t=-2.90956675841230E-01 bound 1 U(i)= 4.26772677264757E+00 bound 2 U(i)= 4.27806563985581E+00 bound 3 U(i)= 4.28945348043329E+00 bound 4 U(i)= 4.30770066462903E+00 bound 5 U(i)= 4.34273541619534E+00 bound 6 U(i)= 4.36360566777129E+00 bound 7 U(i)= 4.33850595867056E+00 bound 8 U(i)= 4.32218494280270E+00 bound 9 U(i)= 4.36144633630214E+00 bound 10 U(i)= 4.41777596631118E+00 bound 11 U(i)= 4.44185832191988E+00 bound 12 U(i)= 4.38516515577116E+00 bound 13 U(i)= 4.34540510486151E+00 bound 14 U(i)= 4.40290775215832E+00 bound 15 U(i)= 4.46766720500928E+00 bound 16 U(i)= 4.49149435175063E+00 bound 17 U(i)= 4.42740903050181E+00 bound 18 U(i)= 4.38328253915401E+00 bound 19 U(i)= 4.44504932373942E+00 bound 20 U(i)= 4.49666620683926E+00 bound 21 U(i)= 4.51685708290466E+00 bound 22 U(i)= 4.48206456327607E+00 bound 23 U(i)= 4.46070846591161E+00 bound 24 U(i)= 4.50587822767189E+00 bound 25 U(i)= 4.79480897520822E+00 bound 26 U(i)= 4.82222434865451E+00 bound 27 U(i)= 4.76675575674362E+00 bound 28 U(i)= 4.73265014804763E+00 bound 29 U(i)= 4.79886293143741E+00 bound 30 U(i)= 5.07849692622230E+00 bound 31 U(i)= 5.13751395293509E+00 bound 32 U(i)= 4.61934859914690E+00 bound 33 U(i)= 4.26189863818886E+00 bound 34 U(i)= 4.59263502884508E+00 bound 35 U(i)= 5.11272237390147E+00 bound 36 U(i)= 5.18300809173569E+00 bound 37 U(i)= 4.30698099810325E+00 bound 38 U(i)= 3.72124997595192E+00 bound 39 U(i)= 4.24070168294358E+00 bound 40 U(i)= 4.82542953033676E+00 bound 41 U(i)= 4.88275931948235E+00 bound 42 U(i)= 3.98121676867779E+00 bound 43 U(i)= 3.38757667045048E+00 bound 44 U(i)= 3.91572578696441E+00 bound 45 U(i)= 4.27418472540301E+00 bound 46 U(i)= 4.30076513790711E+00 bound 47 U(i)= 3.78203268924369E+00 bound 48 U(i)= 3.43646420756666E+00 bound 49 U(i)= 3.74887139893409E+00 bound 50 U(i)= 4.49596610307026E+00 bound 51 U(i)= 4.51138745580790E+00 bound 52 U(i)= 4.41248897198384E+00 bound 53 U(i)= 4.35053032805132E+00 bound 54 U(i)= 4.42942842343222E+00 bound 55 U(i)= 4.77958921814836E+00 bound 56 U(i)= 4.80870356464421E+00 bound 57 U(i)= 4.11283108063101E+00 bound 58 U(i)= 3.67088097908022E+00 bound 59 U(i)= 4.10091844611317E+00 bound 60 U(i)= 4.74531296913746E+00 bound 61 U(i)= 4.77240754831806E+00 bound 62 U(i)= 3.63836155286290E+00 bound 63 U(i)= 2.95121751964930E+00 bound 64 U(i)= 3.61443719416432E+00 bound 65 U(i)= 4.31073033736813E+00 bound 66 U(i)= 4.32345722009482E+00 bound 67 U(i)= 3.20158206805730E+00 bound 68 U(i)= 2.52961109010519E+00 bound 69 U(i)= 3.74754806420756E+00 order= 4, m= 69, point= 0 point= 0, i= 1, ap*U=-4.57485707879493E+08 point= 0, i= 2, ap*U= 4.16854023394755E+08 point= 0, i= 3, ap*U=-2.51616576589800E+08 point= 0, i= 4, ap*U= 1.47341039938835E+08 point= 0, i= 5, ap*U=-1.08296521664239E+08 point= 0, i= 6, ap*U= 2.48271605158905E+08 point= 0, i= 7, ap*U=-2.48251485721417E+08 point= 0, i= 8, ap*U= 1.49273278048067E+08 point= 0, i= 9, ap*U=-1.01589573185812E+08 point= 0, i= 10, ap*U= 1.01478161648333E+08 point= 0, i= 11, ap*U=-2.36791011799739E+08 point= 0, i= 12, ap*U= 2.19439218687606E+08 point= 0, i= 13, ap*U=-1.34950466101034E+08 point= 0, i= 14, ap*U= 8.83552756952946E+07 point= 0, i= 15, ap*U=-8.32069642079968E+07 point= 0, i= 16, ap*U= 1.94571184440533E+08 point= 0, i= 17, ap*U=-1.72566605547283E+08 point= 0, i= 18, ap*U= 1.08700288965116E+08 point= 0, i= 19, ap*U=-6.36422098161725E+07 point= 0, i= 20, ap*U= 6.92402598127694E+07 point= 0, i= 21, ap*U=-1.52949806519043E+08 point= 0, i= 22, ap*U= 1.41073734125837E+08 point= 0, i= 23, ap*U=-9.28976541170008E+07 point= 0, i= 24, ap*U= 3.97393678194772E+07 point= 0, i= 25, ap*U=-5.82041449409658E+05 point= 0, i= 26, ap*U= 1.04062071348266E+06 point= 0, i= 27, ap*U=-1.18096951535336E+06 point= 0, i= 28, ap*U= 1.56289793236189E+06 point= 0, i= 29, ap*U=-8.31672523242414E+05 point= 0, i= 30, ap*U= 4.38694259358367E+04 point= 0, i= 31, ap*U= 7.65124537369394E+04 point= 0, i= 32, ap*U=-1.11728398217950E+05 point= 0, i= 33, ap*U=-2.44018488337840E+04 point= 0, i= 34, ap*U= 2.69115112365609E+04 point= 0, i= 35, ap*U=-8.94764562571179E+03 point= 0, i= 36, ap*U=-2.98172283748848E+04 point= 0, i= 37, ap*U= 3.68139543387297E+04 point= 0, i= 38, ap*U=-2.62750064565970E+02 point= 0, i= 39, ap*U=-1.01788015244549E+04 point= 0, i= 40, ap*U= 7.05231481673949E+03 point= 0, i= 41, ap*U=-5.75750613237775E+03 point= 0, i= 42, ap*U=-4.15356657665226E+03 point= 0, i= 43, ap*U= 1.53117872926510E+03 point= 0, i= 44, ap*U= 4.35663017445922E+03 point= 0, i= 45, ap*U=-9.11820094414762E+03 point= 0, i= 46, ap*U= 7.95502686615496E+03 point= 0, i= 47, ap*U= 1.82831501999239E+03 point= 0, i= 48, ap*U=-6.50631360746438E+03 point= 0, i= 49, ap*U= 4.62796090905335E+03 point= 0, i= 50, ap*U=-1.06186868613510E+05 point= 0, i= 51, ap*U= 3.43202791731947E+05 point= 0, i= 52, ap*U=-3.26343965229680E+05 point= 0, i= 53, ap*U= 3.10793682568213E+04 point= 0, i= 54, ap*U= 5.75253349703136E+04 point= 0, i= 55, ap*U=-3.51936194336308E+03 point= 0, i= 56, ap*U=-5.44445599581768E+04 point= 0, i= 57, ap*U= 5.17969247811983E+04 point= 0, i= 58, ap*U= 1.59496820145228E+02 point= 0, i= 59, ap*U=-1.98445362364775E+02 point= 0, i= 60, ap*U=-4.18100898154505E+03 point= 0, i= 61, ap*U= 2.60574667755360E+04 point= 0, i= 62, ap*U=-2.01156088789158E+04 point= 0, i= 63, ap*U= 2.52774373783340E+03 point= 0, i= 64, ap*U=-3.20149708778896E+02 point= 0, i= 65, ap*U= 3.03785598988541E+03 point= 0, i= 66, ap*U=-8.07633054319246E+03 point= 0, i= 67, ap*U= 5.96393917641485E+03 point= 0, i= 68, ap*U=-1.47273507185351E+03 point= 0, i= 69, ap*U=-1.15673827030102E+04 A matrix and rhs A( 0, 0)=-4.23126958034481E+07 rhs( 0)=-1.79912790734592E+08 vertex 0 x= 2.98005491555891E-03, y= 1.49127039027490E-04 z= 2.97383128525201E-02, t= 1.97754215587208E-01 computed= 4.25198128642824E+00, actual= 4.25198137905924E+00, error = 9.26310015358922E-08 maxerror = 9.26310015358922E-08 end pde44e_nuderiv4d.adb