Commit a684bf2f authored by Lars Stietz's avatar Lars Stietz
Browse files

dtype error solved

parent 76d52ca1
...@@ -15,7 +15,7 @@ coll = CollGaussLobatto(num_nodes,t_interval[0],t_interval[1]) ...@@ -15,7 +15,7 @@ coll = CollGaussLobatto(num_nodes,t_interval[0],t_interval[1])
# initialize Vector # initialize Vector
u_0 = 1 u_0 = 1
l = -5j l = -5j
u = u_0*np.ones(num_nodes+2, dtype = np.complex64) u = u_0*np.ones(num_nodes+2, dtype = complex)
...@@ -96,7 +96,7 @@ if solving == 1: ...@@ -96,7 +96,7 @@ if solving == 1:
# SDC straight forward approach # SDC straight forward approach
k_max = 50 k_max = 100
res_2 = 2 res_2 = 2
res_it = np.empty(0) res_it = np.empty(0)
res_itM = np.empty(0) res_itM = np.empty(0)
...@@ -109,13 +109,13 @@ if solving == 1: ...@@ -109,13 +109,13 @@ if solving == 1:
# LU Q_slow # LU Q_slow
#_,Q_slow,_ = la.lu(np.transpose(coll.Qmat[1:,1:])) #_,Q_fast,_ = la.lu(np.transpose(coll.Qmat[1:,1:]))
#Q_slow = np.transpose(Q_slow) #Q_fast = np.transpose(Q_slow)
# standard Q_slow # standard Q_fast
Q_slow = asmFastMat(coll.delta_m) Q_fast = asmFastMat(coll.delta_m)
u_Mat = sweepMat(u_Mat,coll.delta_m,coll.Qmat,Q_slow,f) u_Mat = sweepMat(u_Mat,coll.delta_m,coll.Qmat,Q_fast,f)
r = res(u,coll.Qmat,u_0) r = res(u,coll.Qmat,u_0)
r_Mat = res(u_Mat,coll.Qmat,u_0) r_Mat = res(u_Mat,coll.Qmat,u_0)
res_2 = np.linalg.norm(r) res_2 = np.linalg.norm(r)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment