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
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])