@@ -217,21 +217,29 @@ template <typename T> struct densvars {
217217 b_43 = pow (b, 4.0 / 3.0 );
218218 }
219219
220- const XCFunctional * parent;
220+ const XCFunctional * parent{ nullptr } ;
221221 double get_param (enum xc_parameter p) const { return parent->settings [p]; }
222222
223- T a, b, gaa, gab, gbb;
224- /* na+nb, na-nb, (grad n)^2, (grad n).(grad s), (grad s)^2 */
225- T n, s, gnn, gns, gss;
226-
227- T tau, taua, taub; // Kinetic energy densities.
228-
229- T lapa, lapb; // Density Laplacians
230-
231- T zeta; // s/n
232- T r_s; // (3/4pi)^1/3*n^(-1/3)
233- T n_m13; // pow(n,-1.0/3.0)
234- T a_43, b_43; // pow(a,4.0/3.0), pow(b,4.0/3.0)
235-
236- T jpaa, jpbb; // square of the alpha and beta paramagnetic current vectors.
223+ T a{static_cast <T>(0 )};
224+ T b{static_cast <T>(0 )};
225+ T gaa{static_cast <T>(0 )};
226+ T gab{static_cast <T>(0 )};
227+ T gbb{static_cast <T>(0 )};
228+ T n{static_cast <T>(0 )}; // / na+nb
229+ T s{static_cast <T>(0 )}; // / na - nb
230+ T gnn{static_cast <T>(0 )}; // / (grad n) ^ 2
231+ T gns{static_cast <T>(0 )}; // / (grad n).(grad s)
232+ T gss{static_cast <T>(0 )}; // / (grad s) ^ 2
233+ T tau{static_cast <T>(0 )}; // / Kinetic energy density.
234+ T taua{static_cast <T>(0 )}; // / Alpha kinetic energy density.
235+ T taub{static_cast <T>(0 )}; // / Beta kinetic energy density.
236+ T lapa{static_cast <T>(0 )}; // / Alpha Laplacian density.
237+ T lapb{static_cast <T>(0 )}; // / Beta Laplacian density.
238+ T zeta{static_cast <T>(0 )}; // / s/n
239+ T r_s{static_cast <T>(0 )}; // / (3/4pi)^1/3*n^(-1/3)
240+ T n_m13{static_cast <T>(0 )}; // / pow(n,-1.0/3.0)
241+ T a_43{static_cast <T>(0 )};
242+ T b_43{static_cast <T>(0 )}; // / pow(a,4.0/3.0), pow(b,4.0/3.0)
243+ T jpaa{static_cast <T>(0 )}; // / square of the alpha paramagnetic current vector.
244+ T jpbb{static_cast <T>(0 )}; // / square of the beta paramagnetic current vector.
237245};
0 commit comments