num_dims = 8
virus_model = VirusInfection()
r = virus_model.run( virus_model.nominal_z )
func_domain = TensorProductDomain( num_dims,
numpy.hstack( ( virus_model.nominal_z.reshape(num_dims,1)*0.9,
virus_model.nominal_z.reshape(num_dims,1)*1.1 ) ) )
plot_2d_interaction_surfaces( virus_model, virus_model.nominal_z,
func_domain, active_dims = None,
num_test_pts_1d = 30 )
num_dims = 6
func_domain = TensorProductDomain( num_dims,
numpy.hstack( ( reaction_model.nominal_z.reshape(num_dims,1)*0.9,
reaction_model.nominal_z.reshape(num_dims,1)*1.1 ) ) )
plot_2d_interaction_surfaces( reaction_model, reaction_model.nominal_z,
func_domain, active_dims = None,
num_test_pts_1d = 30 )
"""
num_dims = 6
func_domain = TensorProductDomain( num_dims, [[0.08,0.12],[0.03,0.04],
[0.08,0.12],[0.8,1.2],
[0.45,0.55],[-0.05,0.05]] )
plot_2d_interaction_surfaces( oscillator_model, oscillator_model.nominal_z,
func_domain, active_dims = None,
num_test_pts_1d = 30 )
pylab.show()
#pylab.show()
# convergence analysis
"""
z = numpy.random.uniform( 0., 1., 10 )
q_old = 0.
for nx in numpy.arange( 1001, 40001, 1000 ):
heat_model.nx = nx
q = heat_model.evaluate( z )
print nx, abs(q-q_old)
q_old = q
"""
# exact solution for k = 1.0
exact_u = x * ( 1.0 - x ) * numpy.exp ( x )
print numpy.absolute( u - exact_u )
#pylab.plot(x,exact_u,'r')
from utilities.visualisation import plot_2d_interaction_surfaces
from utilities.tensor_product_domain import TensorProductDomain
num_dims = 10
r = heat_model.run( heat_model.nominal_z )
func_domain = TensorProductDomain( num_dims,
numpy.hstack( ( heat_model.nominal_z.reshape(num_dims,1)-1.,
heat_model.nominal_z.reshape(num_dims,1)+1. ) ) )
plot_2d_interaction_surfaces( heat_model, heat_model.nominal_z,
func_domain, active_dims = None,
num_test_pts_1d = 10 )
