Commit d89aa7b3 authored by Philipp Arras's avatar Philipp Arras

Fixups

parent 8907d503
......@@ -226,8 +226,10 @@ class FinalAmplitude:
hspace, zeroind) @ azm
from ..operators.harmonic_operators import HarmonicTransformOperator
ht = HarmonicTransformOperator(hspace, space=0)
pd = PowerDistributor(hspace, self._amplitudes[0].target[0], 0)
for i in range(1, len(self._amplitudes)):
ht = HarmonicTransformOperator(ht.target, space=i) @ ht
pd = pd @ PowerDistributor(
pd.domain, self._amplitudes[i].target[0], space=i)
......@@ -240,5 +242,4 @@ class FinalAmplitude:
(i + 1):]).adjoint(self._amplitudes[i]))
A = pd @ a
ht = HarmonicTransformOperator(hspace)
return ht(A*ducktape(hspace, None, prefix + 'xi'))
......@@ -18,9 +18,7 @@
import numpy as np
from ..domain_tuple import DomainTuple
from ..domains.unstructured_domain import UnstructuredDomain
from ..field import Field
from ..sugar import makeDomain
from .linear_operator import LinearOperator
......@@ -36,7 +34,7 @@ class ValueInserter(LinearOperator):
"""
def __init__(self, target, index):
self._domain = makeDomain(UnstructuredDomain(1))
self._domain = DomainTuple.scalar_domain()
self._target = DomainTuple.make(target)
index = tuple(index)
if not all([
......
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