-1

私はC ++にかなり慣れていないので、ポインターの使用方法を学ぼうとしています。座標を作成し、乱数ジェネレーターを使用してランダムな方向に移動する次のファイルがあります。

値 sigmaf_point は、テキスト ファイルから入力されます。

 void methane_coords(double *&sigmaf_point)

 double dummy_int = 1;
 string dummystring;
        string s;

        ifstream Dfile;
        std::stringstream out;

        out << 1;
        s = out.str() + ".TXT";
        Dfile.open (s.c_str());

        if (Dfile.fail())
        {
            return;
        }

        for (int i=0; i<dummy_int; i++)
        {
        Dfile >> sigmaf_point[i];
        }

次に、別の関数で使用します。

double initial_energy(double **coords_fluid, const double *box_size){


// Loop over all pairs of atoms and calculate the LJ energy
double total_energy = 0;

for (int i = 0; i <= n_atoms-1; i++)
{   

        sf1=sigmaf_point(coords_fluid[i][3]);
        ef1=epsilonf_point(coords_fluid[i][3]);

        // Energy fluid-fluid
        for (int j = i+1; j <= n_atoms-1; j++)
        {               
            sf2=sigmaf_point(coords_fluid[j][3]);
            ef2=epsilonf_point(coords_fluid[j][3]);             

            double delta_x = coords_fluid[j][0] - coords_fluid[i][0];
            double delta_y = coords_fluid[j][1] - coords_fluid[i][1];
            double delta_z = coords_fluid[j][2] - coords_fluid[i][2];

            // Apply periodic boundaries
            delta_x = make_periodic(delta_x, box_size[0]);
            delta_y = make_periodic(delta_y, box_size[1]);
            delta_z = make_periodic(delta_z, box_size[2]);

            // Calculate the LJ potential
            s=(sf1+sf2)/2.0;
            e=pow((ef1*ef2),0.5);
            double r = pow((delta_x*delta_x) + (delta_y*delta_y) +
                      (delta_z*delta_z),0.5)/s;

        double e_lj = 4*((1/pow(r,12.0))-(1/pow(r,6.0))/e);

        total_energy = (total_energy + e_lj);
        }
    }

coords_fluid は次のようにメイン ファイルに作成されます。

double **coords_fluid = new double*[5000];

問題は sf1=sigmaf_point(coords_fluid[i][3]); にあります。

sigmaf_point で「式には関数型へのポインタが必要です」というエラーが表示されます。私はこれについて少し混乱しています。変数を呼び出す方法については知っていますが、修正できないようです。

乾杯

4

1 に答える 1

0

まず第一に、ポインターはすでに一種の参照であるため、ポインターへの参照はまったく役に立ちません。したがって、またはに変更しdouble *&ます。それはより速くなります。double *double &

sigmaf_pointまた、関数および配列として使用していることがわかります。どちらですか？の宣言をしていただけますsigmaf_pointか？

それが配列の変更であると仮定すると

sf1 = sigmaf_point(coords_fluid[i][3]);

に

sf1 = sigmaf_point[coords_fluid[i][3]];
于 2012-12-12T12:04:01.500 に答える