Before revision:  We prove  rchol is breakdown free and apply it to solving  linear systems with symmetric diagonally-dominant matrices.
After revision:  We prove  rchol is breakdown free and apply it to solving large sparse linear systems with symmetric diagonally-dominant matrices.