Skip to content
Snippets Groups Projects
Commit 273ca58f authored by Timothy A. Smith's avatar Timothy A. Smith
Browse files

more partial progress

parent 708ac8cc
No related branches found
No related tags found
2 merge requests!29Data fixtures,!28WENO flux bug
......@@ -40,7 +40,13 @@ class FluxDataSingle:
self.direction = self.dirs[direction]
# FIXME:
# - finish adding data from B (one case for all)
# - add data from A (one case for all)
# - rotate fluxes in input for all y,z
# - add rotation routines so y,z data is setup correctly
# (even though values are different, we are still working with x as the
# dominant frame, so need to swap at the end)
single_data = {}
single_data["Case a:x"] = FluxDataSingle(
states_str="2 4 4 4 20,1 1 1 1 5.5",
......@@ -197,17 +203,17 @@ single_data["Case b:z"] = FluxDataSingle(
char_fluxes_neg_str=
weno_weights_pos_str=
weno_weights_pos_str=0.9999965203327451,0.999999959029252,0.999999959029252,0.9999835300215946,0.9999999975371708},{3.2217041402391125*10^-6,3.793560951911702*10^-8,3.793560951911702*10^-8,0.000015247816238787948,2.2803933579018906*10^-9},{2.579631146567429*10^-7,3.0351383606886035*10^-9,3.0351383606886035*10^-9,1.2221621667041045*10^-6,1.824357365679662*10^-10
weno_weights_neg_str=
weno_weights_neg_str=0.999999997585736,0.9999999999716082,0.9999999999716082,0.9999909282503908,0.9999999999997268},{2.2354258683870077*10^-9,2.6288602366075745*10^-11,2.6288602366075745*10^-11,8.39888391917283*10^-6,2.5299566294114733*10^-13},{1.788382117901561*10^-10,2.1030934718853924*10^-12,2.1030934718853924*10^-12,6.728656898188601*10^-7,2.0239658022589593*10^-14
consistent_str=
consistent_str=-2.5,6.899999999999999,4.5,4.5,-26.75
dissipation_pos_str=
dissipation_pos_str=-0.17685508040141606,0.523974175351619,0.49761166908056337,0.49761166908056337,-2.495196349669178
dissipation_neg_str=
dissipation_neg_str=-1.6768551393982627,4.823974372333487,3.9976117689805504,3.9976117689805504,-22.145196359962693
weno_flux=
weno_flux=4.353710219799678,12.247948547685105,8.995223438061114,8.995223438061114,-51.390392709631875
direction="z")
single_data["Case c:x"] = FluxDataSingle(
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment