Read Chapter 11 of the notes and work the following problems
(Don't forget How to do parallel performance analysis.)
Test your code by checking to see if the resulting matrix
U is all ones on and above the diagonal (compute and print
the sum of squares of differences between U_{ij} and 1 for
the appropriate i and j. If you compute L, do the same for
L (but the lower triangle are all ones: L=U' .)
For extra credit, run this with n=2048 on 2, 4, and 8.