StochasticNumerics / mimclib

A software library for UQ methods
GNU General Public License v2.0
6 stars 8 forks source link

Why TOL in tbl_runs-tbl_iters is different from TOLs defined in pythyn #106

Open litvinen opened 7 years ago

litvinen commented 7 years ago

Hello,

In python script I define TOLs

TOLs TOLs=0.005*np.sqrt(2)**-np.arange(0.,8.) array([ 0.005 , 0.00353553, 0.0025 , 0.00176777, 0.00125 , 0.00088388, 0.000625 , 0.00044194])

But, TOLs in DB are different See please below TOL={0.32, 0.08 , 0.04}. Why? q2

haji-ali commented 7 years ago

How do you set the TOLs that you define?

litvinen commented 7 years ago

my echo.py script

!/usr/bin/python

import numpy as np import argparse

if name == "main": parser = argparse.ArgumentParser(add_help=True) parser.register('type', 'bool', lambda v: v.lower() in ("yes", "true", "t", "1")) parser.add_argument("-tries", type=int, action="store", default=0, help="Number of realizations") parser.add_argument("-mimc_min_dim", type=int, action="store", default=1, help="MIMC dim")

args, unknowns = parser.parse_known_args()

dim = args.mimc_min_dim
base = "run.py -mimc_TOL {TOL} -mimc_max_TOL 0.5  -mimc_min_dim 1 -qoi_seed {seed} \
-mimc_theta 0.5 -mimc_M0 3 -mimc_h0inv 3 \
-mimc_w 2 -mimc_s 4 -mimc_gamma 2 -mimc_beta 2 \
-mimc_bayes_fit_lvls 3 -mimc_moments 4 \
        -mimc_bayesian True -mimc_reuse_samples True ".format(TOL="{TOL}", seed="{seed}")
base += " ".join(unknowns)

if args.tries == 0:
    cmd_single = "python " + base + " -mimc_verbose 10 -db False "
    print(cmd_single.format(seed=0, bayesian=False, TOL=0.001))
else:
    cmd_multi = " python " + base + " -mimc_verbose 10 -db True -db_tag {tag} > a2.txt"
    TOLs = 0.005*np.sqrt(2)**-np.arange(0.,8.)
    for TOL in TOLs:
            print cmd_multi.format(bayesian=True,
                               tag="new1AugD", TOL=TOL,
                               seed=np.random.randint(2**32-1))