generate Q1 PoU test data with dg_p1

parent 0b6b1cd8
......@@ -24,7 +24,7 @@ set_log_levels({
make_problem = lambda : init_problem((1, 1))
base_config = {
'local_space_type': 'cg_p1',
'local_space_type': 'dg_p1',
'domain': domain,
'problem': make_problem()['id'],
'mu': 1,
......@@ -158,7 +158,8 @@ def expected_results_fixed_h_refined_H(global_grid_refines, num_H_refs_less, ini
'error': [0.18812850466547035, 0.10496923582575818, 0.04812058367396117, 1.1495276927521529e-15],
'estimate': [0.8043649753996472, 0.4426426990913186, 0.21654695989615283, 4.156337533325682e-15]}
else:
assert False, f'missing expected results for global_grid_refines={global_grid_refines}, num_H_refs_less={num_H_refs_less}, initial_macro_grid={initial_macro_grid}'
assert False, f"""missing expected results for global_grid_refines={global_grid_refines},
num_H_refs_less={num_H_refs_less}, initial_macro_grid={initial_macro_grid}"""
def test_fixed_H_refined_h_single_subdomain():
......@@ -177,9 +178,13 @@ def test_fixed_h_refined_H():
if __name__ == '__main__':
r = run_fixed_H_refined_h()
for kk, vv in r.items():
print(kk, vv)
run_fixed_H_refined_h(num_h_refs=8, num_macro_elements=[1, 1])
print('')
print('')
run_fixed_H_refined_h(num_h_refs=7, num_macro_elements=[2, 2])
print('')
print('')
run_fixed_H_refined_h(num_h_refs=6, num_macro_elements=[4, 4])
print('')
print('')
r = run_fixed_h_refined_H()
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment