请大家帮满看看这个程序为什么停止运行

jszjvictor

1. //general scalar equation: d(rho*phi)/dt=d/dx(rho*D*(dphi/dx))
   
2. //W=West(i); C=centre(i+1); E=East(i+2)
   
3. //w=west(i+1/2); e=east(i+3/2)
   
4. 
   
5. #include<iostream>                                   
   
6. #include<fstream>                                      //for files
   
7. using namespace std;
   
8. 
   
9. void Square(int n, double dx, double a, double b, double x[], double y[])
   
10. {
    
11.      int i;
    
12.      for(i=0; i<n; i++)
    
13.     {
    
14.              x[i]=(i+0.5)*dx;                          //define initial x value for each cell, mid-point
    
15.                  if (i>a&&i<b)                             //define initial y value for each cell   
    
16.                  {
    
17.                      y[i]=1.0;
    
18.                  }                                
    
19.                  else
    
20.                  {
    
21.                      y[i]=0.0;     
    
22.                  }                        
    
23.     }
    
24. }
    
25. 
    
26. double Diffusion(int i, double dx, double y[], double rho[], double D[])
    
27. {
    
28.      double yW, yC, yE, rhoW, rhoC, rhoE, dW, dC, dE, phiW, phiC, phiE, phiw, phie, rhodw, rhode;
    
29.      yW=y[i]; yC=y[i+1]; yE=y[i+2];                               //y at W, C and E
    
30.      rhoW=rho[i]; rhoC=rho[i+1]; rhoE=rho[i+2];                   //density at W, C and E
    
31.      dW=D[i]; dC=D[i+1]; dE=D[i+2];                               //diffusivity at W, C and E
    
32.      phiW=yW/rhoW; phiC=yC/rhoC; phiE=yE/rhoE;                    //y divided by density to get phi
    
33.      phiw=phiC-phiW; phie=phiE-phiC;                              //dleta phi
    
34.      rhodw=(rhoW*dW+rhoC*dC)/2; rhode=(rhoE*dE+rhoC*dC)/2;        //density*diffusivity at w and e
    
35.      return((rhode*phie-rhodw*phiw)/(dx*dx));                     //diffusion term
    
36. }
    
37. 
    
38. void Output(int n, double x[], double y[])
    
39. {
    
40.     int i;
    
41.     ofstream outdata;                                            //libray for output                     
    
42.     outdata.open("output.dat");                                  //open the file
    
43.     if(outdata.is_open())                                        //check if the file is open
    
44.     {
    
45.              for(i=0; i<n; i++)                                  //for each point
    
46.              {
    
47.                       outdata<<x[i]<<" "<<y[i]<<endl;            //write the value of x and y to file
    
48.              }
    
49.     outdata.close();                                             //close the file
    
50.     }
    
51.     else                                                         //if _open() return falsch
    
52.     {
    
53.              for(i=0; i<n; i++)
    
54.              {
    
55.                       cout<<x[i]<<" "<<y[i]<<endl;
    
56.              }
    
57.     }
    
58. }
    
59. 
    
60. int main()
    
61. {
    
62.     double l=1.0;                                      //total length
    
63.     double tpmax=1.0;                                  //total time
    
64.     double dt=0.01;                                    //time step interval
    
65.     const int n=100;                                   //define cell number
    
66.     double dx, a, b, tp;
    
67.     double x[n], y[n], rho[n], D[n], RHS[n], ybuffer[n];
    
68.     dx=l/n;                                            //calculate the length of each cell
    
69.     a=0.3*n;                                           //define the critical point
    
70.     b=0.7*n;                                           //define the critical point
    
71.     int i;
    
72.     for(i=0; i<n; i++)                                 //define the density and diffusivity values of air
    
73.     {
    
74.              rho[i]=1.225;
    
75.              D[i]=0.00002;
    
76.     }
    
77.                                                         
    
78.     Square(n, dx, a, b, &x[n], &y[n]);                               //switch on square
    
79.    
    
80.     for(tp=0.0; tp<tpmax; tp=tp+dt)
    
81.     {
    
82.               for(i=0; i<(n-1); i++)
    
83.               {
    
84.                        RHS[i+1]=Diffusion(i, dx, &y[n], &rho[n], &D[n]);                //switch on diffusion                                 
    
85.                        ybuffer[i+1]=y[i+1]+dt*RHS[i+1];                                 //explicit time integration to get y
    
86.                        y[i+1]=ybuffer[i+1];
    
87.               }
    
88.     }
    
89.    
    
90.     Output(n, &x[n], &y[n]);                                         //results output
    
91.    
    
92.     return 0;                                                        // finish
    
93. }

复制代码

没有报错，但是在.exe中会停止运行

bruceteen

Square(n, dx, a, b, &x[n], &y[n]);
--->
Square(n, dx, a, b, x, y);

littledick

bruceteen 发表于 2013-01-04 09:50
Square(n, dx, a, b, &x[n], &y[n]);
--->
Square(n, dx, a, b, x, y);

我觉得LZ问的可能不是这个。

jszjvictor

回复 2# bruceteen


    谢谢，现在能运行完了

jszjvictor

回复 3# littledick

谢谢，能运行了，但可能还有其他问题，但最起码不会停止了
   

jszjvictor

回复 3# littledick

谢谢，能运行了，但可能还有其他问题，但最起码不会停止了
   

jszjvictor

回复 3# littledick

谢谢，能运行了，但可能还有其他问题，但最起码不会停止了
   

jszjvictor

回复 3# littledick

谢谢，能运行了，但可能还有其他问题，但最起码不会停止了