Delete comment and add dtype keyword
@@ -308,9 +308,8 @@ class Surface():
# If we have complex dielectric constants we need to initialize Precon with
# complex type else it'll be float.
- #if type(self.E_hat) == complex:
if numpy.iscomplexobj(self.E_hat):
- self.Precond = numpy.zeros((4, self.N), type(self.E_hat))
+ self.Precond = numpy.zeros((4, self.N), dtype=type(self.E_hat))
else:
self.Precond = numpy.zeros((4, self.N))
# Stores the inverse of the block diagonal (also a tridiag matrix)
