Exemple #1
0
def interpolate(C, a):

    #####################################
    # interpolation
    scalar1 = regression.convert2(C, a)
    scalar1.name = "source term for" + a.name + "and" + b.name
    scalar1.show2()
    scalar1.show3()
    vect1   = regression.convert(C, a)
    vect1.show()
    tic()
    a1  = sl.interpolate2(a, vect1)
    toc()
    a1.show()
    return a1, vect1, scalar1
    print "Rsquared:", Rsquared
    print "shiiba coeffs:", C
    print "upWindCorrectionTerm: (", len(upWindCorrectionTerm), ")", upWindCorrectionTerm,

##################################
# Shiiba regression results
print "\n\n\n=============== Shiiba regression results ===================\n\n"

Rsquared
phi0
phi1
C
Rsquared

#####################################
# interpolation

vect1 = regression.convert(C, a)
vect1.show()
tic()
a1 = sl.interpolate2(a, vect1)
toc()
a1.show()

diff = b - a1
diff.show()
print diff.matrix.max()
print diff.matrix.min()
print diff.matrix.mean()
print diff.matrix.var()