Adapt common blocks structure to cluster case

src/include/Commons.h

@@ -214,6 +214,8 @@ public:
 	//! \brief QUESTION: definition?
 	std::complex<double> *vj0;
 	//! \brief QUESTION: definition?
+	std::complex<double> *vj;
+	//! \brief QUESTION: definition?
 	std::complex<double> *vyhj;
 	//! \brief QUESTION: definition?
 	std::complex<double> *vyj0;

src/libnptm/Commons.cpp

@@ -99,6 +99,7 @@ C1_AddOns::C1_AddOns(C4 *c4) {
 	nlemt = 2 * c4->nlem;
 	vh = new complex<double>[(nsph * nsph - 1) * c4->litpo];
 	vj0 = new complex<double>[nsph * c4->lmtpo];
+	vj = new complex<double>[1]; // QUESTION: is 1 really enough for a general case?
 	vyhj = new complex<double>[nsph * (nsph - 1) * c4->litpo];
 	vyj0 = new complex<double>[nsph * c4->lmtpo];
 	am0m = new complex<double>*[nlemt];