Skip to content

Commit

Permalink
Updated expected output for fokker_planck_time_evolution_tests.jl
Browse files Browse the repository at this point in the history 
Previous output was affected by a default `vperp_bc = "periodic"` which
had unintended effects in `reconcile_element_boundaries_centered!()`,
which is used by `derivative!()`, etc.
  • Loading branch information
@johnomotani
johnomotani committed Nov 24, 2023
1 parent 1396824 commit 0905d65
Showing 1 changed file with 36 additions and 24 deletions.
60 changes: 36 additions & 24 deletions test/fokker_planck_time_evolution_tests.jl
View file
Original file line number Diff line number Diff line change
Expand Up @@ -36,39 +36,51 @@ end

const expected =
expected_data(
vec([-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 ]),
vec([0.155051025721682 0.644948974278318 1.000000000000000 1.500000000000000 2.000000000000000 2.500000000000000 3.000000000000000 ]),
[-3.0, -2.5, -2.0, -1.5, -1.0, -0.5, 0.0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0],
[0.155051025721682, 0.644948974278318, 1.000000000000000, 1.500000000000000, 2.000000000000000, 2.500000000000000, 3.000000000000000],
# Expected phi:
vec([-1.267505494648937 -1.275240686416660 ]),
[-1.267505494648937, -1.275683298550937],
# Expected n_charged:
vec([0.281533032234007 0.279363721076024 ]),
[0.2815330322340072, 0.2792400986636072],
# Expected upar_charged:
vec([0.0 0.0 ]),
[0.0, 0.0],
# Expected ppar_charged:
vec([0.179822802480489 0.147624463307870 ]),
[0.17982280248048935, 0.14891126175332367],
# Expected pperp_charged
vec([0.143401466675068 0.158821758223731 ]),
[0.14340146667506784, 0.1581377822859991],
# Expected qpar_charged
vec([0.0 0.0 ]),
[0.0, 0.0],
# Expected v_t_charged
vec([1.051172608301042 1.053709630977256 ]),
[1.0511726083010418, 1.0538509291794658],
# Expected dSdt
vec([0.0 0.000008786074305 ]),
[0.0, 1.1853081348031516e-5],
# Expected f_charged:
[ 0.000000000000000 ; 0.000619960016181 ; 0.005882016862627 ; 0.033848669972256 ; 0.118138103172003 ; 0.229579465209362 ; 0.000000000000000 ; 0.229579465209362 ; 0.118138103172003 ; 0.033848669972256 ; 0.005882016862627 ; 0.000619960016181 ; 0.000000000000000 ;;
0.000000000000000 ; 0.000478053009971 ; 0.004535640674379 ; 0.026100809957764 ; 0.091096642266611 ; 0.177029407552679 ; 0.000000000000000 ; 0.177029407552679 ; 0.091096642266611 ; 0.026100809957764 ; 0.004535640674379 ; 0.000478053009971 ; 0.000000000000000 ;;
0.000000000000000 ; 0.000266581711212 ; 0.002529256854782 ; 0.014554868262370 ; 0.050799175561231 ; 0.098718764270694 ; 0.000000000000000 ; 0.098718764270694 ; 0.050799175561231 ; 0.014554868262370 ; 0.002529256854782 ; 0.000266581711212 ; 0.000000000000000 ;;
0.000000000000000 ; 0.000076376939017 ; 0.000724644221386 ; 0.004170039574837 ; 0.014554207474836 ; 0.028283399503664 ; 0.000000000000000 ; 0.028283399503664 ; 0.014554207474836 ; 0.004170039574837 ; 0.000724644221386 ; 0.000076376939017 ; 0.000000000000000 ;;
0.000000000000000 ; 0.000013272321882 ; 0.000125924283949 ; 0.000724644221264 ; 0.002529142026698 ; 0.004914917865936 ; 0.000000000000000 ; 0.004914917865936 ; 0.002529142026698 ; 0.000724644221264 ; 0.000125924283949 ; 0.000013272321882 ; 0.000000000000000 ;;
0.000000000000000 ; 0.000001398892434 ; 0.000013272321882 ; 0.000076376939004 ; 0.000266569608421 ; 0.000518028531855 ; 0.000000000000000 ; 0.000518028531855 ; 0.000266569608421 ; 0.000076376939004 ; 0.000013272321882 ; 0.000001398892434 ; 0.000000000000000 ;;
0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ;;;
0.000000000000000 ; 0.000138780288865 ; 0.005791061761509 ; 0.027815242205125 ; 0.081445074775077 ; 0.152354938126957 ; 0.186355844666422 ; 0.152354938126957 ; 0.081445074775077 ; 0.027815242205125 ; 0.005791061761509 ; 0.000138780288865 ; 0.000000000000000 ;;
0.000000000000000 ; 0.000041802336103 ; 0.004577414398302 ; 0.022242056573360 ; 0.065733717950559 ; 0.123790379285355 ; 0.151688746487743 ; 0.123790379285355 ; 0.065733717950559 ; 0.022242056573360 ; 0.004577414398302 ; 0.000041802336103 ; 0.000000000000000 ;;
0.000000000000000 ; -0.000102715219225 ; 0.002768824456286 ; 0.013936837294540 ; 0.042320490360991 ; 0.081223176170832 ; 0.100027475688989 ; 0.081223176170832 ; 0.042320490360991 ; 0.013936837294540 ; 0.002768824456286 ; -0.000102715219225 ; 0.000000000000000 ;;
0.000000000000000 ; -0.000164682482058 ; 0.000767345547797 ; 0.004535511738255 ; 0.014877992136829 ; 0.029791500386470 ; 0.037097641503996 ; 0.029791500386470 ; 0.014877992136829 ; 0.004535511738255 ; 0.000767345547797 ; -0.000164682482058 ; 0.000000000000000 ;;
0.000000000000000 ; -0.000091767217558 ; 0.000022353567868 ; 0.000693446611185 ; 0.002889786006700 ; 0.006284514040367 ; 0.007979739613849 ; 0.006284514040367 ; 0.002889786006700 ; 0.000693446611185 ; 0.000022353567868 ; -0.000091767217558 ; 0.000000000000000 ;;
0.000000000000000 ; -0.000032024960906 ; -0.000089243114443 ; -0.000201682249005 ; -0.000175449913808 ; 0.000111463339173 ; 0.000290773816557 ; 0.000111463339173 ; -0.000175449913808 ; -0.000201682249005 ; -0.000089243114443 ; -0.000032024960906 ; 0.000000000000000 ;;
0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ; 0.000000000000000 ])
[0.0 0.0 0.0 0.0 0.0 0.0 0.0;
0.0006199600161806666 0.00047805300997075977 0.0002665817112117718 7.637693901737056e-5 1.3272321881722645e-5 1.3988924344690309e-6 0.0;
0.005882016862626724 0.0045356406743786385 0.002529256854781707 0.0007246442213864763 0.00012592428394890537 1.3272321881722645e-5 0.0;
0.03384866997225574 0.026100809957763767 0.01455486826237011 0.004170039574837177 0.000724644221263874 7.637693900444835e-5 0.0;
0.11813810317200342 0.09109664226661075 0.05079917556123135 0.01455420747483572 0.0025291420266981487 0.0002665696084208068 0.0;
0.22957946520936198 0.17702940755267918 0.0987187642706944 0.0282833995036643 0.004914917865936108 0.0005180285318549189 0.0;
0.0 0.0 0.0 0.0 0.0 0.0 0.0;
0.22957946520936204 0.1770294075526792 0.0987187642706944 0.0282833995036643 0.004914917865936108 0.0005180285318549189 0.0;
0.11813810317200349 0.0910966422666108 0.050799175561231376 0.01455420747483573 0.0025291420266981487 0.0002665696084208068 0.0;
0.03384866997225574 0.026100809957763767 0.01455486826237011 0.004170039574837177 0.000724644221263874 7.637693900444835e-5 0.0;
0.005882016862626724 0.0045356406743786385 0.002529256854781707 0.0007246442213864763 0.00012592428394890537 1.3272321881722645e-5 0.0;
0.0006199600161806666 0.00047805300997075977 0.0002665817112117718 7.637693901737056e-5 1.3272321881722645e-5 1.3988924344690309e-6 0.0;
0.0 0.0 0.0 0.0 0.0 0.0 0.0;;;
0.0 0.0 0.0 0.0 0.0 0.0 0.0;
0.0001712743622973216 7.105465094508053e-5 -7.829380680167827e-5 -0.00015364081956318698 -9.097098213067502e-5 -3.311284120491419e-5 0.0;
0.005883280697248667 0.004667594200766182 0.002855965521103658 0.0008138347136178689 2.44260649525292e-5 -9.753249634264602e-5 0.0;
0.02792209301450194 0.022385716644538384 0.01413535091105969 0.004677801530322722 0.0007105315221401102 -0.00022400635166536323 0.0;
0.08117458037332098 0.06563459159004267 0.04247673844050208 0.015087784332275832 0.0029056314178876035 -0.00023019804543218203 0.0;
0.15133793170654106 0.12313903060106579 0.08111673445361306 0.029975277983613262 0.00626735398468981 7.553501812465833e-6 0.0;
0.18493902160817713 0.15073513412904313 0.09976414473955808 0.037251581926306565 0.007941836186495122 0.00016196175024033304 0.0;
0.15133793170654092 0.12313903060106571 0.08111673445361306 0.02997527798361324 0.006267353984689816 7.553501812469816e-6 0.0;
0.081174580373321 0.06563459159004267 0.042476738440502065 0.015087784332275821 0.002905631417887614 -0.0002301980454321778 0.0;
0.027922093014501933 0.022385716644538384 0.014135350911059698 0.004677801530322729 0.0007105315221401184 -0.00022400635166536134 0.0;
0.005883280697248667 0.004667594200766184 0.002855965521103663 0.0008138347136178759 2.4426064952530956e-5 -9.753249634264635e-5 0.0;
0.0001712743622973275 7.105465094508572e-5 -7.829380680167411e-5 -0.00015364081956318568 -9.097098213067551e-5 -3.311284120491447e-5 0.0;
0.0 0.0 0.0 0.0 0.0 0.0 0.0])
###########################################################################################
# to modify the test, with a new expected f, print the new f using the following commands
# in an interative Julia REPL. The path is the path to the .dfns file.
Expand Down

0 comments on commit 0905d65

Please sign in to comment.