aboutsummaryrefslogtreecommitdiff
path: root/boltzgen/kernel/generator.py
blob: 8d0de37c7e4eb3fd0f6ee14eaf95b03464c13ee4 (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
import sympy

from mako.template import Template
from pathlib import Path

import kernel.target.cl
import kernel.target.cpp

class Generator:
    def __init__(self, descriptor, moments, collision, boundary = ''):
        self.descriptor = descriptor
        self.moments    = moments
        self.collision  = collision
        self.boundary   = boundary

    def kernel(self, target, precision, layout, geometry, extras = []):
        template_path = Path(__file__).parent/('template/basic.' + target + '.mako')
        if not template_path.exists():
            raise Exception("Target '%s' not supported" % target)

        layout_impl = eval("kernel.target.%s.%s" % (target, layout))
        if layout_impl is None:
            raise Exception("Target '%s' doesn't support layout '%s'" % (target, layout))
        else:
            layout_impl = layout_impl(self.descriptor, geometry)

        if geometry.dimension() != self.descriptor.d:
            raise Exception('Geometry dimension must match descriptor dimension')

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

            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 = {
                'single': 'float',
                'double': 'double'
            }.get(precision),

            boundary_src = Template(self.boundary).render(
                descriptor = self.descriptor,
                geometry   = geometry,
                float_type = precision
            ),

            extras = extras
        )