Commit 29803a51 authored by Ruestig, Julian (jruestig)'s avatar Ruestig, Julian (jruestig) 📡

Clearer naming in correlated_field mf

parent f8796588
Pipeline #45515 passed with stages
in 9 minutes and 16 seconds
......@@ -119,15 +119,15 @@ def MfCorrelatedField(target, amplitudes, name='xi'):
ht = ht2 @ ht1
psp = [aa.target[0] for aa in amplitudes]
pd0 = PowerDistributor(hsp, psp[0], 0)
pd1 = PowerDistributor(pd0.domain, psp[1], 1)
pd0 = PowerDistributor(hsp, power_space=psp[0], space=0)
pd1 = PowerDistributor(pd0.domain, power_space=psp[1], space=1)
pd = pd0 @ pd1
domain0 = ContractionOperator(pd.domain, 1).adjoint
domain1 = ContractionOperator(pd.domain, 0).adjoint
d = [domain0, domain1]
blowup0 = ContractionOperator(pd.domain, 1).adjoint
blowup1 = ContractionOperator(pd.domain, 0).adjoint
d = [blowup0, blowup1]
a = [dd @ amplitudes[ii] for ii, dd in enumerate(d)]
a = [BUp @ amplitudes[ii] for ii, BUp in enumerate(d)]
a = reduce(mul, a)
A = pd @ a
# For `vol` see comment in `CorrelatedField`
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment