Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
ift
NIFTy
Commits
505b3ab6
Commit
505b3ab6
authored
May 19, 2020
by
Martin Reinecke
Browse files
Merge branch 'integration_operator' of gitlab.mpcdf.mpg.de:ift/nifty into integration_operator
parents
40a9d52f
cba6ac90
Pipeline
#75197
passed with stages
in 25 minutes and 10 seconds
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
test/test_operators/test_integration.py
View file @
505b3ab6
...
...
@@ -24,25 +24,25 @@ from ..common import setup_function, teardown_function
def
test_integration_operator
():
x1
=
ift
.
RGSpace
((
9
,),
distances
=
2.
)
# FIXME: it works for distances = 1
#
x2 = ift.RGSpace((2, 12), distances=(0.3,))
x1
=
ift
.
RGSpace
((
9
,),
distances
=
2.
)
x2
=
ift
.
RGSpace
((
2
,
12
),
distances
=
(
0.3
,))
dom1
=
ift
.
makeDomain
(
x1
)
#
dom2 = ift.makeDomain((x1, x2))
dom2
=
ift
.
makeDomain
((
x1
,
x2
))
f1
=
ift
.
from_random
(
'normal'
,
dom1
)
#
f2 = ift.from_random('normal', dom2)
f2
=
ift
.
from_random
(
'normal'
,
dom2
)
op1
=
ift
.
ScalingOperator
(
dom1
,
1
)
#
op2 = ift.ScalingOperator(dom2, 1)
op2
=
ift
.
ScalingOperator
(
dom2
,
1
)
res1
=
f1
.
integrate
()
res2
=
op1
.
integrate
()(
f1
)
assert_allclose
(
res1
.
val
,
res2
.
val
)
#
res3 = f2.integrate()
#
res4 = op2.integrate()(f2)
#
assert_allclose(res3.val, res4.val)
#
res5 = f2.integrate(spaces=1)
#
res6 = op2.integrate(spaces=1)(f)
#
assert_allclose(res5.val, res6.val)
#
for op in [op1, op2]:
#
ift.extra.consistency_check(op, domain_dtype=np.float64,
#
target_dtype=np.float64)
#
ift.extra.consistency_check(op, domain_dtype=np.complex128,
#
target_dtype=np.comple
y
128)
res3
=
f2
.
integrate
()
res4
=
op2
.
integrate
()(
f2
)
assert_allclose
(
res3
.
val
,
res4
.
val
)
res5
=
f2
.
integrate
(
spaces
=
1
)
res6
=
op2
.
integrate
(
spaces
=
1
)(
f
2
)
assert_allclose
(
res5
.
val
,
res6
.
val
)
for
op
in
[
op1
,
op2
]:
ift
.
extra
.
consistency_check
(
op
,
domain_dtype
=
np
.
float64
,
target_dtype
=
np
.
float64
)
ift
.
extra
.
consistency_check
(
op
,
domain_dtype
=
np
.
complex128
,
target_dtype
=
np
.
comple
x
128
)
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment