Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F60985185
pade_utils.py
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Fri, May 3, 19:31
Size
1 KB
Mime Type
text/x-python
Expires
Sun, May 5, 19:31 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
17449343
Attached To
R6746 RationalROMPy
pade_utils.py
View Options
# Copyright (C) 2018 by the RROMPy authors
#
# This file is part of RROMPy.
#
# RROMPy is free software: you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# RROMPy is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with RROMPy. If not, see <http://www.gnu.org/licenses/>.
#
import
numpy
as
np
from
rrompy.utilities.base.types
import
Np1D
,
Np2D
,
Tuple
from
rrompy.utilities.warning_manager
import
warn
machineEps
=
np
.
finfo
(
float
)
.
eps
def
checkRobustTolerance
(
ev
:
Np1D
,
E
:
int
,
tol
:
float
)
->
dict
:
N
=
len
(
ev
)
-
1
ts
=
tol
*
np
.
linalg
.
norm
(
ev
)
diff
=
N
-
np
.
sum
(
np
.
abs
(
ev
)
>=
ts
)
if
diff
<=
0
:
return
{}
Enew
=
E
-
diff
Nnew
=
min
(
N
,
Enew
)
if
Nnew
==
N
:
strN
=
""
else
:
strN
=
"N from {} to {} and "
.
format
(
N
,
Nnew
)
warn
((
"Smallest {} eigenvalues below tolerance.
\n
Reducing {}E from {} to "
"{}."
)
.
format
(
diff
+
1
,
strN
,
E
,
Enew
))
newParameters
=
{
"N"
:
Nnew
,
"E"
:
Enew
}
return
newParameters
def
checkPolyfitRank
(
xs
:
Np1D
,
ys
:
Np2D
,
E
:
int
,
rcond
:
float
,
w
:
Np1D
=
None
)
->
Tuple
[
bool
,
Np2D
]:
G
=
np
.
polyfit
(
xs
,
ys
,
deg
=
E
,
w
=
w
,
rcond
=
rcond
,
full
=
True
)
if
G
[
2
]
<
E
+
1
:
return
False
,
G
[
2
]
-
1
return
True
,
G
[
0
]
def
checkFitRank
(
xs
:
Np1D
,
ys
:
Np2D
,
deg
:
int
,
rcond
:
float
,
fit
:
callable
,
w
:
Np1D
=
None
)
->
Tuple
[
bool
,
Np2D
]:
G
=
fit
(
xs
,
ys
,
deg
,
w
=
w
,
rcond
=
rcond
,
full
=
True
)
if
G
[
1
][
1
]
<
deg
+
1
:
return
False
,
G
[
1
][
1
]
-
1
return
True
,
G
[
0
]
Event Timeline
Log In to Comment