public static void GetCATMatrix(Illuminant source, Illuminant dest, ChromaticAdaptionMethod method = ChromaticAdaptionMethod.BradFord)
{
GetMatrix(method, ref Ma, ref MaPrime);
double[] prb_source = Matrix.Dot(Ma, source.GetXYZ());
double[] prb_dest = Matrix.Dot(Ma, dest.GetXYZ());
double[,] prb = Matrix.Diagonal(3, Elementwise.Divide(prb_dest, prb_source));
M = Matrix.Dot(Matrix.Dot(MaPrime, prb), Ma);
}
