pymatj / pygtm Goto Github PK

Hi again,

a question regarding the sorting of the qs:

this statement is used to detect birefringes for a layer and to decide if eq (13) or (15) is used for qualculation of the Cps.

I am wondering if the birefringes criterium for a single layer can't be done by checking for off-diagonal elements in the Euler-transformed epsilon 3x3 matrix?

in the best case one saves to calculate Cp_t1 and Cp_t2 and it would be also more simple and intuitive in my opinion.

Best

Daniel

Hi @pyMatJ ,

I am currently checking the GTM core code I stumbled a bit over this code here:

Shouldn't the if-else-case be within the for loop .
As I understood it, you need to sort according to the four cases:

1. complex number with imaginary part >= 0: transmode[0]
2. real number with >= 0: transmode[1]
3. complex number with imaginary part < 0: reflmode[0]
4. real number with < 0: reflmode[1]

Maybe you can check and tell me the idea of the current code.

Thanks

Daniel