1)

2)

for i = 1, number_of_tree_nodes do

3)

Z_i = node of tree;

4)

for all Z_j child of Z_i do

5)

prepare update of Z_i by j-th child;

6)

;

7)

end

8)

Calculate LL^T decomposition of Z_i;

9)

end