all matrices (K,F,C) inside the code are now sparse

better than passing around structured arrays … all the matrices are
first constructed locally, and in the routine AssembleArlequin they are
passed into global matrices with the full size. the monte carlo
matrices are also directly mounted at the full size so they can be
directly summed to K in the SolveArelquin routine