Symmetric positive definite~(SPD) matrices have shown important value and applications in statistics and machine learning, such as FMRI analysis and traffic prediction. Previous works on spd matrices mostly focus on discriminative models, where predictions are made directly on $E(X|y)$