def __init__(self, N, d): PointSet.__init__(self, N, d) num_triangles = self.N * (self.N - 1) * (self.N - 2) / 6 self.num_angles = int(3 * num_triangles) self.theta = np.empty([ self.num_angles, ]) self.theta_tensor = np.empty([N, N, N]) self.corners = np.empty([self.num_angles, 3])
def __init__(self, N, d): from scipy import special PointSet.__init__(self, N, d) self.T = self.N * (self.N - 1) * (self.N - 2) / 6 self.M = int(3 * self.T) self.theta = np.empty([ self.M, ]) self.theta_tensor = np.empty([N, N, N]) self.corners = np.empty([self.M, 3]) self.abs_angles = np.empty([self.N, self.N]) self.C = 0 self.A = np.empty((self.C, self.M)) self.b = np.empty((self.C, 1))