Question

以下代码描述了2个谐波振荡器。它们最初是非耦合和独立的，我只看其中一个，这是一个机械振荡器。另一个振荡器的变量已被声明并强制为0.我首先进行300,000次迭代，并为80个不同频率的wdm（机械驱动器w）执行此操作。

//find frequency response
int index_A;
double wdm_1;
wdm_1=wm-2*3.142*2e5;
double wdm_2;
wdm_2=wm+2*3.142*2e5;
double wdm_prec=2*3.142*5e3;
index_A=(wdm_2-wdm_1)/wdm_prec;
printf("%d \n", index_A);

However, my code runs till 58 frequencies and gives this error:

奇怪的是，如果我分别以40个频率的两个部分运行代码并将文件附加在一起，那么它可以正常工作。

Resonance peak in frequency domain

此外，当我将v0的大小减小到3时，代码正常运行。但是，我稍后会需要其他变量。

int j=0;
    for (j=0; j<= index_A ; j++){
        wdm=wdm_1+j*wdm_prec;
        printf("%d \n",j);
        v0[0] = 0;
        v0[1] = 0;
        v0[2] = 0;
        v0[3] = 0;
        v0[4] = 0;
        v0[5] = 0;
        v0[6]= wdm;

            for (i=0; i< n ; i++){
                if (cabs(xa)>=1){
                        printf("Breaking Loop \n");
                        break;
                        }

                v1 = rk4vec_ameya_complex_1 ( tau, 7, v0, dtau, rk4vec_f_ameya_complex_1 );
                memcpy(v0, v1, 4 * sizeof ( double complex ) );
                tau=tau+dtau;
            }

        fprintf(f1, "%g, %g \n", wdm/(2*3.142), cabs(v1[2]) );
    }
printf("Completed");
fclose(f1);

在迭代时间结束时，我将频率wdm的值和位移x = v1 [2]的最后一个值保存在一个文件中，然后继续用另一个频率进行时间迭代。因此，我的文件包含频率响应。

fprintf(f1, "%g, %g \n", wdm/(2*3.142), cabs(v1[2]) );

我使用了来自people.sc.fsu.edu/~jburkardt/c_src/rk4/rk4.html的runge kutta（rk4.c）并使用

修改了它的复杂数据类型

#include <complex.h>

以下是runge-kutta-4必须解决的函数：

/******************************************************************************/

double complex *rk4vec_f_ameya_complex_1 ( double t, int n, double complex u[] )

/******************************************************************************/

{
    double complex drive_m;
    double complex drive_c;
    double x;
    double xrf0_1;
    double complex *uprime;
    uprime = ( double * ) malloc ( 7 * sizeof ( double complex ) );

    //Check if memory unavailable
    if(uprime==NULL){
        printf("No memory available \n");
        return 0;
    }

    ///////////////////Second Order////////////////////////
    xrf0_1=xrf0*(1-exp(-0.2*t));
    drive_m=(xrf0*cexp(I*((u[6]-wm)/gammac)*t)/(2*wm*gammac));
    uprime[2]=u[3];
    uprime[3]=(wm/gammac)*(drive_m-u[3]*(2*I+2*gammam/wm)-u[2]*(2*I*gammam/gammac));

    return uprime;
    free(uprime);
}

如果我使用malloc（）导致内存不足，请提出任何解决方案。

Answer 1

一个明显的问题 - 您从uprime返回后免费rk4vec_f_ameya_complex_1，我希望您收到无法访问的代码警告。另外，uprime总是大小相同，为什么你会使用malloc？如果你只是把它变成一个大小为7的数组，那么你就不会遇到malloc的问题，你的代码可能运行得更快。

为什么以下代码内存不足？

1 个答案: