diff --git a/corrector/data/checkpoint/0/u.h5 b/corrector/data/checkpoint/0/u.h5 new file mode 100644 index 0000000..ce66d52 Binary files /dev/null and b/corrector/data/checkpoint/0/u.h5 differ diff --git a/corrector/data/checkpoint/1/u.h5 b/corrector/data/checkpoint/1/u.h5 new file mode 100644 index 0000000..7aa2643 Binary files /dev/null and b/corrector/data/checkpoint/1/u.h5 differ diff --git a/corrector/data/checkpoint/10/u.h5 b/corrector/data/checkpoint/10/u.h5 new file mode 100644 index 0000000..bb04d8c Binary files /dev/null and b/corrector/data/checkpoint/10/u.h5 differ diff --git a/corrector/data/checkpoint/11/u.h5 b/corrector/data/checkpoint/11/u.h5 new file mode 100644 index 0000000..040fb55 Binary files /dev/null and b/corrector/data/checkpoint/11/u.h5 differ diff --git a/corrector/data/checkpoint/12/u.h5 b/corrector/data/checkpoint/12/u.h5 new file mode 100644 index 0000000..4fdc91d Binary files /dev/null and b/corrector/data/checkpoint/12/u.h5 differ diff --git a/corrector/data/checkpoint/13/u.h5 b/corrector/data/checkpoint/13/u.h5 new file mode 100644 index 0000000..64f4044 Binary files /dev/null and b/corrector/data/checkpoint/13/u.h5 differ diff --git a/corrector/data/checkpoint/14/u.h5 b/corrector/data/checkpoint/14/u.h5 new file mode 100644 index 0000000..51d5f1d Binary files /dev/null and b/corrector/data/checkpoint/14/u.h5 differ diff --git a/corrector/data/checkpoint/15/u.h5 b/corrector/data/checkpoint/15/u.h5 new file mode 100644 index 0000000..5706cde Binary files /dev/null and b/corrector/data/checkpoint/15/u.h5 differ diff --git a/corrector/data/checkpoint/16/u.h5 b/corrector/data/checkpoint/16/u.h5 new file mode 100644 index 0000000..e4752d1 Binary files /dev/null and b/corrector/data/checkpoint/16/u.h5 differ diff --git a/corrector/data/checkpoint/17/u.h5 b/corrector/data/checkpoint/17/u.h5 new file mode 100644 index 0000000..4711665 Binary files /dev/null and b/corrector/data/checkpoint/17/u.h5 differ diff --git a/corrector/data/checkpoint/18/u.h5 b/corrector/data/checkpoint/18/u.h5 new file mode 100644 index 0000000..a2c2e35 Binary files /dev/null and b/corrector/data/checkpoint/18/u.h5 differ diff --git a/corrector/data/checkpoint/19/u.h5 b/corrector/data/checkpoint/19/u.h5 new file mode 100644 index 0000000..6ab74cb Binary files /dev/null and b/corrector/data/checkpoint/19/u.h5 differ diff --git a/corrector/data/checkpoint/2/u.h5 b/corrector/data/checkpoint/2/u.h5 new file mode 100644 index 0000000..d0a1e4d Binary files /dev/null and b/corrector/data/checkpoint/2/u.h5 differ diff --git a/corrector/data/checkpoint/20/u.h5 b/corrector/data/checkpoint/20/u.h5 new file mode 100644 index 0000000..4c0c7ba Binary files /dev/null and b/corrector/data/checkpoint/20/u.h5 differ diff --git a/corrector/data/checkpoint/21/u.h5 b/corrector/data/checkpoint/21/u.h5 new file mode 100644 index 0000000..57fdc84 Binary files /dev/null and b/corrector/data/checkpoint/21/u.h5 differ diff --git a/corrector/data/checkpoint/22/u.h5 b/corrector/data/checkpoint/22/u.h5 new file mode 100644 index 0000000..0e4c805 Binary files /dev/null and b/corrector/data/checkpoint/22/u.h5 differ diff --git a/corrector/data/checkpoint/23/u.h5 b/corrector/data/checkpoint/23/u.h5 new file mode 100644 index 0000000..f3a3d85 Binary files /dev/null and b/corrector/data/checkpoint/23/u.h5 differ diff --git a/corrector/data/checkpoint/24/u.h5 b/corrector/data/checkpoint/24/u.h5 new file mode 100644 index 0000000..a447287 Binary files /dev/null and b/corrector/data/checkpoint/24/u.h5 differ diff --git a/corrector/data/checkpoint/3/u.h5 b/corrector/data/checkpoint/3/u.h5 new file mode 100644 index 0000000..6644e2a Binary files /dev/null and b/corrector/data/checkpoint/3/u.h5 differ diff --git a/corrector/data/checkpoint/4/u.h5 b/corrector/data/checkpoint/4/u.h5 new file mode 100644 index 0000000..d531212 Binary files /dev/null and b/corrector/data/checkpoint/4/u.h5 differ diff --git a/corrector/data/checkpoint/5/u.h5 b/corrector/data/checkpoint/5/u.h5 new file mode 100644 index 0000000..7cd846b Binary files /dev/null and b/corrector/data/checkpoint/5/u.h5 differ diff --git a/corrector/data/checkpoint/6/u.h5 b/corrector/data/checkpoint/6/u.h5 new file mode 100644 index 0000000..d5b83bb Binary files /dev/null and b/corrector/data/checkpoint/6/u.h5 differ diff --git a/corrector/data/checkpoint/7/u.h5 b/corrector/data/checkpoint/7/u.h5 new file mode 100644 index 0000000..b015247 Binary files /dev/null and b/corrector/data/checkpoint/7/u.h5 differ diff --git a/corrector/data/checkpoint/8/u.h5 b/corrector/data/checkpoint/8/u.h5 new file mode 100644 index 0000000..fa30796 Binary files /dev/null and b/corrector/data/checkpoint/8/u.h5 differ diff --git a/corrector/data/checkpoint/9/u.h5 b/corrector/data/checkpoint/9/u.h5 new file mode 100644 index 0000000..8a85132 Binary files /dev/null and b/corrector/data/checkpoint/9/u.h5 differ diff --git a/corrector/input/corrector.yaml b/corrector/input/corrector.yaml new file mode 100644 index 0000000..870ed40 --- /dev/null +++ b/corrector/input/corrector.yaml @@ -0,0 +1,57 @@ +mesh: './mesh/volunteer.h5' +density: 1.06 +dynamic_viscosity: 0.038 + +io: + write_path: './results' + write_xdmf: True + write_timeseries: false + write_checkpoint: True + log: false + +postprocessing: + enabled: false + debug_xdmf: false + pressure_drop: + method: boundaries # slices, boundaries, spheres + boundaries: [2,3] + slices: + - + spheres: + - + timeseries: + # FOR ALL THAT IS LOVED & HOLY, do not put the .h5 file extension here! + +measurements: +- 'all' +- './data/' + +dt: 0.0343 +pressure_checkpoints: + + +methods: +- COR + +COR: + time_schemes: + - implicit + fix_const: 2 # 0: false, 1: substract zero_point, 2: nullspace, 3: both + convection_schemes: + - standard + zero_point: + boundary_type: 'dirichlet' + theta: 1.0 + solver: + type: 'newton' + max_iteration: 30 + tol: 1e-7 + mode: 'auto' + reynolds_continuation: + enabled: true + Nfrac: {20: 'normal' , 50: [1,2,3,4,5,6] } + fe_space: P2P1 + stabilization_parameter: 0.01 + compute_norms: false + linear_solver: + method: mumps diff --git a/corrector/mesh/volunteer.h5 b/corrector/mesh/volunteer.h5 new file mode 100644 index 0000000..9c2048d Binary files /dev/null and b/corrector/mesh/volunteer.h5 differ diff --git a/corrector/run_estimate.py b/corrector/run_estimate.py new file mode 100644 index 0000000..48aea11 --- /dev/null +++ b/corrector/run_estimate.py @@ -0,0 +1,30 @@ +import sys +import logging +import os.path +from dolfin import * +from dPdirectestim.dPdirectestim import * + + +logging.getLogger().setLevel(logging.INFO) +#set_log_level(PROGRESS) + +parameters['form_compiler']['optimize'] = True +parameters['form_compiler']['cpp_optimize'] = True +parameters['form_compiler']['cpp_optimize_flags'] = \ + '-O3 -ffast-math -march=native' + +if len(sys.argv) > 1 and os.path.exists(sys.argv[1]): + inpfile = sys.argv[1] + print('Found input file ' + inpfile) +else: + raise Exception('No input file given') + +estimator = DPDirectEstim(inpfile) + +estimator.estimate() + + +if '2019' in dolfin.__version__: + list_timings(TimingClear.clear, [TimingType.wall]) +else: + list_timings(TimingClear_clear, [TimingType_wall])