Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F90042997
funnel_output_engine.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
Mon, Oct 28, 18:58
Size
1 KB
Mime Type
text/x-python
Expires
Wed, Oct 30, 18:58 (2 d)
Engine
blob
Format
Raw Data
Handle
21998263
Attached To
R6746 RationalROMPy
funnel_output_engine.py
View Options
import
numpy
as
np
import
fenics
as
fen
import
mshr
from
rrompy.hfengines.linear_problem
import
ScatteringProblemEngine
from
rrompy.solver.fenics
import
fenics2Sparse
class
FunnelOutputEngine
(
ScatteringProblemEngine
):
def
__init__
(
self
,
k0
:
float
,
n
:
int
):
super
()
.
__init__
(
mu0
=
[
k0
])
mesh
=
mshr
.
generate_mesh
(
mshr
.
Circle
(
fen
.
Point
(
0.
,
0.
),
5.
)
-
mshr
.
Rectangle
(
fen
.
Point
(
-
5.
,
-
5.
),
fen
.
Point
(
-
1.
,
5.
))
-
mshr
.
Rectangle
(
fen
.
Point
(
-
5.
,
-
5.
),
fen
.
Point
(
5.
,
0.
))
-
mshr
.
Rectangle
(
fen
.
Point
(
-
1.
,
1.
),
fen
.
Point
(
0.
,
5.
)),
n
)
self
.
V
=
fen
.
FunctionSpace
(
mesh
,
"P"
,
1
)
eps
=
1e-4
self
.
DirichletBoundary
=
(
lambda
x
,
on_boundary
:
on_boundary
and
x
[
0
]
<
-
1.
+
eps
)
self
.
NeumannBoundary
=
(
lambda
x
,
on_boundary
:
on_boundary
and
x
[
0
]
>=
-
1.
+
eps
and
(
x
[
0
]
<
eps
or
x
[
1
]
<=
eps
))
self
.
RobinBoundary
=
"REST"
y
=
fen
.
SpatialCoordinate
(
mesh
)[
1
]
self
.
DirichletDatum
=
fen
.
sin
(
.
5
*
np
.
pi
*
y
)
/
(
1.
+
y
**
2.
)
self
.
autoSetDS
()
l2R0
=
fen
.
inner
(
self
.
u
,
self
.
v
)
*
self
.
ds
(
1
)
L2R0
=
fenics2Sparse
(
l2R0
,
{},
None
,
-
1
)
bcR
=
np
.
where
(
np
.
abs
(
L2R0
.
dot
(
np
.
ones
(
L2R0
.
shape
[
1
])))
>
1e-10
)[
0
]
bcR
=
bcR
[
np
.
argsort
(
self
.
V
.
tabulate_dof_coordinates
()[
bcR
,
1
])]
self
.
_C
=
np
.
zeros
((
len
(
bcR
),
L2R0
.
shape
[
1
]))
self
.
_C
[
np
.
arange
(
len
(
bcR
)),
bcR
]
=
1.
self
.
outputNormMatrix
=
L2R0
[
bcR
][:,
bcR
]
Event Timeline
Log In to Comment