Python stops working with max_order over 31

I am using pyroomacoustics purely to generate room impulse responses. I would like to generate more reverberant rooms, of at least one second.

The highest max_order that does not trigger “Python stopped working” is 31. Same behavior on two different computers. I could not find any documentation concerning this. Is it an issue, or am I doing something wrong? I attached a small code snippet that reproduces the error on my machines.

import pyroomacoustics as pra
import numpy as np

def compute_rir(order):
    fromPos = np.zeros((3))
    toPos = np.ones((3,1))
    roomSize = np.array([3,3,3])
    room = pra.ShoeBox(roomSize, fs = 1000, absorption=0.95, max_order=order)
    room.add_source(fromPos)
    mics = pra.MicrophoneArray(toPos, room.fs)
    room.add_microphone_array(mics)
    room.compute_rir()

compute_rir(31)
print("Completed 31")
compute_rir(32)
print("Completed 32")

[Pvjesper]

Thank you! I only have my laptop with me, but I tried it for a few orders, up to 95, and it works well. So I think it is really solved.

[fakufaku]

@Pvjesper The latest version update includes a complete rewrite of the simulation core in C++ which seems to solve the problem. I have added a test for order over 31 that seems to run without problem. https://github.com/LCAV/pyroomacoustics/actions/runs/132155203 Can you try to update to v0.4.0 and let me know it this solves your problem ?