multi_field.py 3.34 KB
Newer Older
Martin Reinecke's avatar
Martin Reinecke committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
from ..field import Field
import numpy as np
from .multi_domain import MultiDomain


class MultiField(object):
    def __init__(self, val):
        """
        Parameters
        ----------
        val : dict
        """
        self._val = val

    def __getitem__(self, key):
        return self._val[key]

    def keys(self):
        return self._val.keys()

    def items(self):
        return self._val.items()

    def values(self):
        return self._val.values()

    @property
    def domain(self):
        return MultiDomain({key: val.domain for key, val in self._val.items()})

    def _check_domain(self, other):
        if other.domain != self.domain:
            raise ValueError("domains are incompatible.")

    def vdot(self, x):
        result = 0.
        self._check_domain(x)
        for key, sub_field in self.items():
            result += sub_field.vdot(x[key])
        return result

    def lock(self):
        for v in self.values():
            v.lock()
        return self

    def copy(self):
        return MultiField({key: val.copy() for key, val in self.items()})

Martin Reinecke's avatar
Martin Reinecke committed
50
51
52
53
54
55
56
57
    @staticmethod
    def build_dtype(dtype, domain):
        if isinstance(dtype, dict):
            return dtype
        if dtype is None:
            dtype = np.float64
        return {key: dtype for key in domain.keys()}

Martin Reinecke's avatar
Martin Reinecke committed
58
59
    @staticmethod
    def zeros(domain, dtype=None):
Martin Reinecke's avatar
Martin Reinecke committed
60
61
        dtype = self.build_dtype(dtype, domain)
        return MultiField({key: Field.zeros(dom, dtype=dtype[key])
Martin Reinecke's avatar
cleanup    
Martin Reinecke committed
62
                           for key, dom in domain.items()})
Martin Reinecke's avatar
Martin Reinecke committed
63
64
65

    @staticmethod
    def ones(domain, dtype=None):
Martin Reinecke's avatar
Martin Reinecke committed
66
67
        dtype = self.build_dtype(dtype, domain)
        return MultiField({key: Field.ones(dom, dtype=dtype[key])
Martin Reinecke's avatar
cleanup    
Martin Reinecke committed
68
                           for key, dom in domain.items()})
Martin Reinecke's avatar
Martin Reinecke committed
69
70
71

    @staticmethod
    def empty(domain, dtype=None):
Martin Reinecke's avatar
Martin Reinecke committed
72
73
        dtype = self.build_dtype(dtype, domain)
        return MultiField({key: Field.empty(dom, dtype=dtype[key])
Martin Reinecke's avatar
cleanup    
Martin Reinecke committed
74
                           for key, dom in domain.items()})
Martin Reinecke's avatar
Martin Reinecke committed
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89

    def norm(self):
        """ Computes the L2-norm of the field values.

        Returns
        -------
        norm : float
            The L2-norm of the field values.
        """
        return np.sqrt(np.abs(self.vdot(x=self)))

    def __neg__(self):
        return MultiField({key: -val for key, val in self.items()})

    def conjugate(self):
Martin Reinecke's avatar
cleanup    
Martin Reinecke committed
90
91
        return MultiField({key: sub_field.conjugate()
                           for key, sub_field in self.items()})
Martin Reinecke's avatar
Martin Reinecke committed
92
93
94
95
96
97
98
99
100
101
102
103


for op in ["__add__", "__radd__", "__iadd__",
           "__sub__", "__rsub__", "__isub__",
           "__mul__", "__rmul__", "__imul__",
           "__div__", "__rdiv__", "__idiv__",
           "__truediv__", "__rtruediv__", "__itruediv__",
           "__floordiv__", "__rfloordiv__", "__ifloordiv__",
           "__pow__", "__rpow__", "__ipow__",
           "__lt__", "__le__", "__gt__", "__ge__", "__eq__", "__ne__"]:
    def func(op):
        def func2(self, other):
Martin Reinecke's avatar
cleanup    
Martin Reinecke committed
104
105
106
107
108
109
110
111
            if isinstance(other, MultiField):
                self._check_domain(other)
                result_val = {key: getattr(sub_field, op)(other[key])
                              for key, sub_field in self.items()}
            else:
                result_val = {key: getattr(val, op)(other)
                              for key, val in self.items()}
            return MultiField(result_val)
Martin Reinecke's avatar
Martin Reinecke committed
112
113
        return func2
    setattr(MultiField, op, func(op))