powered by
Default reference point
default_ref_pt(p)
A diagonal matrix of the desired dimension
the dimension