[scripts] minor

parent ec29b16e
......@@ -82,12 +82,12 @@ def compute_on_single_refinemet(cfg, make_pou, pou_id):
def run_fixed_H_refined_h(
num_h_refs=3,
macro_grid=[3, 3],
num_macro_elements=[3, 3],
make_pou=lambda dd: ContinuousLagrangePartitionOfUnity(dd),
pou_id='Q1PoU',
DATASET_ID=FILENAME + '.fixed_H_refined_h'):
config = dict(base_config)
config['num_macro_elements'] = macro_grid
config['num_macro_elements'] = num_macro_elements
config['num_h_refs'] = num_h_refs
new_dataset(DATASET_ID, **config)
results = {'num_subdomains': [], 'num_DoFs': [], 'error': [], 'estimate': []}
......@@ -136,19 +136,19 @@ def run_fixed_h_refined_H(
return results
def expected_results_fixed_H_refined_h(num_h_refs, macro_grid):
if num_h_refs == 3 and macro_grid == [1, 1]:
def expected_results_fixed_H_refined_h(num_h_refs, num_macro_elements):
if num_h_refs == 3 and num_macro_elements == [1, 1]:
return {'num_subdomains': [1, 1, 1, 1],
'num_DoFs': [4, 9, 25, 81],
'error': [2.2827713033690556e-17, 0.17829090592609537, 0.18394776571648075, 0.18812850466547035],
'estimate': [1.379885348450665e-16, 0.7255116728902438, 0.7943998932316314, 0.8043649753996472]}
elif num_h_refs == 3 and macro_grid == [3, 3]:
elif num_h_refs == 3 and num_macro_elements == [3, 3]:
return {'num_subdomains': [9, 9, 9, 9],
'num_DoFs': [36, 81, 225, 729],
'error': [2.437083428529774e-16, 0.06342823397024529, 0.07209377002856301, 0.07411333407007065],
'estimate': [1.2525548119768789e-15, 0.28110033383280547, 0.30527579252513726, 0.3023629356335546]}
else:
assert False, f'missing expected results for num_h_refs={num_h_refs}, macro_grid={macro_grid}'
assert False, f'missing expected results for num_h_refs={num_h_refs}, num_macro_elements={num_macro_elements}'
def expected_results_fixed_h_refined_H(global_grid_refines, num_H_refs_less, initial_macro_grid):
......
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