Commit 56860b32 authored by Matthijs's avatar Matthijs
Browse files

change should not be complex dtype, which would happen if u=complex. Fixed...

change should not be complex dtype, which would happen if u=complex. Fixed this by not specifying dtype
parent ecc617cf
......@@ -116,7 +116,7 @@ class SimpleAlgorithm:
p = u.shape[2]
q = u.shape[3]
change = zeros(maxiter + 1, dtype=u.dtype)
change = zeros(maxiter + 1)
change[0] = 999
######################################
......
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