compute propagator for each block of block-diagonal system matrix

odetoolbox/system_of_shapes.py

@@ -39,7 +39,7 @@ def get_block_diagonal_blocks(A):
     assert A.shape[0] == A.shape[1], "matrix A should be square"
 
     N = A.shape[0]
-    
+
     # A_mirrored = (A + A.T).applyfunc(_is_zero)   # make the matrix symmetric so we only have to check one triangle
     A_mirrored = (A + A.T) != 0   # make the matrix symmetric so we only have to check one triangle