- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
请参见我的以下代码片段(floatalloc2用于分配数据类型为float的2D连续数组,如果有兴趣,请参见附录):
#include <stdio.h>
#include <stdlib.h>
#include <mpi.h>
int main(int argc, char *argv[])
{
float **p=NULL, **buffer=NULL;
int it, nt=3, i, j, k, NP, MYID, nx=1, nz=2, nsrc=3, isrc;
MPI_Init ( &argc, &argv );
MPI_Comm_size ( MPI_COMM_WORLD, &NP );
MPI_Comm_rank ( MPI_COMM_WORLD, &MYID );
p = floatalloc2(nx,nz);
memset(p[0],0,nz*nx*sizeof(float));
for (it=0; it<nt; it++){
for (isrc=MYID; isrc<nsrc; isrc+=NP){
for (j=0; j<nz; j++){
for (i=0; i<nx; i++){
p[j][i] += 1.5 + (float)(isrc) + (float)(j);
}
}
}
for (k=0;k<nsrc-1;k++){
if (MYID==k){
buffer = floatalloc2(nx,nz);
memset(buffer[0],0,nz*nx*sizeof(float));
buffer = p;
}else{
buffer = floatalloc2(nx,nz);
memset(buffer[0],0,nz*nx*sizeof(float));
}
MPI_Barrier(MPI_COMM_WORLD);
MPI_Bcast(&buffer[0][0],nx*nz,MPI_FLOAT,k,MPI_COMM_WORLD);
MPI_Barrier(MPI_COMM_WORLD);
for (j=0; j<nz; j++){
for (i=0; i<nx; i++){
printf("it=%d,k=%d,Node %d,p[%d][%d]=%f\n",it,k,MYID,j,i,p[j][i]);
}
}
free(*buffer);free(buffer); /*w/o this line is ok while not ok with this line */
}
}
MPI_Finalize();
exit(0);
}
free(*buffer);free(buffer);
,那是相当不错的。但是,如果使用
free
,则结果是错误的。那么我的代码有什么问题呢?
/*@out@*/ void *sf_alloc (size_t n, size_t size )
/*< output-checking allocation >*/
{
void *ptr;
size *= n;
if (0>=size) sf_error("%s: illegal allocation (%d bytes)",__FILE__,size);
ptr = malloc (size);
if (NULL == ptr)
sf_error ("%s: cannot allocate %lu bytes:", __FILE__,size);
return ptr;
}
/*@out@*/ float *sf_floatalloc (size_t n)
/*< float allocation >*/
{
float *ptr;
ptr = (float*) sf_alloc (n,sizeof(float));
return ptr;
}
/*@out@*/ float **floatalloc2 (size_t n1 , size_t n2 )
/*< float 2-D allocation, out[0] points to a contiguous array >*/
{
size_t i2;
float **ptr;
ptr = (float**) sf_alloc (n2,sizeof(float*));
ptr[0] = sf_floatalloc (n1*n2);
for (i2=1; i2 < n2; i2++) {
ptr[i2] = ptr[0]+i2*n1;
}
return ptr;
}
最佳答案
假设我们使用float **p = floatalloc2(columns, rows)
。然后,为p
浮点指针(rows
)分配float *
,为p[0]
浮点分配columns*rows
。如果该库不提供floatfree2()
或free2()
函数,则可以使用free(*p); free(p);
以此顺序释放为2D数组分配的内存。这也是OP所做的。
内循环
for (k=0;k<nsrc-1;k++){
if (MYID==k){
buffer = floatalloc2(nx,nz);
memset(buffer[0],0,nz*nx*sizeof(float));
buffer = p;
}else{
buffer = floatalloc2(nx,nz);
memset(buffer[0],0,nz*nx*sizeof(float));
}
...
free(*buffer);free(buffer); /*w/o this line is ok while not ok with this line */
}
buffer = p;
行替换了前两行分配和初始化的2D数组;这本质上是内存泄漏。我想这个想法是将
p
的内容复制到它,但是我不确定。
free(*buffer); free(buffer);
行最终释放了
p
迭代中
k == MYID
描述的2D数组。
buffer
为
p
的别名。
for (k = 0; k < nsrc-1; k++) {
if (MYID == k) {
buffer = p;
} else {
buffer = floatalloc2(nx,nz);
memset(buffer[0], 0, nz*nx*sizeof(float));
}
...
if (buffer != p) {
free(*buffer);
free(buffer);
}
}
k == MYID
时,
k == MYID
实际上指向
buffer
;否则,将为每个迭代动态分配和释放它。显然,当
p
为
buffer
别名时,我们不会释放它,因为那样会释放
p
。
p
之前添加
free(*p); free(p);
,以释放为
MPI_Finalize()
分配的内存,但这并不是绝对必要的,因为该过程将要退出。 (但是,如果您使用例如Valgrind查找内存泄漏,或者希望表明您(程序员)确实正确跟踪动态分配,则很有用;在这种情况下,添加注释可能会很有用。)
float **p, **buffer;
p = floatalloc2(nx, nz);
memset(p[0], 0, nz*nx*sizeof p[0][0]);
buffer = floatalloc2(nx, nz);
memset(buffer[0], 0, nz*nx*sizeof buffer[0][0]);
for (it=0; it < nt; it++) {
for (isrc=MYID; isrc < nsrc; isrc+=NP) {
for (j=0; j<nz; j++){
for (i=0; i<nx; i++){
p[j][i] += 1.5 + (float)(isrc) + (float)(j);
}
}
}
for (k=0; k<nsrc-1; k++) {
float **data;
if (MYID == k) {
data = p;
} else {
data = buffer;
memset(buffer[0], 0, nz*nx*sizeof buffer[0][0]);
}
MPI_Barrier(MPI_COMM_WORLD);
MPI_Bcast(&(data[0][0]), nx*nz, MPI_FLOAT, k, MPI_COMM_WORLD);
MPI_Barrier(MPI_COMM_WORLD);
for (j=0; j<nz; j++) {
for (i=0; i<nx; i++) {
printf("it=%d,k=%d,Node %d,data[%d][%d]=%f, p[%d][%d]=%f\n",
it, k, MYID, j, i, data[j][i], p[j][i]);
}
}
}
}
free(*buffer);
free(buffer);
free(*p);
free(p);
p
和
sizeof p[0][0]
是评估为
sizeof buffer[0][0]
/
p
每个元素的大小(以字符为单位)的语句,并且实际上不检查任何内存。即使编译器
buffer
或
p == NULL
,这些也是允许的并且是安全的,因为编译器仅检查表达式的类型(在
p[0] == NULL
运算符的右侧)。为了提醒自己,我从未在
sizeof
运算符右侧的变量引用周围使用括号。 (当您使用诸如
sizeof
之类的类型时,必须加上括号。然后,我将其写为
float
。)这可以帮助我记住
sizeof (float)
是一个运算符,尽管它的行为不像一个函数。可能看起来像一个函数。
sizeof
只是对
float **data;
或
p
中实际数据的引用(或别名);我们需要使用的只是分配它,
buffer
或
data = p;
。因为它是指向完全相同的数据的引用或别名,所以我们根本不对它进行抄送。如果这样做,我们最终将释放原始数据
data = buffer;
或
free()
。
p
,其中
buffer
是元素,并且
float **data
和
data[row][column]
都是指针)不是很有效。在实践中,使用结构代替。例如:
typedef struct {
int rows;
int cols;
float *data;
} float2d;
#define FLOAT2D_INIT { 0, 0, NULL }
static inline void float2d_alloc(float2d *m, const int rows, const int cols)
{
do {
if (!m) {
fprintf(stderr, "float2d_alloc(): No matrix specified.\n");
break;
}
if (rows < 1 || cols < 1) {
fprintf(stderr, "float2d_alloc(): Invalid matrix size (%d rows, %d cols).\n", rows, cols);
break;
}
m->data = malloc((size_t)rows * (size_t)cols * sizeof m->data[0]);
if (!m->data) {
fprintf(stderr, "float2d_alloc(): Not enough memory available.\n");
break;
}
m->rows = rows;
m->cols = cols;
return;
} while (0);
/* MPI_Abort(MPI_COMM_WORLD, 1); */
exit(1);
}
static inline void float2d_free(float2d *m)
{
if (m) {
free(m->data);
m->rows = 0;
m->cols = 0;
m->data = NULL;
}
}
data[row]
,则要使用
data
来访问行
float2d p;
,列
r
上的元素。所有数据在内存中都是连续的;有
c
个元素,总大小为
p.data[r*p.cols + c]
。要引用行
p.cols*p.rows
上的浮点数数组,可以使用
p.cols*p.rows*sizeof p.data[0]
(等效于
r
)。
关于c - 在循环中分配和释放内存(C + MPI),我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/51788667/
我在具有 2CPU 和 3.75GB 内存 (https://aws.amazon.com/ec2/instance-types/) 的 c3.large Amazon EC2 ubuntu 机器上运
我想通过用户空间中的mmap-ing并将地址发送到内核空间从用户空间写入VGA内存(视频内存,而不是缓冲区),我将使用pfn remap将这些mmap-ed地址映射到vga内存(我将通过 lspci
在 Mathematica 中,如果你想让一个函数记住它的值,它在语法上是很轻松的。例如,这是标准示例 - 斐波那契: fib[1] = 1 fib[2] = 1 fib[n_]:= fib[n] =
我读到动态内存是在运行时在堆上分配的,而静态内存是在编译时在堆栈上分配的,因为编译器知道在编译时必须分配多少内存。 考虑以下代码: int n; cin>>n; int a[n]; 如果仅在运行期间读
我是 Python 的新手,但我之前还不知道这一点。我在 for 循环中有一个基本程序,它从站点请求数据并将其保存到文本文件但是当我检查我的任务管理器时,我发现内存使用量只增加了?长时间运行时,这对我
我正在设计一组数学函数并在 CPU 和 GPU(使用 CUDA)版本中实现它们。 其中一些函数基于查找表。大多数表占用 4KB,其中一些占用更多。基于查找表的函数接受一个输入,选择查找表的一两个条目,
读入一个文件,内存被动态分配给一个字符串,文件内容将被放置在这里。这是在函数内部完成的,字符串作为 char **str 传递。 使用 gdb 我发现在行 **(str+i) = fgetc(aFil
我需要证实一个理论。我正在学习 JSP/Java。 在查看了一个现有的应用程序(我没有写)之后,我注意到一些我认为导致我们的性能问题的东西。或者至少是其中的一部分。 它是这样工作的: 1)用户打开搜索
n我想使用memoization缓存某些昂贵操作的结果,这样就不会一遍又一遍地计算它们。 两个memoise和 R.cache适合我的需要。但是,我发现缓存在调用之间并不可靠。 这是一个演示我看到的问
我目前正在分析一些 javascript shell 代码。这是该脚本中的一行: function having() { memory = memory; setTimeout("F0
我有一种情况,我想一次查询数据库,然后再将整个数据缓存在内存中。 我得到了内存中 Elasticsearch 的建议,我用谷歌搜索了它是什么,以及如何在自己的 spring boot 应用程序中实现它
我正在研究 Project Euler (http://projecteuler.net/problem=14) 的第 14 题。我正在尝试使用内存功能,以便将给定数字的序列长度保存为部分结果。我正在
所以,我一直在做 Java 内存/注意力游戏作业。我还没有达到我想要的程度,它只完成了一半,但我确实让 GUI 大部分工作了......直到我尝试向我的框架添加单选按钮。我认为问题可能是因为我将 JF
我一直在尝试使用 Flask-Cache 的 memoize 功能来仅返回 statusTS() 的缓存结果,除非在另一个请求中满足特定条件,然后删除缓存。 但它并没有被删除,并且 Jinja 模板仍
我对如何使用 & 运算符来减少内存感到非常困惑。 我可以回答下面的问题吗? clase C{ function B(&$a){ $this->a = &$a; $thi
在编写代码时,我遇到了一个有趣的问题。 我有一个 PersonPOJO,其 name 作为其 String 成员之一及其 getter 和 setter class PersonPOJO { priv
在此代码中 public class Base { int length, breadth, height; Base(int l, int b, int h) { l
Definition Structure padding is the process of aligning data members of the structure in accordance
在 JavaScript Ninja 的 secret 中,作者提出了以下方案,用于在没有闭包的情况下内存函数结果。他们通过利用函数是对象这一事实并在函数上定义一个属性来存储过去调用函数的结果来实现这
我正在尝试找出 map 消耗的 RAM 量。所以,我做了以下事情;- Map cr = crPair.collectAsMap(); // 200+ entries System.out.printl
我是一名优秀的程序员,十分优秀!