Skip to content

Commit

Permalink
Tweak parameters of kinetic electron test and update expected results
Browse files Browse the repository at this point in the history 
Better electron timestep size for kinetic electron test
  • Loading branch information
@johnomotani
johnomotani committed Sep 30, 2024
1 parent a169682 commit e75d6ac
Showing 1 changed file with 69 additions and 69 deletions.
138 changes: 69 additions & 69 deletions moment_kinetics/test/kinetic_electron_tests.jl
View file
Original file line number Diff line number Diff line change
Expand Up @@ -108,10 +108,10 @@ kinetic_input["timestepping"] = OptionsDict("type" => "PareschiRusso2(2,2,2)",
)

kinetic_input["electron_timestepping"] = OptionsDict("nstep" => 5000000,
"dt" => 1.0e-5,
"dt" => 5.0e-6,
"nwrite" => 10000,
"nwrite_dfns" => 100000,
"decrease_dt_iteration_threshold" => 1000,
"decrease_dt_iteration_threshold" => 5000,
"increase_dt_iteration_threshold" => 0,
"cap_factor_ion_dt" => 10.0,
"initialization_residual_value" => 1.0e10,
Expand Down Expand Up @@ -172,72 +172,72 @@ function run_test()

# Regression test
# Benchmark data generated in serial on Linux
expected_Ez = [-0.5990683230706185 -0.6042082363495851;
-0.4944296396481284 -0.49692371894536586;
-0.30889032954504736 -0.3090990586904173;
-0.2064830747303776 -0.20700297720010077;
-0.21232457328748663 -0.2132748045598696;
-0.18233875912042674 -0.18276920923500758;
-0.16711429522309232 -0.1674324272230308;
-0.16920776495088916 -0.16937992443371716;
-0.1629417555658927 -0.16309341722744303;
-0.16619150334079993 -0.16633546753735795;
-0.15918194883360942 -0.15931554370144113;
-0.14034706409006803 -0.140469880250037;
-0.12602184032280567 -0.12613381924054493;
-0.10928716440800472 -0.10938345602505639;
-0.07053969674257217 -0.0706024520856333;
-0.0249577746169536 -0.024980098134854842;
-2.8327303308330514e-15 -1.599033453711614e-10;
0.024957774616960776 0.02498009782733815;
0.07053969674257636 0.07060245115760132;
0.10928716440799909 0.10938345732933795;
0.1260218403227975 0.1261338225947928;
0.1403470640900294 0.14046988178255268;
0.1591819488336015 0.15931556545456152;
0.16619150334082114 0.1663353993955267;
0.16294175556587748 0.16309307445724816;
0.16920776495090983 0.1693805039915145;
0.1671142952230893 0.1674318780154963;
0.1823387591204167 0.18277420263305205;
0.21232457328753865 0.21326329266495697;
0.20648307473037922 0.20700517064938181;
0.3088903295450278 0.3091144991453789;
0.4944296396481271 0.49684270193048663;
0.5990683230705801 0.6040141042995336]
expected_vthe = [27.08122333732766 27.083668406411196;
27.087128061238488 27.08840157326006;
27.090525010446868 27.090443986816897;
27.091202856161452 27.0914901864659;
27.09265674296987 27.093297466503625;
27.093298138334738 27.09337068853881;
27.094377689895747 27.094548022524926;
27.09501542767647 27.095170446421935;
27.095227831625575 27.095304545176944;
27.095420218946682 27.09555512096241;
27.095754478126825 27.095876494374046;
27.096054218271775 27.096188914603825;
27.096199500698383 27.096294431476554;
27.09632238748948 27.096423453543142;
27.096502792691805 27.096594041947167;
27.096597492028636 27.096694147970585;
27.096610989303674 27.096702959927107;
27.096597492397745 27.096694148339555;
27.096502794930903 27.096594044186332;
27.096322390449956 27.09642345650393;
27.096199499205674 27.096294429984052;
27.09605421760595 27.096188913937898;
27.095754438597055 27.095876454845936;
27.09542019655419 27.095555098545283;
27.095228009815475 27.095304723869976;
27.095015217848847 27.09517023619458;
27.094377437638478 27.09454777080713;
27.093294828184774 27.093367377705533;
27.092639150183448 27.09327987116632;
27.0912092735745 27.091496606764487;
27.09048496370012 27.090403937882265;
27.08714601914595 27.08841951855733;
27.08144246136634 27.08388753119234]
expected_Ez = [-0.5990683230706185 -0.604849806235434;
-0.4944296396481284 -0.49739671491727844;
-0.30889032954504736 -0.30924318765687464;
-0.2064830747303776 -0.20682475071884582;
-0.21232457328748663 -0.21299072376949116;
-0.18233875912042674 -0.18256905463006085;
-0.16711429522309232 -0.1673112962636778;
-0.16920776495088916 -0.1693227707158167;
-0.1629417555658927 -0.16304933113558318;
-0.16619150334079993 -0.16629539618289285;
-0.15918194883360942 -0.1592799009526323;
-0.14034706409006803 -0.140437217833422;
-0.12602184032280567 -0.12610387949683538;
-0.10928716440800472 -0.10935785133612701;
-0.07053969674257217 -0.07058573063123225;
-0.0249577746169536 -0.024974174596810936;
-2.8327303308330514e-15 -1.441401377024236e-10;
0.024957774616960776 0.02497417427570905;
0.07053969674257636 0.07058572965952663;
0.10928716440799909 0.10935785264749627;
0.1260218403227975 0.12610388283669527;
0.1403470640900294 0.1404372197714126;
0.1591819488336015 0.15927992284761766;
0.16619150334082114 0.1662953275454769;
0.16294175556587748 0.1630489871826757;
0.16920776495090983 0.1693233489685909;
0.1671142952230893 0.16731075590341918;
0.1823387591204167 0.1825740389953209;
0.21232457328753865 0.21297925141919793;
0.20648307473037922 0.20682690396901446;
0.3088903295450278 0.30925854110074175;
0.4944296396481271 0.49731601862961966;
0.5990683230705801 0.6046564647413697]
expected_vthe = [27.08102229345079 27.08346736523219;
27.087730258479823 27.089003820908527;
27.091898844901323 27.09181784480061;
27.092455021687254 27.092742387764524;
27.09350739287911 27.094148133125078;
27.093817059011126 27.093889601910092;
27.09443981315218 27.094610141036807;
27.09484177005478 27.094996783801374;
27.094985914811055 27.0950626278904;
27.095122128675094 27.09525702879687;
27.09536357532887 27.09548558966323;
27.095582117080163 27.095716810823177;
27.09568783962135 27.09578276803757;
27.0957775472326 27.095878610625554;
27.095909169276535 27.09600041573683;
27.095978269355648 27.096074922150624;
27.095988166679223 27.096080134292468;
27.095978269713978 27.096074922508883;
27.095909171602027 27.096000418062378;
27.09577755035281 27.095878613746088;
27.095687838236376 27.095782766652857;
27.09558211622511 27.095716809968053;
27.09536353456768 27.09548554890375;
27.095122105596843 27.095257005693973;
27.094986093051983 27.09506280663278;
27.094841563692096 27.094996577040796;
27.094439553087433 27.094609881510113;
27.093813728418613 27.09388627063591;
27.093489818175936 27.094130555874184;
27.09246140309467 27.092748772044477;
27.09185903467811 27.09177803239964;
27.08774827015981 27.089021820036553;
27.081240668889404 27.0836857414255]

if expected_Ez == nothing
# Error: no expected input provided
Expand All @@ -257,7 +257,7 @@ function run_test()
# Iteration counts are fairly inconsistent, but it's good to check that they at
# least don't unexpectedly increase by an order of magnitude.
# Expected iteration count is from a serial run on Linux.
expected_electron_advance_linear_iterations = 10695
expected_electron_advance_linear_iterations = 11394
@test electron_advance_linear_iterations < 2.0 * expected_electron_advance_linear_iterations
if !(electron_advance_linear_iterations < 2.0 * expected_electron_advance_linear_iterations)
println("electron_advance_linear_iterations=$electron_advance_linear_iterations was greater than twice the expected $expected_electron_advance_linear_iterations.")
Expand Down

0 comments on commit e75d6ac

Please sign in to comment.