Commit 10aefd34 authored by Aytekin Gel's avatar Aytekin Gel
Browse files

PSUADE based deterministic calibration files 

Exp_n21.dat is the experimental dataset (in this case the 21 sample based analytical solution) used to guide the calibration

evalRSM_QUAD.c is the C based surrogate model exported by PSUADE initially and then revised to perform residual calculation

psuade.in is the script to run the optimization to find the set of five parameters that minimize the residuals.
parent 4594dba1

Case1_ParticleSettling/DeterministicCalibration/PSUADE/3_Calibration/Deterministic/Exp_n21.dat

0 → 100644
+22 −0
Original line number Diff line number Diff line 
%eval_id  x1          y2         y2 variances

1	0.25	0.088980934	0.001

2	0.24	0.088593673	0.001

3	0.23	0.088014071	0.001

4	0.22	0.087231143	0.001

5	0.21	0.0862336	0.001

6	0.2	0.085009842	0.001

7	0.19	0.08354796	0.001

8	0.18	0.081835729	0.001

9	0.17	0.079860606	0.001

10	0.16	0.077609728	0.001

11	0.15	0.075069909	0.001

12	0.14	0.072227635	0.001

13	0.13	0.069069065	0.001

14	0.12	0.065580023	0.001

15	0.11	0.061746	0.001

16	0.1	0.057552148	0.001

17	0.09	0.052983279	0.001

18	0.08	0.048023859	0.001

19	0.07	0.042658011	0.001

20	0.06	0.036869505	0.001

21	0.05	0.03064176	0.001

Case1_ParticleSettling/DeterministicCalibration/PSUADE/3_Calibration/Deterministic/evalRSM_QUAD.c

0 → 100644
+252 −0

File added.

Preview size limit exceeded, changes collapsed.

Case1_ParticleSettling/DeterministicCalibration/PSUADE/3_Calibration/Deterministic/psuade.in

0 → 100644
+34 −0
Original line number Diff line number Diff line 
PSUADE

INPUT

   dimension = 5

   variable 1 t1:P_0  =   4.8779439400000002e-01   2.0000224679999999e+01

   variable 2 t2:beta  =   2.0000728950000002e+00   5.0000000000000000e+00

   variable 3 t3:StatWeight  =   2.9636698739999998e+00   2.0007980220000000e+01

   variable 4 t4:ep_gstar  =   3.5015843499999999e-01   4.9974982499999998e-01

   variable 5 t5:VelfacCoeff  =   5.0113757199999998e-01   9.9986000399999997e-01

END

OUTPUT

   dimension = 1

   variable 1 Y1

END

METHOD

   sampling = LPTAU

   num_samples = 50

END

APPLICATION

   driver = ./evalRSM_QUAD

   opt_driver = ./evalRSM_QUAD

END

ANALYSIS

   analyzer output_id  = 1

#  analyzer rstype = MARS

## bobyqa: bound-constrained optimization

   optimization method = bobyqa

   optimization num_local_minima = 10

#  optimization print_level = 0

#  optimization num_fmin = 0

#  optimization output_id = 0

   optimization max_feval = 10000

   optimization deltax = 1.0e-6

END

END