added std.err estimate in sepdpsolve
work in progess: smepdpsolve