不报错，也没警告，但想要的东西出不来

这是Gillespie算法，代码没显示错误和警告，但是跑出来只有时间序列，没有状态的数据
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Users\Lenovo\source\repos\GA_XYZW(example1)\Debug\GA_XYZW(example1).exe”。已加载符号。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\ntdll.dll”。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\kernel32.dll”。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\KernelBase.dll”。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\vcruntime140d.dll”。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\msvcp140d.dll”。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\ucrtbased.dll”。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\ucrtbased.dll”。
“GA_XYZW(example1).exe”(Win32): 已卸载“C:\Windows\SysWOW64\ucrtbased.dll”
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\kernel.appcore.dll”。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\msvcrt.dll”。
“GA_XYZW(example1).exe”(Win32): 已加载“C:\Windows\SysWOW64\rpcrt4.dll”。
线程 0x4380 已退出，返回值为 0 (0x0)。
线程 0x8274 已退出，返回值为 0 (0x0)。
线程 0x7ae8 已退出，返回值为 0 (0x0)。
程序“[38464] GA_XYZW(example1).exe”已退出，返回值为 0 (0x0)。




#include <iostream>
#include <cmath>
#include <fstream>
#include <string>
#include <stdlib.h>
using namespace std;

class System_States
{
    public:
    long Reaction, Reactions;
    double States[4] = {0};
    double A[6] = {0};
    double C[6] = { 0.1,0.03,0.1,0.01,0.2,0.001 };
    double A0=0 ;
    double Tau=0;

    System_States()
    {
        States[0] = 1.0;
        States[1] = 1.0;
        States[2] = 1.0;
        States[3] = 1.0;
        Reactions = 6;
    };

    void Propensity_function()
    {
        
        A[0] = C[0] * States[1];
        A[1] = C[1] * States[2];
        A[2] = C[2] * States[1] * (States[1] - 1.) / 2.;
        A[3] = C[3] * States[3];
        A[4] = C[4] * States[0] * States[1];
        A[5] = C[5] * States[1] * (States[1] - 1.) / 2.;
        A0 = A[0] + A[1] + A[2] + A[3] + A[4] + A[5];
    }

    void Update_Tau()
    {
        double R1;
        R1 = 1.0 * rand() / RAND_MAX;
        Tau = log(1.0 / R1) / A0;
    }

    void Update_Reaction()
    {
        double R2;
        double R2A0;
        R2 = 1.0 * rand() / RAND_MAX;
        R2A0 = R2 * A0;

        if (A[0] >= R2A0)
        {
            Reaction = 1;
        }
        else
        {
            double sum = A[0];
            for (int i = 1;i < Reactions; i++)
            {
                sum = sum + A[i];
                if (sum >=R2A0)
                {
                    Reaction = i + 1;
                    break;
                }
            }
        }

    }
    void Update_States()
    {
        switch (Reaction)
        {
            case 1:
            States[1] = States[1] - 1;
            States[2] = States[2] + 1;
            break;
            case 2:
            States[1] = States[1] + 1;
            States[2] = States[2] - 1;
            break;
            case 3:
            States[1] = States[1] - 2;
            States[3] = States[3] + 1;
            break;
            case 4:
            States[1] = States[1] + 2;
            States[3] = States[3] - 1;
            break;
            case 5:
            States[0] = States[0] - 1;
            States[1] = States[1] + 1;
            break;
            case 6:
            States[0] = States[0] + 1;
            States[1] = States[1] - 1;
            default:
                break;
        }

    }
};

void Gillespie(System_States States, long Maxsteps, long Maxtime, long SampleNumber)
{
    double T= 0.0;
    long step = 0;
    ofstream fout("data." + to_string(SampleNumber));
    while (T <= Maxtime && step <= Maxsteps)
    {
        fout << T << "" << States.States[0] << "" << States.States[1] << "" << States.States[2] << "" << States.States[3]  << "\n";
        States.Propensity_function();
        States.Update_Tau();
        States.Update_Reaction();
        States.Update_States();
        T= T+ States.Tau;
        step = step + 1;
    }
    fout.close();

};
int main()
{
    using namespace std;
    srand((unsigned int)time(NULL));
    System_States s1;
    Gillespie(s1, 1000, 100000, 1);
        return 0;

}


        
        





   
        
   

[此贴子已经被作者于2021-11-6 09:47编辑过]