我循环了一百万次。在循环中,我调用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并行多次执行此函数,我想确保我尽可能节省内存。
答案 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
。
答案很长:请发布您的代码,以便人们可以告诉您是否做错了什么以及如何解决它以提高正确性和/或性能(广泛的主题)。
最后的答案是:使用提供的代码,无论是在循环体内还是外部声明double
和u1
都没有区别。一个好的编译器可能会生成相同的代码。
通过仅测试一次奇数情况,您可以稍微改进一下代码:
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;变量,所以(实际上)只有一个变量实例,即便如此,只有在函数执行时才会出现。
因此,调用该函数的次数并不重要。