如果我运行一个循环一百万次,我是否必须担心在每次迭代中声明双精度?

时间:2018-03-24 19:15:46

标签: c memory-management

我循环了一百万次。在循环中,我调用C函数进行一些数学运算(确切地说,从各种分布生成随机变量)。作为该函数的一部分,我声明了几个double变量来保存转换的一部分。一个例子:

void getRandNorm(double *randnorm, double mean, double var, int n) 
{
    // Declare variables
    double u1;
    double u2;
    int arrptr = 0;
    double sigma = sqrt(var); // the standard deviation

    while (arrptr < n) { 
        // Generate two uniform random variables
        u1 = rand() / (double)RAND_MAX;
        u2 = rand() / (double)RAND_MAX;

        // Box-Muller transform
        randnorm[arrptr] = sqrt(-2*log(u1))*cos(2*pi*u2)*sigma+mean;
        arrptr++;
        if (arrptr < n) { // for an odd n, we cannot add off the end
            randnorm[arrptr] = sqrt(-2*log(u2))*cos(2*pi*u1)*sigma+mean;
            arrptr++;
        }
    }   
}

调用循环:

iter = 1000000 // or something
for (i = 0; i < iter; i++) {

    // lots of if statements
    getRandNorm(sample1, truemean1, truevar1, n);

    // some more analysis
}

我正在努力加快运行时间。在我看来,我不知道我声明的所有这些double变量发生了什么。我假设为一百万个循环中的每个循环分配一个新的8字节内存块用于double。所有这些内存位置会发生什么?它们在C函数中声明;他们能活下来吗?在脚本退出之前,它们是否仍处于锁定状态?

这个问题的上下文是将这个C程序包装成python函数。如果我要从python并行多次执行此函数,我想确保我尽可能节省内存。

3 个答案:

答案 0 :(得分:3)

如果你在谈论这样的事情:

for(int i=0;i<100000;i++){
    double d = 5;
    // some other stuff here
}

d仅由编译器分配一次。它主要等同于在for循环之上声明它,除了范围没有延伸到目前为止。

但是,如果您正在做这样的事情:

for(int i=0;i<1000000;i++){
    double *d = malloc(sizeof(double));
    free(d);
}

然后是的,您将分配double 100万次,但它可能会重新使用内存进行后续分配。最后,如果你在我的第二个例子中没有释放内存,你将泄漏16-32MB的内存。

答案 1 :(得分:2)

简短的回答是:,如果你在C中的循环中声明这些module ApplicationCable class Connection < ActionCable::Connection::Base identified_by :current_user def connected self.current_user = find_verified_user logger.add_tags "ActionCable", "User #{current_user.id}" end protected def find_verified_user if current_user = env['warden'].user current_user else reject_unauthorized_connection end end end end class ChatroomsChannel < ApplicationCable::Channel def subscribed current_user.chatrooms.each do |chatroom| stream_from "chatrooms:#{chatroom.id}" end end def unsubscribed stop_all_streams end def send_message(data) Rails.logger.info data @chatroom = Chatroom.find(data["chatroom_id"]) message = @chatroom.messages.create(body: data["body"], user: current_user) MessageRelayJob.perform_later(message) end end 变量应该没关系。通过double变量,我假设你的意思是变量类型为double

答案很长:请发布您的代码,以便人们可以告诉您是否做错了什么以及如何解决它以提高正确性和/或性能(广泛的主题)。

最后的答案是:使用提供的代码,无论是在循环体内还是外部声明doubleu1都没有区别。一个好的编译器可能会生成相同的代码。

通过仅测试一次奇数情况,您可以稍微改进一下代码:

u2

注意:void getRandNorm(double *randnorm, double mean, double var, int n, double pi) { // Declare variables double u1, u2; double sigma = sqrt(var); // the standard deviation int arrptr, odd; odd = n & 1; // check if n is odd n -= odd; // make n even for (arrptr = 0; arrptr < n; arrptr += 2) { // Generate two uniform random variables u1 = rand() / (double)RAND_MAX; u2 = rand() / (double)RAND_MAX; // Box-Muller transform randnorm[arrptr + 0] = sqrt(-2*log(u1)) * cos(2*pi*u2) * sigma + mean; randnorm[arrptr + 1] = sqrt(-2*log(u2)) * cos(2*pi*u1) * sigma + mean; } if (odd) { u1 = rand() / (double)RAND_MAX; u2 = rand() / (double)RAND_MAX; randnorm[arrptr++] = sqrt(-2*log(u1)) * cos(2*pi*u2) * sigma + mean; } } 用于对称,编译器不会为此添加生成任何代码。

答案 2 :(得分:1)

关于你的问题:如果我运行一个循环一百万次,我是否必须担心在每次迭代中声明双精度?

变量正在堆栈中声明。所以他们消失了#39;当函数退出时。下一次执行该功能会重新创建&#39;变量,所以(实际上)只有一个变量实例,即便如此,只有在函数执行时才会出现。

因此,调用该函数的次数并不重要。