import subprocess as sp
import numpy as np
import k3d

def getFacets3D(filename):
    exe = ["./getFacets3D", filename]
    p = sp.Popen(exe, stdout=sp.PIPE, stderr=sp.PIPE)
    stdout, stderr = p.communicate()
    temp1 = stderr.decode("utf-8")

    lines = temp1.strip().split('\n')

    all_points = []
    all_indices = []
    current_index = 0

    for line in lines:
        if line.strip() == '':
            # Move to next plane
            if len(all_points) > current_index:
                # Calculate indices for the previous plane
                n = len(all_points) - current_index
                indices = np.vstack((np.zeros(n-2, dtype=np.uint32) + current_index,
                                     np.arange(1, n-1, dtype=np.uint32) + current_index,
                                     np.arange(2, n, dtype=np.uint32) + current_index)).T.flatten()
                all_indices.append(indices)
                current_index = len(all_points)
        else:
            # Convert each line into a tuple of floats (x, y, z) and add it to all_points
            point = tuple(map(float, line.split()))
            all_points.append(point)

    # Adding indices for the last plane if it exists
    if len(all_points) > current_index:
        n = len(all_points) - current_index
        indices = np.vstack((np.zeros(n-2, dtype=np.uint32) + current_index,
                             np.arange(1, n-1, dtype=np.uint32) + current_index,
                             np.arange(2, n, dtype=np.uint32) + current_index)).T.flatten()
        all_indices.append(indices)

    # Convert lists to numpy arrays
    all_points = np.array(all_points, dtype=np.float32)
    all_indices = np.concatenate(all_indices).astype(np.uint32)

    return all_points, all_indices