Skip to content

Commit a7ea5e0

Browse files
Tiny refactor (#570)
1 parent 66a3d6d commit a7ea5e0

File tree

1 file changed

+4
-6
lines changed

1 file changed

+4
-6
lines changed

pySDC/helpers/spectral_helper.py

Lines changed: 4 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -312,14 +312,13 @@ def get_conv(self, name, N=None):
312312
'''
313313
N = N if N else self.N
314314
sp = self.sparse_lib
315-
xp = self.xp
316315

317316
def get_forward_conv(name):
318317
if name == 'T2U':
319-
mat = (sp.eye(N) - sp.diags(xp.ones(N - 2), offsets=+2)).tocsc() / 2.0
318+
mat = (sp.eye(N) - sp.eye(N, k=2)).tocsc() / 2.0
320319
mat[:, 0] *= 2
321320
elif name == 'D2T':
322-
mat = sp.eye(N) - sp.diags(xp.ones(N - 2), offsets=+2)
321+
mat = sp.eye(N) - sp.eye(N, k=2)
323322
elif name[0] == name[-1]:
324323
mat = self.sparse_lib.eye(self.N)
325324
else:
@@ -585,9 +584,8 @@ def get_Dirichlet_recombination_matrix(self):
585584
'''
586585
N = self.N
587586
sp = self.sparse_lib
588-
xp = self.xp
589587

590-
return sp.eye(N) - sp.diags(xp.ones(N - 2), offsets=+2)
588+
return sp.eye(N) - sp.eye(N, k=2)
591589

592590

593591
class UltrasphericalHelper(ChebychevHelper):
@@ -632,7 +630,7 @@ def get_S(self, lmbda):
632630

633631
if lmbda == 0:
634632
sp = scipy.sparse
635-
mat = ((sp.eye(N) - sp.diags(np.ones(N - 2), offsets=+2)) / 2.0).tolil()
633+
mat = ((sp.eye(N) - sp.eye(N, k=2)) / 2.0).tolil()
636634
mat[:, 0] *= 2
637635
else:
638636
sp = self.sparse_lib

0 commit comments

Comments
 (0)