aboutsummaryrefslogtreecommitdiff
path: root/boltzgen/kernel/generator.py
blob: 92ceb7228e07f7ab23ec744038ca3014860714d5 (plain)
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
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
import sympy

from mako.template import Template
from pathlib import Path

from . import memory

class Generator:
    def __init__(self, descriptor, moments, collision, target, precision, index, layout):
        self.descriptor = descriptor
        self.moments    = moments
        self.collision  = collision
        self.target     = target
        self.float_type = eval("memory.precision.%s" % target).get_float_type(precision)

        try:
            self.index_impl = eval("memory.index.%s" % index)
        except AttributeError:
            raise Exception("There is no cell indexing scheme '%s'" % index) from None

        try:
            self.layout_impl = eval("memory.layout.%s" % layout)
        except AttributeError:
            raise Exception("There is no layout '%s'" % layout) from None

    def instantiate(self, template, geometry, extras = []):
        template_path = Path(__file__).parent/("template/%s.%s.mako" % (template, self.target))
        if not template_path.exists():
            raise Exception("Target '%s' doesn't provide '%s'" % (self.target, template))

        return Template(filename = str(template_path)).render(
            descriptor = self.descriptor,
            geometry   = geometry,
            index      = self.index_impl(geometry),
            layout     = self.layout_impl(self.descriptor, self.index_impl, geometry),

            moments_subexpr    = self.moments[0],
            moments_assignment = self.moments[1],
            collision_subexpr    = self.collision[0],
            collision_assignment = self.collision[1],
            ccode = sympy.ccode,

            float_type = self.float_type,

            extras = extras
        )

    def kernel(self, geometry, functions, extras = []):
        if geometry.dimension() != self.descriptor.d:
            raise Exception('Geometry dimension must match descriptor dimension')

        return "\n".join(map(lambda f: self.instantiate(f, geometry, extras), functions))

    def custom(self, geometry, source):
        return Template(text = source).render(
            descriptor = self.descriptor,
            geometry   = geometry,
            index      = self.index_impl(geometry),
            layout     = self.layout_impl(self.descriptor, self.index_impl, geometry),

            moments_subexpr    = self.moments[0],
            moments_assignment = self.moments[1],
            collision_subexpr    = self.collision[0],
            collision_assignment = self.collision[1],
            ccode = sympy.ccode,

            float_type = self.float_type,
        )