def value(self, acquisition_direction): accum = complex(0, 0) for l in range(self.order + 1): if l%2 == 0: for m in range(-l, l + 1): accum += self.__coefficients[self.__index(l,m), 0]*sph_harm(l, m, acquisition_direction) return accum
def __x(self, acquisition_directions, order): rows = [] for acquisition_direction in acquisition_directions: row = [] for l in range(order + 1): for m in range(-l, (l + 1)): row.append(sph_harm(l,m,acquisition_direction)) rows.append(row) return np.matrix(rows)