def Corner3D():
"Calculate the variances for the corner point generation"
a = [ [ 2, 3.0/2 , 0.5 ],
[ 2, 2.0 , 0.5 ],
[ 2, 3.0/2 , 1 ]]
a=numarray.array(a)
f=PlainKolmogorovSpec2D
r=numarray.array( [ f(1), f(2.0**0.5), f(3.0**0.5)])
print r
return la.solve_linear_equations(a,r)
def Corner3DV2():
"Experimetal -- look for better initialisation"
a = [ [ 2, 3.0/2 , 0.5 ],
[ 2, 3.0 , 0.5 ],
[ 0, 3.0/2 , 1 ]]
a=numarray.array(a)
f=PlainKolmogorovSpec2D
r=numarray.array( [ f(1), f(2.0**0.5), f(3.0**0.5)])
print r
return la.solve_linear_equations(a,r)
def Corner3DV3():
"Test only"
# This corresponds to even term on edges and odd terms on face
# diagonals and body diagonals.
a = [ [ 1.0 , 3.0/2 , 1.0/2 ],
[ 3.0/2, 2.0 , 1.0/2 ],
[ 3.0/2, 3.0/2 , 1.0 ]]
a=numarray.array(a)
f=PlainKolmogorovSpec2D
r=numarray.array( [ f(1), f(2.0**0.5), f(3.0**0.5)])
print r
return la.solve_linear_equations(a,r)
[ 1,0,0,1,1,0 ],
[ 0,1,0,1,0,1 ],
[ 0,0,1,0,1,1 ],
[ 0,1,1,1,1,0 ],
[ 1,0,1,1,0,1 ] ]
rhs = [9, 1 ,1 ,1 , 4 , 4]
#import sys; sys.path.extend(["/import/appcs/bn204/p/bnprog-devel-main/"+x for x in ["bin", "lib"] ])
#import pybnlib
p2=numarray.array([3,-1, -1, -1])
# Symetric cube:
p3=numarray.array([7,-1, -1, -1, -1 ,-1, -1 , -1])
a,r=GenLinSystem(p3)
x=la.solve_linear_equations(a,r)
r2=la.mlab.innerproduct(a,x)
#ps=vtk.vtkPoints()
#ps.InsertPoint(0,1.1,1,1)
#et=vtk.vtkEdgeTable()
