[python3] implements workaround for #171

src/pymor/operators/numpy.py

@@ -814,8 +814,14 @@ def _apply_inverse(matrix, V, options=None):
                     raise InversionError('bicgstab failed with error code {} (illegal input or breakdown)'.
                                          format(info))
     elif options['type'] == 'bicgstab_spilu':
-        ilu = spilu(matrix, drop_tol=options['spilu_drop_tol'], fill_factor=options['spilu_fill_factor'],
-                    drop_rule=options['spilu_drop_rule'], permc_spec=options['spilu_permc_spec'])
+        # workaround for https://github.com/pymor/pymor/issues/171
+        try:
+            ilu = spilu(matrix, drop_tol=options['spilu_drop_tol'], fill_factor=options['spilu_fill_factor'],
+                        drop_rule=options['spilu_drop_rule'], permc_spec=options['spilu_permc_spec'])
+        except TypeError as t:
+            logger = getLogger('pymor.operators.numpy._apply_inverse')
+            logger.error("ignoring drop_rule and permc_spec in ilu factorization")
+            ilu = spilu(matrix, drop_tol=options['spilu_drop_tol'], fill_factor=options['spilu_fill_factor'])
         precond = LinearOperator(matrix.shape, ilu.solve)
         for i, VV in enumerate(V):
             R[i], info = bicgstab(matrix, VV, tol=options['tol'], maxiter=options['maxiter'], M=precond)