epicf/ParticleSourceBox.py at master · dumbman/epicf · GitHub

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
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
from Vec3d import Vec3d
from ParticleSource import ParticleSource


class ParticleSourceBox(ParticleSource):

    def __init__(self):
        super().__init__()
        self.xleft = None
        self.xright = None
        self.ytop = None
        self.ybottom = None
        self.znear = None
        self.zfar = None


    @classmethod
    def init_from_config(cls, conf, this_source_config_part, sec_name):
        new_obj = cls()
        new_obj.read_particles_and_source_pars_from_config(
            conf, this_source_config_part, sec_name)
        new_obj.geometry_type = "box"
        new_obj.check_correctness_of_box_config_fields(
            conf, this_source_config_part)
        new_obj.set_box_parameters_from_config(this_source_config_part)
        new_obj.generate_initial_particles()
        return new_obj


    @classmethod
    def init_from_h5(cls, h5_source_group):
        new_obj = cls()
        new_obj.read_particles_and_source_pars_from_h5(h5_source_group)
        new_obj.geometry_type = "box"
        new_obj.read_hdf5_box_parameters(h5_source_group)
        return new_obj


    def check_correctness_of_box_config_fields(self, conf, this_source_config_part):
        self.x_right_ge_zero(conf, this_source_config_part)
        self.x_right_le_particle_source_x_left(conf, this_source_config_part)
        self.x_left_le_grid_x_size(conf, this_source_config_part)
        self.y_bottom_ge_zero(conf, this_source_config_part)
        self.y_bottom_le_particle_source_y_top(conf, this_source_config_part)
        self.y_top_le_grid_y_size(conf, this_source_config_part)
        self.z_near_ge_zero(conf, this_source_config_part)
        self.z_near_le_particle_source_z_far(conf, this_source_config_part)
        self.z_far_le_grid_z_size(conf, this_source_config_part)


    def set_box_parameters_from_config(self, this_source_config_part):
        self.xleft = this_source_config_part.getfloat("box_x_left")
        self.xright = this_source_config_part.getfloat("box_x_right")
        self.ytop = this_source_config_part.getfloat("box_y_top")
        self.ybottom = this_source_config_part.getfloat("box_y_bottom")
        self.znear = this_source_config_part.getfloat("box_z_near")
        self.zfar = this_source_config_part.getfloat("box_z_far")


    def read_hdf5_box_parameters(self, this_source_h5_group):
        self.xleft = this_source_h5_group.attrs["box_x_left"]
        self.xright = this_source_h5_group.attrs["box_x_right"]
        self.ytop = this_source_h5_group.attrs["box_y_top"]
        self.ybottom = this_source_h5_group.attrs["box_y_bottom"]
        self.zfar = this_source_h5_group.attrs["box_z_far"]
        self.znear = this_source_h5_group.attrs["box_z_near"]


    def x_right_ge_zero(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_x_right") < 0:
            raise ValueError("Expect box_x_right >= 0")


    def x_right_le_particle_source_x_left(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_x_right") > \
           this_source_config_part.getfloat("box_x_left"):
            raise ValueError("Expect box_x_right <= box_x_left")


    def x_left_le_grid_x_size(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_x_left") > \
           conf["SpatialMesh"].getfloat("grid_x_size"):
            raise ValueError("Expect box_x_left <= grid_x_size")


    def y_bottom_ge_zero(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_y_bottom") < 0:
            raise ValueError("Expect box_y_bottom >= 0")


    def y_bottom_le_particle_source_y_top(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_y_bottom") > \
           this_source_config_part.getfloat("box_y_top"):
            raise ValueError("Expect box_y_bottom <= box_y_top")


    def y_top_le_grid_y_size(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_y_top") > \
           conf["SpatialMesh"].getfloat("grid_y_size"):
            raise ValueError("Expect box_y_top <= grid_y_size")


    def z_near_ge_zero(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_z_near") < 0:
            raise ValueError("Expect box_z_near >= 0")


    def z_near_le_particle_source_z_far(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_z_near") > \
           this_source_config_part.getfloat("box_z_far"):
            raise ValueError("Expect box_z_near <= box_z_far")


    def z_far_le_grid_z_size(self, conf, this_source_config_part):
        if this_source_config_part.getfloat("box_z_far") > \
           conf["SpatialMesh"].getfloat("grid_z_size"):
            raise ValueError("Expect box_z_far <= grid_z_size")


    def write_hdf5_source_parameters(self, this_source_h5_group):
        super().write_hdf5_source_parameters(this_source_h5_group)
        this_source_h5_group.attrs.create("box_x_left", self.xleft)
        this_source_h5_group.attrs.create("box_x_right", self.xright)
        this_source_h5_group.attrs.create("box_y_top", self.ytop)
        this_source_h5_group.attrs.create("box_y_bottom", self.ybottom)
        this_source_h5_group.attrs.create("box_z_far", self.zfar)
        this_source_h5_group.attrs.create("box_z_near", self.znear)


    def uniform_position_in_source(self):
        return self.uniform_position_in_cube()


    def uniform_position_in_cube(self):
        p = Vec3d(self.random_in_range(self.xleft, self.xright),
                  self.random_in_range(self.ybottom, self.ytop),
                  self.random_in_range(self.znear, self.zfar))
        return p


    @classmethod
    def is_box_source(cls, conf_sec_name):
        return 'ParticleSourceBox' in conf_sec_name