I was working around with the decomposition of the multidimensional linear canonical transform (which is not even continuous w.r.t. the parameters) into a few fractional Fourier transforms (and other related operators that are real analytic w.r.t. the parameters); this requires decomposition of the LCT matrix, and long story short, it is necessary to be able to decompose an arbitrary real matrix $B$ into a product $D¹P¹SP²D²$, with $D¹$ and $D²$ being diagonal, $P¹$ and $P²$ being positive definite, $S$ being symmetric, and all being real. 

So far I could not find any combination of known mutliplicative matrix decomposition methods that arrives at this result; I need help.

Note: Here $D¹$ and $D²$ are completely indepedent, the $¹$ and $²$ having absolutely no meanings, and so on.