/// <summary>
 /// Create a new square sparse identity matrix where each diagonal value is set to One.
 /// </summary>
 public static DiagonalMatrix CreateIdentity(int order)
 {
     return(new DiagonalMatrix(DiagonalMatrixStorage <Complex> .OfValue(order, order, One)));
 }