Wrong oriented bounding box

Hello

I am getting a wrong oriented bounding box for a TopoDS_Shell consisting of two rectangular faces in the XY plane rotated by 23 degrees about the Z-axis. The long side of the overall shape is 2 and short is 1 but I get 0.75 and 0.5.The body was imported using the STEP file importer and I used AddOBB function with theIsTriangulationUsed=False and theIsOptimal=True. The step file was created in an external program (Cubit) and eported.

If I take a single 2x1 rectangle, rotate it by the same amount, split it using OCC functions, it gives the right answer. If I export this internally split rectangle to a STEP file and import it back into OCC and compute the OBB, I get the right answer.

I used pythonOCC 7.7.0 but I have a feeling it is coming from the C++ OCC library

A picture of the model, a step file and a log file are attached