fix issue with msort() using quick sort algorithm

NOTE: this is not a functional commit.
The trend of refractive indices vs thickness of overlay of LPFG seems
wrong. Besides, there are some refractive indices with real part equals
to zero and image part is positive. Do not know what is the reason.
add TODO item
refractor organization of libs, with references of CAMFR