Question

使用递归函数处理1-D谐波振荡器，代码在fortran中工作并尝试转换为C ++以学习新语法。

到目前为止，我已修复了通过谷歌搜索的所有通知错误，但仍未获得结果。该计划似乎永远不会发布结果。

我认为可能的答案是，由于我的for循环或被调用函数中的语法问题，可能根本没有计算或者有一些无限循环？但我没有看到它，也没有发现它们的错误。

有关此计划无法正常运作的建议吗？

＆＃13;

//
//  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;
}

＆＃13;

Answer 1

您的for循环未终止。而不是：

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

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

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

这个想法是for循环的第二部分不是限制;它是一个谓词，应该评估true（保持循环再循环一次）或false（退出循环）。因为只要期望真/假值，C / C ++就会将非零整数值视为true，所以循环只会继续。

您需要在代码中的多个位置进行更正。

一维谐振子递归函数，可能的无限循环，c ++语法

1 个答案: