Do you have any suggestions how I can speed up the “f_img_update” function in the code below? I tried adding numba @jit decorator, but without luck. I guess that numpy vectorization might work, but I am not enough educated to make it work on my own :/

I have two numpy arrays – one containing coordinates of 3D points and the second with triangles composed from those points. I need to calculate top surface 2D area and volume of this triangles (tetrahedrons?) with base height of the lowest point.

I have a list of lists of indices like so:
outer_list = [[0], [1, 2], [1, 2], [3, 5], [4], [3, 5]]
Each inner list contains the index number of itself within the outer list.

Aim: Write a Python function that takes as input a matrix X of shape (M, p) where each row is a different vector x and compute the expression above for each of these rows, thus obtaining an output of shape (M, p).