SubSparseMatrix class
make CCD
well, there's something need to be changed: CDcomp is not Hermite in BCS calculations
Seems everything involving Operators.C is working