diff --git a/moment_kinetics/test/kinetic_electron_tests.jl b/moment_kinetics/test/kinetic_electron_tests.jl index 6e58393e2..723da1be2 100644 --- a/moment_kinetics/test/kinetic_electron_tests.jl +++ b/moment_kinetics/test/kinetic_electron_tests.jl @@ -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, @@ -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 @@ -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.")