小编jan*_*nty的帖子

我为LBM编写了一个C++代码来解决扩散方程,我在Ubuntu中用g ++编译它.编译得很好,但是当我尝试执行程序时,它并没有结束.我认为该程序有某种无限循环,但我无法解决在哪里？提前谢谢你们.

这是我的代码:

#include <iostream>
#include <fstream>

using namespace std;

int main(){
    const int Lx = 100;
    const int Ly = 100;
    int sMod = 8; //speed model 
    float f[sMod][Lx+1][Ly+1];
    float rho[Lx+1][Ly+1];
    float feq, sum;
    float x[Lx+1];
    float y[Ly+1];
    float csq, alpha, omega,Tw;
    float w[9] = {4./9., 1./9., 1./9., 1./9., 1./9., 1./36., 1./36., 1./36., 1./36.};
    int i, j, dt, dx, k, dy; 

    FILE * mFile;
    dt = 1.0;
    dx = 1.0;
    dy = dx;

    x[0] = 0.0;
    for(i = 1; …