コード例 #1
0
ファイル: linesink.py プロジェクト: mbakker7/timml
 def initialize(self):
     self.Ncp = self.order + 1
     self.z1 = self.x1 + 1j * self.y1
     self.z2 = self.x2 + 1j * self.y2
     self.L = np.abs(self.z1 - self.z2)
     self.thetaNormOut = np.arctan2(self.y2 - self.y1,
                                    self.x2 - self.x1) - np.pi / 2.0
     self.cosnorm = np.cos(self.thetaNormOut) * np.ones(self.Ncp)
     self.sinnorm = np.sin(self.thetaNormOut) * np.ones(self.Ncp)
     #
     self.xc, self.yc = controlpoints(self.Ncp, self.z1, self.z2, eps=0)
     if self.zcinout is not None:
         self.xcin, self.ycin = controlpoints(self.Ncp, self.zcinout[0],
                                              self.zcinout[1], eps=0)
         self.xcout, self.ycout = controlpoints(self.Ncp, self.zcinout[2],
                                                self.zcinout[3], eps=0)
     else:
         self.xcin, self.ycin = controlpoints(self.Ncp, self.z1, self.z2,
                                              eps=1e-6)
         self.xcout, self.ycout = controlpoints(self.Ncp, self.z1, self.z2,
                                                eps=-1e-6)
     if self.aq is None:
         self.aq = self.model.aq.find_aquifer_data(self.xc, self.yc)
     if self.addtomodel:
         self.aq.add_element(self)
     self.parameters = np.empty((self.Nparam, 1))
     self.parameters[:,
     0] = self.Qls / self.L  # Not sure if that needs to be here
コード例 #2
0
 def initialize(self):
     LineSinkHoBase.initialize(self)
     self.xcin, self.ycin = controlpoints(self.Ncp-1, self.z1, self.z2, eps=1e-6, include_ends=True)
     self.xcout, self.ycout = controlpoints(self.Ncp-1, self.z1, self.z2, eps=-1e-6, include_ends=True)
     if self.aqin is None:
         self.aqin = self.model.aq.find_aquifer_data(self.xcin[0], self.ycin[0])
     if self.aqout is None:
         self.aqout = self.model.aq.find_aquifer_data(self.xcout[0], self.ycout[0])
コード例 #3
0
 def initialize(self):
     self.Ncp = self.order + 1
     self.z1 = self.x1 + 1j * self.y1
     self.z2 = self.x2 + 1j * self.y2
     self.L = np.abs(self.z1 - self.z2)
     self.thetaNormOut = np.arctan2(self.y2 - self.y1,
                                    self.x2 - self.x1) - np.pi / 2.0
     self.cosnorm = np.cos(self.thetaNormOut) * np.ones(self.Ncp)
     self.sinnorm = np.sin(self.thetaNormOut) * np.ones(self.Ncp)
     #
     self.xc, self.yc = controlpoints(self.Ncp, self.z1, self.z2, eps=0)
     if self.zcinout is not None:
         self.xcin, self.ycin = controlpoints(self.Ncp,
                                              self.zcinout[0],
                                              self.zcinout[1],
                                              eps=0)
         self.xcout, self.ycout = controlpoints(self.Ncp,
                                                self.zcinout[2],
                                                self.zcinout[3],
                                                eps=0)
     else:
         self.xcin, self.ycin = controlpoints(self.Ncp,
                                              self.z1,
                                              self.z2,
                                              eps=1e-6)
         self.xcout, self.ycout = controlpoints(self.Ncp,
                                                self.z1,
                                                self.z2,
                                                eps=-1e-6)
     if self.aq is None:
         self.aq = self.model.aq.find_aquifer_data(self.xc, self.yc)
     if self.addtomodel:
         self.aq.add_element(self)
     self.parameters = np.empty((self.Nparam, 1))
     self.parameters[:,
                     0] = self.Qls / self.L  # Not sure if that needs to be here