C++更改實例變量

Question

我正在編寫一個程序來計算n體問題。我遇到了一個問題，試圖改變我的代碼，所以它很容易適應任何數量的機構。函數corr（）似乎存在一個問題，在函數corr（）被使用之後，對該函數中的類粒子實例的某些變量所做的更改似乎會丟失。我如何解決這個問題？

#include <cstdlib> 
#include <iostream> 
#include <cmath> 
#include <fstream> 
#define h 1000.0 
#define G 6.67384*pow(10.0,-11) 

using namespace std; 


class particle{ 
     public: 
     double kx1,kx2,kx3,kx4, kv1, kv2, kv3, kv4; 
     double ky1, ky2, ky3, ky4, kvy1, kvy2, kvy3, kvy4; 
     double x,y,vx,vy,m; 


     double dist(particle aap){ 
      double dx = x - aap.x; 
      double dy = y - aap.y; 
      return sqrt(pow(dx,2.0)+pow(dy,2.0)); 
      } 

     double g(double x1, double y1,particle aap){ 
      return G*aap.m*(aap.x-x1)/pow(dist(aap),3.0); 
      } 

     double p(double x1, double y1, particle aap){ 
      return G*aap.m*(aap.y-y1)/pow(dist(aap),3.0); 
     } 


     void update(){   //zet het object 1 stap vooruit 
      x = x + (1/6.0)*(kx1+2*kx2+2*kx3+kx4); 
      vx = vx + (1/6.0)*(kv1+2*kv2+2*kv3+kv4); 
      y = y + (1/6.0)*(ky1+2*ky2+2*ky3+ky4); 
      vy = vy + (1/6.0)*(kvy1+2*kvy2+2*kvy3+kvy4); 
      } 

    void create(double x1, double y1, double vx1, double vy1, double m1){ 
         x = x1; 
         y = y1; 
         vx = vx1; 
         vy = vy1; 
         m =m1; 
         } 

    bool operator ==(particle &other){ 
      if(x == other.x && y == other.y && vx == other.vx && vy == other.vy){ 
       return true; 
       } 
       } 


     }; 



particle zon, maan, aarde; 

void set(){ 
    zon.create(1, 1, -2, 1, 2*pow(10.0,30)); 
    aarde.create(1.5*pow(10.0,11), 0, 2, 29780, 6*pow(10.0,24)); 
    maan.create(aarde.x + 1, aarde .y + 3.844399*pow(10.0,8), aarde.vx + -1022.0, aarde.vy + 1, 7.3347*pow(10.0,22)); 
    } 



double xforce(double x1, double y1, particle aap){  //kracht in de x-richting 

     particle bodies[] = {zon, aarde, maan}; 

     double fx; 
     for (int i = 0; i < 3; i++){ 
      if (bodies[i].x == aap.x && bodies[i].y == aap.y && bodies[i].vx == aap.vx && bodies[i].vy == aap.vy){;} 

      else{ 
       fx += aap.g(x1,y1,bodies[i]); 
       } 
       } 
     return fx; 
     } 

double yforce(double x1, double y1, particle aap){ //kracht in de y-richting 

     particle bodies[] = {zon, aarde, maan}; 

     double fy; 
     for (int i = 0; i <= 3; i++){ 
      if (bodies[i].x == aap.x && bodies[i].y == aap.y && bodies[i].vx == aap.vx && bodies[i].vy == aap.vy) {;} 

      else{ 
       fy += aap.p(x1,y1,bodies[i]); 
       } 
       } 
     return fy; 
     } 


void corr(particle body){ 
    body.kx1 = h*body.vx; 
    body.kv1 = h*xforce(body.x, body.y, body); 
    body.ky1 = h*body.vy; 
    body.kvy1 =h*yforce(body.x, body.y, body); 

    body.kx2 = h*(body.vx + 0.5*body.kv1); 
    body.kv2 = h*xforce(body.x + 0.5*body.kx1, body.y + 0.5*body.ky1, body); 
    body.ky2 = h*(body.vy + 0.5*body.kvy1); 
    body.kvy2 = h*yforce(body.x + 0.5*body.kx1, body.y + 0.5*body.ky1, body); 

    body.kx3 = h*(body.vx+ 0.5*body.kv2); 
    body.kv3 = h*xforce(body.x + 0.5*body.kx2, body.y + 0.5*body.ky2, body); 
    body.ky3 = h*(body.vy+ 0.5*body.kvy2); 
    body.kvy3 = h*yforce(body.x + 0.5*body.kx2, body.y + 0.5*body.ky2,body); 

    body.kx4 = h*(body.vx+body.kv3); 
    body.kv4 = h*xforce(body.x+ body.kx3, body.y + body.ky3, body); 
    body.ky4 = h*(body.vy + body.kvy3); 
    body.kvy4 = h*yforce(body.x + body.kx3, body.y + body.ky3, body); 
    } 







void bereken(){ 
    set(); 
    ofstream file; 
    file.open("3body.txt"); 
    for(int i =0; i <=30000; i++){ 
      corr(maan); 
      corr(zon); 
      corr(aarde); 
      zon.update(); 
      aarde.update(); 
      maan.update(); 
      file << i*h <<" "<< zon.x << " "<< zon.y << " "<< zon.vx<< " "<< zon.vy <<" "<< aarde.x << " " << aarde.y <<" "<< aarde.vx <<" " << aarde.vy <<" "<< maan.x<<" "<<maan.y<<"\n"; 
      } 
    file.close(); 
    } 


int main() 
{ 
    set(); 
    bereken(); 
    system("pause"); 
    return 0; 

} 

Answer 1

只需使用引用：

void corr(particle& body)

相同的修訂必須在其他地方被應用。你現在有一個代碼，COPIES你的對象，做一些計算就可以了，然後就刪除該臨時對象......通過使用引用（&）你有效地「共享」的主被叫之間的對象。

在你的類的方法你也應該利用這段時間提高效率，傳遞這個對象的引用，但 - 這是簡單的通過引用（通常爲指針一樣），但是複製這個巨大的物體需要時間。

Answer 2

弗雷迪肖邦提出的建議的工作和固定功能更正件（）。但是現在我發現xforce（）和yforce（）函數無法正常工作。他們應該返回一個特定粒子上的總力量，這是他們在第一步中完成的，但在此之後，他們似乎會返回錯誤的值。任何人都有一個想法，爲什麼發生？