def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
imp = try_imports(['from pysparse import jdsym, itsolvers, precon',
'from pysparse.eigen import jdsym;'
' from pysparse import itsolvers, precon'],
'cannot import pysparse eigensolvers!')
self.imp = imp
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
imp = try_imports(['from pysparse import jdsym, itsolvers, precon',
'from pysparse.eigen import jdsym;'
' from pysparse import itsolvers, precon'],
'cannot import pysparse eigensolvers!')
self.imp = imp
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
try:
import scipy.linalg.lapack as llapack
except ImportError:
import scipy.lib.lapack as llapack
self.llapack = llapack
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
import scipy.linalg as sla
self.sla = sla
aux = try_imports(['import scipy.sparse.linalg as ssla'],
'cannot import scipy sparse eigenvalue solvers!')
self.ssla = aux['ssla']
def __init__( self, conf, **kwargs ):
EigenvalueSolver.__init__( self, conf, **kwargs )
try:
import symeig
self.symeig = symeig.symeig
except:
self.symeig = None
output( 'cannot import symeig, required by %s solver' % self.name )
raise
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
try:
import scipy.linalg.lapack as llapack
except ImportError:
import scipy.lib.lapack as llapack
self.llapack = llapack
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
import scipy.linalg as sla
self.sla = sla
aux = try_imports(['import scipy.sparse.linalg as ssla'],
'cannot import scipy sparse eigenvalue solvers!')
self.ssla = aux['ssla']
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
imp = try_imports(
[
"from pysparse import jdsym, itsolvers, precon",
"from pysparse.eigen import jdsym;" " from pysparse import itsolvers, precon",
],
"cannot import pysparse eigensolvers!",
)
self.imp = imp
def __init__(self, conf, comm=None, context=None, **kwargs):
if comm is None:
init_slepc_args()
from petsc4py import PETSc as petsc
from slepc4py import SLEPc as slepc
EigenvalueSolver.__init__(self,
conf,
petsc=petsc,
slepc=slepc,
comm=comm,
context=context,
**kwargs)
def process_conf(conf, kwargs):
"""
Missing items are set to default values.
Example configuration, all items::
solver_2 = {
'name' : 'lobpcg',
'kind' : 'eig.scipy_lobpcg',
'i_max' : 20,
'n_eigs' : 5,
'eps_a' : None,
'largest' : True,
'precond' : None,
'verbosity' : 0,
}
"""
get = make_get_conf(conf, kwargs)
common = EigenvalueSolver.process_conf(conf)
return Struct(i_max=get('i_max', 20),
n_eigs=get('n_eigs', None),
eps_a=get('eps_a', None),
largest=get('largest', True),
precond=get('precond', None),
verbosity=get('verbosity', 0)) + common
def process_conf( conf ):
"""
Missing items are set to default values.
Example configuration, all items:
solver_2 = {
'name' : 'eigen1',
'kind' : 'eig.pysparse',
'i_max' : 150,
'eps_a' : 1e-5,
'tau' : -10.0,
'method' : 'qmrs',
'verbosity' : 0,
'strategy' : 1,
}
"""
get = conf.get_default_attr
i_max = get( 'i_max', 100 )
eps_a = get( 'eps_a', 1e-5 )
tau = get( 'tau', 0.0 )
method = get( 'method', 'qmrs' )
verbosity = get( 'verbosity', 0 )
strategy = get( 'strategy', 1 )
common = EigenvalueSolver.process_conf( conf )
return Struct( **locals() ) + common
def process_conf(conf, kwargs):
"""
Missing items are set to default values.
Example configuration, all items::
solver_2 = {
'name' : 'eigen1',
'kind' : 'eig.pysparse',
'i_max' : 150,
'eps_a' : 1e-5,
'tau' : -10.0,
'method' : 'qmrs',
'verbosity' : 0,
'strategy' : 1,
}
"""
get = make_get_conf(conf, kwargs)
common = EigenvalueSolver.process_conf(conf)
return Struct(i_max=get('i_max', 100),
n_eigs=get('n_eigs', 5),
eps_a=get('eps_a', 1e-5),
tau=get('tau', 0.0),
method=get('method', 'qmrs'),
verbosity=get('verbosity', 0),
strategy=get('strategy', 1)) + common
def process_conf(conf, kwargs):
"""
Missing items are set to default values.
Example configuration, all items::
solver_2 = {
'name' : 'lobpcg',
'kind' : 'eig.scipy_lobpcg',
'i_max' : 20,
'n_eigs' : 5,
'eps_a' : None,
'largest' : True,
'precond' : None,
'verbosity' : 0,
}
"""
get = make_get_conf(conf, kwargs)
common = EigenvalueSolver.process_conf(conf)
return Struct(i_max=get('i_max', 20),
n_eigs=get('n_eigs', None),
eps_a=get('eps_a', None),
largest=get('largest', True),
precond=get('precond', None),
verbosity=get('verbosity', 0)) + common
def process_conf(conf, kwargs):
"""
Missing items are set to default values.
Example configuration, all items::
solver_2 = {
'name' : 'eigen1',
'kind' : 'eig.pysparse',
'i_max' : 150,
'eps_a' : 1e-5,
'tau' : -10.0,
'method' : 'qmrs',
'verbosity' : 0,
'strategy' : 1,
}
"""
get = make_get_conf(conf, kwargs)
common = EigenvalueSolver.process_conf(conf)
return Struct(i_max=get('i_max', 100),
eps_a=get('eps_a', 1e-5),
tau=get('tau', 0.0),
method=get('method', 'qmrs'),
verbosity=get('verbosity', 0),
strategy=get('strategy', 1)) + common
def process_conf( conf ):
"""
Missing items are set to default values.
Example configuration, all items::
solver_2 = {
'name' : 'lobpcg',
'kind' : 'eig.scipy_lobpcg',
'i_max' : 20,
'n_eigs' : 5,
'eps_a' : None,
'largest' : True,
'precond' : None,
'verbosity' : 0,
}
"""
get = conf.get_default_attr
i_max = get( 'i_max', 20 )
n_eigs = get( 'n_eigs', None )
eps_a = get( 'eps_a', None )
largest = get( 'largest', True )
precond = get( 'precond', None )
verbosity = get( 'verbosity', 0 )
common = EigenvalueSolver.process_conf( conf )
return Struct( **locals() ) + common
def process_conf(conf, kwargs):
"""
Missing items are set to default values.
Example configuration, all items::
solver_20 = {
'name' : 'eigen2',
'kind' : 'eig.sgscipy',
'force_n_eigs' : True,
}
"""
get = make_get_conf(conf, kwargs)
common = EigenvalueSolver.process_conf(conf)
return Struct(force_n_eigs=get('force_n_eigs', False)) + common
def process_conf(conf, kwargs):
"""
Missing items are set to default values.
Example configuration, all items::
solver_20 = {
'name' : 'eigen2',
'kind' : 'eig.sgscipy',
'force_n_eigs' : True,
}
"""
get = make_get_conf(conf, kwargs)
common = EigenvalueSolver.process_conf(conf)
return Struct(force_n_eigs=get('force_n_eigs', False)) + common
def __init__(self, conf, comm=None, context=None, **kwargs):
import matlab.engine as me
EigenvalueSolver.__init__(self, conf, me=me, context=context, **kwargs)
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
from scipy.sparse.linalg.eigen import lobpcg
self.lobpcg = lobpcg
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
import scipy.lib.lapack as llapack
self.llapack = llapack
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
from scipy.sparse.linalg.eigen import lobpcg
self.lobpcg = lobpcg
def __init__(self, conf, **kwargs):
EigenvalueSolver.__init__(self, conf, **kwargs)
