def getSolution(self, inputWave): A = self.get_fullAi() B = self.get_fullBi() * spheroidal.get_Bin(inputWave, self.particle, self.nmax) x = -scipy.linalg.solve(A, B) return (x[0:self.nmax + 1], x[self.nmax + 1:])
def getMatrixSolution(self, inputWave): b_in = spheroidal.get_Bin(inputWave, self.particle, self.nmax) b_int = scipy.linalg.solve(self.get_BSm(), -b_in) b_sca = asarray(self.get_BRm() * b_int) return (b_sca, b_int)
def getSolution(self,inputWave): A = self.get_fullA() B = self.get_fullB() * spheroidal.get_Bin(inputWave, self.particle, self.nmax) x = -scipy.linalg.solve(A, B) return (x[0:self.nmax + 1], x[self.nmax + 1:])