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:])
Example #2
0
 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:])
 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)