C++函数每次调用返回的vector值始终相同问题排查

我尝试在main函数中调用函数更新vector的值，再通过辅助函数从中获取两个特定值。但运行代码时，gibbsInfo（目标vector）的值在main的for循环每次迭代中都完全相同。我认为自己对main中声明的vector<vector<double>>的作用域及其与gibbsInformation函数返回的vector<vector<double>>的关系存在根本性误解。以下是相关代码，恳请提供排查建议：

#include<iostream>
#include<cmath>
#include<vector>

using namespace std;

vector<vector<double>> gibbsInformation(double, double);
vector<double> commonTangent(vector<vector<double>>, double);

int main(){

    for(int j = 1000; j <= 1500; j += 25){
        vector<vector<double>> gibbsInfo = gibbsInformation(0.5, j);
        vector<double> coords = commonTangent(gibbsInfo, j);

        cout << "(" << coords[0] << ", " << coords[2] << ")," << endl;
        cout << "(" << coords[1] << ", " << coords[2] << ")," << endl;

        for(int i = 0; i < gibbsInfo[0].size(); i++){
            cout << "Xa Value: " << gibbsInfo[0][i] << " Gibbs Value: " << gibbsInfo[1][i] << endl;
        }
    }

    return 0;
}

vector<vector<double>> gibbsInformation(double omega, double temp){

    vector<vector<double>> gibbs(2, vector<double>(200,0.0));

    int index = 0;

    for(double x = -0.5; x <= 1.5; x += 0.01){
        gibbs[0][index] = x;
        gibbs[1][index] = x + (0.9 * (1 - x)) + (omega * x * (1 - x)) + (((1000 - (int) temp) / 1000) * x * (1 - x)) + pow(0.5 - x,4);
        index++;
    }

    return gibbs;
}

vector<double> commonTangent(vector<vector<double>> gibbs, double temp){

    vector<double> coordinates(3, 0.0);

    coordinates[2] = temp;//set the third entry of coordinates to the temperature
    double minDistance = 1000;

    for(int i = 1; i < gibbs[1].size(); i++){

        //find the equation for the line
        double slope = (gibbs[1][i] - gibbs[1][i - 1]) / 0.01;

        if(i + 25 < gibbs[1].size()){//translates down by a distance of 0.25 away from the current point in the x direction
            for(int l = i + 25; l < gibbs[1].size(); l++){
                if(abs(gibbs[1][l] - (gibbs[1][i] + (slope * 0.01 * (l - i)))) < minDistance){
                    minDistance = abs(gibbs[1][l] - (gibbs[1][i] + (slope * 0.01 * (l - i))));
                    coordinates[0] = gibbs[0][i];
                    coordinates[1] = gibbs[0][l];
                }
            }
        }
    }

    return coordinates;
}

排查建议 #

• 核心问题：整数除法导致温度项失效
  看gibbsInformation函数里的这行计算：

  (((1000 - (int) temp) / 1000) * x * (1 - x))
  

  这里(int)temp是1000到1500的整数，1000 - (int)temp的结果是0到-500的整数，和1000做整数除法时，结果永远是0（C++中整数除法会截断小数部分）。这就导致温度变化完全不会影响吉布斯自由能的计算结果，所以每次返回的gibbsInfo自然完全相同。

• 作用域误解澄清
  你在main循环里每次声明的gibbsInfo都是全新的局部变量，每次都会接收gibbsInformation返回的新vector，作用域逻辑是完全正确的——问题根本不在作用域，而是计算逻辑里的整数除法陷阱。

• 修复与优化点

  1. 修复温度项计算：把整数除法改成浮点除法，比如把1000改成1000.0，同时去掉没必要的(int)强转：

     (((1000.0 - temp) / 1000.0) * x * (1 - x))
     

     这样温度变化就能正确反映到吉布斯自由能的计算中。
  2. 重置最小值变量：commonTangent函数里的minDistance只初始化了一次，后续调用时会沿用之前的最小值，导致结果错误。建议把double minDistance = 1000;移到外层for循环内部，或者每次调用函数时重新初始化。

内容的提问来源于stack exchange，提问作者Nathaniel Hess