一维谐波振荡器递归函数，可能Infinit循环，C++语法

Question

使用递归函数工作在一维谐波振荡器上，代码在fortran中工作并尝试转换为C++以学习新的语法。

我修复了迄今为止通过google搜索到的所有通知错误，但仍然没有得到结果。该计划看似永远运行，没有公布结果。

我认为可能的答案可能是它没有被计算出来或有一些infinit循环，因为我的循环或调用函数中的语法问题？但我没有看到它，也没有发现它们的错误。

有关为何该程序无法正常工作的任何建议？

// 
 
// main.cpp 
 
// 1-d HO 
 
// 
 
// Created by Grant Metheny on 3/2/16. 
 
// Copyright (c) 2016 Grant Metheny C++ Codes. All rights reserved. 
 
// 
 

 
#include <iostream> 
 
#include <vector> 
 
#include <string> 
 
#include <fstream> 
 
#include <cmath> 
 
#include <math.h> 
 

 
using namespace std; 
 

 
int i = 0; 
 
int n = 0; 
 
double x = 1.; 
 
double xmax = 5; 
 
double imax = 1000; 
 
double wavefunc = 0; 
 
double fact = 1; 
 
double hpol = 1; 
 
double wf0 = 1; 
 
double wf1 = 1; 
 
double wf2 = 1; 
 
double wf3 = 1; 
 
double wf4 = 1; 
 
double wf5 = 1; 
 
double wf6 = 1; 
 

 
double wavefunction(int n, double x) 
 
{ 
 
    return wavefunc = pow(2.0,-(n*.5)) * pow(M_PI,.25) * exp(-(.5*pow(x,2.0))); 
 
} 
 

 
double factorial(int n) 
 
{ 
 
    for (i = 0; i <= n; i++) 
 
     
 
     if (i == 0) 
 
      fact = 1.; 
 
     else 
 
      fact = fact * i; 
 
    
 
    return fact; 
 
} 
 

 
double hermite(int n, double x) 
 
{ 
 
    for (i = 0; i <= n; i++) 
 
     if (i==1) 
 
      hpol = 1.0; 
 
     else if (n==1) 
 
      hpol = 2*x; 
 
     else 
 
      hpol = 2*x*hermite(n-1,x) - 2*(n-1)*hermite(n-2,x); 
 
    
 
    return hpol; 
 
} 
 

 
double dx = 2*xmax/imax; 
 

 

 
int main(int argc, const char * argv[]) { 
 
    
 
    
 
    for (i=0; i <= imax; i++) { 
 
    
 
    x = 5. - dx*i; 
 
    
 
    n = 0; 
 
    wf0 = hermite(n,x) * wavefunction(n,x) * pow(factorial(n),(-(.5))); 
 
    
 
    n = 1; 
 
    wf1 = hermite(n,x) * wavefunction(n,x) * pow(factorial(n),(-(.5))); 
 
    
 
    n = 2; 
 
    wf2 = hermite(n,x) * wavefunction(n,x) * pow(factorial(n),(-(.5))); 
 
    
 
    n = 3; 
 
    wf3 = hermite(n,x) * wavefunction(n,x) * pow(factorial(n),(-(.5))); 
 
    
 
    n = 4; 
 
    wf4 = hermite(n,x) * wavefunction(n,x) * pow(factorial(n),(-(.5))); 
 
    
 
    n = 5; 
 
    wf5 = hermite(n,x) * wavefunction(n,x) * pow(factorial(n),(-(.5))); 
 
    
 
    n = 6; 
 
    wf6 = hermite(n,x) * wavefunction(n,x) * pow(factorial(n),(-(.5))); 
 
    
 
    cout <<"I="<< i <<"X="<< x <<"WF0="<< wf0<<"WF1=" << wf1; // wf2, wf3, wf4, wf5, wf6 
 
    
 
} 
 
    return 0; 
 
}

Answer 1

你for循环不终止。取而代之的是：

for (i = 0; n; i++) 

（这显然是用Fortran启发），你需要这样做：

for (i = 0; i < n; i++) 

的想法是，for循环的第二部分不是限制;它是一个谓词，应该评估为true（以保持循环再进行一次迭代）或false（退出循环）。因为无论何时预期为真/假值，C/C++都会将非零整数值视为true，因此循环会继续进行。

您需要在代码中的几个位置进行修正。