Fix codi

sdpm/src/sdpmSolver.cpp

@@ -547,12 +547,12 @@ void Solver::post(size_t imd, size_t ifq, sdpmVectorXcd const &etau, sdpmVectorX
         for (size_t i = 0; i < elems.size(); ++i)
         {
             // phi_t, phi_x, phi_x(0), u^2
-            sdpmComplex phit = -sdpmComplex(0, 1) * omega * emu(_rows[elems[i]]);
+            sdpmComplex phit = -sdpmComplex(0., 1.) * omega * emu(_rows[elems[i]]);
             sdpmComplex phix = u1[i](0) - (uc0[i](0) + uc1[i](0) * sdpmComplex(0., 1.) * omega);
-            sdpmDouble phix0 = _u[elems[i]->getId() - 1](0) - ui(0);
-            sdpmComplex cp1 = -sdpmComplex(2, 0) * phit + mi * mi * (2 * phix0 * phix + phix0 * phit);
+            sdpmComplex phix0(_u[elems[i]->getId() - 1](0) - ui(0), 0.);
+            sdpmComplex cp1 = -sdpmComplex(2., 0.) * phit + sdpmComplex(mi * mi, 0.) * (sdpmComplex(2., 0.) * phix0 * phix + phix0 * phit);
             for (size_t j = 0; j < 3; ++j)
-                cp1 -= 2 * _u[elems[i]->getId() - 1](j) * u1[i](j);
+                cp1 -= sdpmComplex(2 * _u[elems[i]->getId() - 1](j), 0.) * u1[i](j);
             // cp
             _cp1Re[imd][ifq][elems[i]->getId() - 1] = cp1.real();
             _cp1Im[imd][ifq][elems[i]->getId() - 1] = cp1.imag();
@@ -579,7 +579,7 @@ void Solver::post(size_t imd, size_t ifq, sdpmVectorXcd const &etau, sdpmVectorX
         b->setUnsteadyMomentCoeff(imd, ifq, bcm1);
 
         // compute mean load coefficient
-        sdpmVector3d cf0(0., 0., 0.);
+        sdpmVector3cd cf0 = sdpmVector3cd::Zero();
         std::vector<sdpmVector3d> loads = b->getLoads();
         for (size_t ii = 0; ii < loads.size(); ++ii)
             cf0 += loads[ii];