Unexpected behavior of CRL + crystal about shadow3 HOT 2 OPEN

Thanks Sergey for your post. I'm looking to it and it is hard to interpret... or to identify what's wrong.
To be continued...

from shadow3.

There is indeed a bug in shadow3.
It happens that the lens switch polarization of about half of the intensity.
This makes not sense.
I will look at that in detail and fix shadow3, but this will take some time.

For the time being, you can patch the beam after the lens to put all the intensity as s-polarized using a script:

where the script is:

import numpy

out_object = in_object_1.duplicate()
intensity = out_object._beam.getshonecol(23,nolost=0)
print(intensity)

out_object._beam.rays[:,6] = numpy.sqrt(intensity)
out_object._beam.rays[:,7] = 0.0
out_object._beam.rays[:,8] = 0.0

out_object._beam.rays[:,15] = 0.0
out_object._beam.rays[:,16] = 0.0
out_object._beam.rays[:,17] = 0.0

this changes the output from before:

to now:

Thanks a lot Sergey for reporting this problem!

from shadow3.