コード例 #1
0
ファイル: fiber_cb_8p.py プロジェクト: simvisage/spirrid
                          x = np.linspace(-20.1, 20.5, 100),
                          Lr = np.linspace(0.1, 20.0, 50)
                          ),
             theta_vars = dict(tau = RV('uniform', 0.7, 1.0),
                          l = RV('uniform', 5.0, 10.0),
                          D_f = 26e-3,
                          E_f = 72e3,
                          theta = 0.0,
                          xi = RV('weibull_min', scale = 0.017, shape = 8, n_int = 10),
                          phi = 1.0,
                          Ll = 50.0,
 #                              Lr = 1.0
                          ),
             n_int = 5)
 
 e_arr = make_ogrid(s.evar_lst)
 n_e_arr = [ e / np.max(np.fabs(e)) for e in e_arr ]
 
 max_mu_q = np.max(np.fabs(s.mu_q_arr))
 n_mu_q_arr = s.mu_q_arr / max_mu_q
 n_std_q_arr = np.sqrt(s.var_q_arr) / max_mu_q
 
 #===========================================================================
 # Prepare plotting 
 #===========================================================================
 tdir = tempfile.mkdtemp()
 n_img = n_mu_q_arr.shape[0]
 fnames = [os.path.join(tdir, 'x%02d.jpg' % i) for i in range(n_img) ]
 
 f = m.figure(1, size = (1000, 500), fgcolor = (0, 0, 0),
              bgcolor = (1., 1., 1.))
コード例 #2
0
ファイル: ogrid_test.py プロジェクト: kelidas/scratch
# This software is provided without warranty under the terms of the BSD
# license included in the Spirrid top directory "licence.txt" and may be
# redistributed only under the conditions described in the aforementioned
# license.
#
# Thanks for using Simvisage open source!
#
#-------------------------------------------------------------------------------

from spirrid import make_ogrid
import numpy as np

evar_names = ['a', 'b', 'c']
evar_lst = [np.arange(1, 5, 1), np.arange(5, 10, 1), np.arange(10, 15, 1)]

e_ogrid = make_ogrid(evar_lst)

def f(a, b, c,):
    return a * b * c


def mu_q(*e):
    eargs = dict(zip(evar_names, e))
    res = f(**eargs)
    return res

otypes = [ float for i in range(len(evar_lst))]
mu_q_vec = np.vectorize(mu_q, otypes=[float])

print mu_q_vec(*e_ogrid)