Hi dear Greg,
Thank you for your response.
I initially ran this test with Python and then on Linux and got same results. The configuration of the Python code is similar to what Taoning provided in his workshop.
I have also uploaded all the inputs(ssky, rad files, grid points, RGB mat and spectrum mat) and simulated outputs (-pM and -pY values for spectrum mat and RGB mat) here .
my_scene = pr.Scene(case_name, surfaces=[floor, wall, ceiling, window], materials=[mat], sources=[sky, spectra_sky])
params = pr.SamplingParameters(
I=True,
ab=3,
ad=4096,
aa=0,
lw=0.01,
co=True,
cs=20,
)
grid_rays = "\n".join(" ".join(map(str, row)) for row in grid)
result_p = pr.rtrace(grid_rays.encode(), my_scene.octree, params=params.args()+['-pY'], header=False, outform='f')
result_m = pr.rtrace(grid_rays.encode(), my_scene.octree, params=params.args()+['-pM'], header=False, outform='f')