diff --git a/src/USER-EWALDN/ewald_n.cpp b/src/USER-EWALDN/ewald_n.cpp index c0668944b3..affca321c3 100644 --- a/src/USER-EWALDN/ewald_n.cpp +++ b/src/USER-EWALDN/ewald_n.cpp @@ -93,13 +93,14 @@ void EwaldN::init() //dielectric = force->dielectric; mumurd2e = dielectric = 1.0; + int tmp; Pair *pair = force->pair; - int *ptr = pair ? (int *) pair->extract("ewald_order") : NULL; - double *cutoff = pair ? (double *) pair->extract("cut_coul") : NULL; + int *ptr = pair ? (int *) pair->extract("ewald_order",tmp) : NULL; + double *cutoff = pair ? (double *) pair->extract("cut_coul",tmp) : NULL; if (!(ptr||cutoff)) error->all("KSpace style is incompatible with Pair style"); int ewald_order = ptr ? *((int *) ptr) : 1<<1; - int ewald_mix = ptr ? *((int *) pair->extract("ewald_mix")) : GEOMETRIC; + int ewald_mix = ptr ? *((int *) pair->extract("ewald_mix",tmp)) : GEOMETRIC; memset(function, 0, EWALD_NFUNCS*sizeof(int)); for (int i=0; i<=EWALD_NORDER; ++i) // transcribe order if (ewald_order&(1<ntypes; if (function[1]) { // geometric 1/r^6 - double **b = (double **) force->pair->extract("B"); + double **b = (double **) force->pair->extract("B",tmp); delete [] B; B = new double[n+1]; bytes += (n+1)*sizeof(double); for (int i=0; i<=n; ++i) B[i] = sqrt(fabs(b[i][i])); } if (function[2]) { // arithmetic 1/r^6 - double **epsilon = (double **) force->pair->extract("epsilon"); - double **sigma = (double **) force->pair->extract("sigma"); + double **epsilon = (double **) force->pair->extract("epsilon",tmp); + double **sigma = (double **) force->pair->extract("sigma",tmp); if (!(epsilon&&sigma)) error->all("epsilon or sigma reference not set by pair style in ewald/n"); double eps_i, sigma_i, sigma_n, *bi = B = new double[7*n+7];