- r - 以节省内存的方式增长 data.frame
- ruby-on-rails - ruby/ruby on rails 内存泄漏检测
- android - 无法解析导入android.support.v7.app
- UNIX 域套接字与共享内存(映射文件)
我正在实现求解线性系统的多线程 Jordan-Gauss 方法,我发现在两个线程上运行只比在单线程上运行少大约 15% 的时间,而不是理想的 50%。所以我写了一个简单的程序来复制这个。在这里,我创建了一个 2000x2000 矩阵,并为每个线程提供了 2000/THREADS_NUM 行,以便用它们进行一些计算。
#include <stdlib.h>
#include <stdio.h>
#include <pthread.h>
#include <time.h>
#ifndef THREADS_NUM
#define THREADS_NUM 1
#endif
#define MATRIX_SIZE 2000
typedef struct {
double *a;
int row_length;
int rows_number;
} TWorkerParams;
void *worker_thread(void *params_v)
{
TWorkerParams *params = (TWorkerParams *)params_v;
int row_length = params->row_length;
int i, j, k;
int rows_number = params->rows_number;
double *a = params->a;
for(i = 0; i < row_length; ++i) // row_length is always the same
{
for(j = 0; j < rows_number; ++j) // rows_number is inverse proportional
// to the number of threads
{
for(k = i; k < row_length; ++k) // row_length is always the same
{
a[j*row_length + k] -= 2.;
}
}
}
return NULL;
}
int main(int argc, char *argv[])
{
// The matrix is of size NxN
double *a =
(double *)malloc(MATRIX_SIZE * MATRIX_SIZE * sizeof(double));
TWorkerParams *params =
(TWorkerParams *)malloc(THREADS_NUM * sizeof(TWorkerParams));
pthread_t *workers = (pthread_t *)malloc(THREADS_NUM * sizeof(pthread_t));
struct timespec start_time, end_time;
int rows_per_worker = MATRIX_SIZE / THREADS_NUM;
int i;
if(!a || !params || !workers)
{
fprintf(stderr, "Error allocating memory\n");
return 1;
}
for(i = 0; i < MATRIX_SIZE*MATRIX_SIZE; ++i)
a[i] = 4. * i; // just an example matrix
// Initializtion of matrix is done, now initialize threads' params
for(i = 0; i < THREADS_NUM; ++i)
{
params[i].a = a + i * rows_per_worker * MATRIX_SIZE;
params[i].row_length = MATRIX_SIZE;
params[i].rows_number = rows_per_worker;
}
// Get start time
clock_gettime(CLOCK_MONOTONIC, &start_time);
// Create threads
for(i = 0; i < THREADS_NUM; ++i)
{
if(pthread_create(workers + i, NULL, worker_thread, params + i))
{
fprintf(stderr, "Error creating thread\n");
return 1;
}
}
// Join threads
for(i = 0; i < THREADS_NUM; ++i)
{
if(pthread_join(workers[i], NULL))
{
fprintf(stderr, "Error creating thread\n");
return 1;
}
}
clock_gettime(CLOCK_MONOTONIC, &end_time);
printf("Duration: %lf msec.\n", (end_time.tv_sec - start_time.tv_sec)*1e3 +
(end_time.tv_nsec - start_time.tv_nsec)*1e-6);
return 0;
}
这里是我如何编译它的:
gcc threads_test.c -o threads_test1 -lrt -pthread -DTHREADS_NUM=1 -Wall -Werror -Ofast
gcc threads_test.c -o threads_test2 -lrt -pthread -DTHREADS_NUM=2 -Wall -Werror -Ofast
现在,当我运行时,我得到:
./threads_test1
Duration: 3695.359552 msec.
./threads_test2
Duration: 3211.236612 msec.
这意味着双线程程序的运行速度比单线程快 13%,即使线程之间没有同步并且它们不共享任何内存也是如此。我找到了这个答案:https://stackoverflow.com/a/14812411/5647501并认为这可能是处理器缓存的一些问题,所以我添加了填充,但结果仍然保持不变。我更改了我的代码如下:
typedef struct {
double *a;
int row_length;
int rows_number;
volatile char padding[64 - 2*sizeof(int)-sizeof(double)];
} TWorkerParams;
#define VAR_SIZE (sizeof(int)*5 + sizeof(double)*2)
#define MEM_SIZE ((VAR_SIZE / 64 + 1) * 64 )
void *worker_thread(void *params_v)
{
TWorkerParams *params = (TWorkerParams *)params_v;
volatile char memory[MEM_SIZE];
int *row_length = (int *)(memory + 0);
int *i = (int *)(memory + sizeof(int)*1);
int *j = (int *)(memory + sizeof(int)*2);
int *k = (int *)(memory + sizeof(int)*3);
int *rows_number = (int *)(memory + sizeof(int)*4);
double **a = (double **)(memory + sizeof(int)*5);
*row_length = params->row_length;
*rows_number = params->rows_number;
*a = params->a;
for(*i = 0; *i < *row_length; ++*i) // row_length is always the same
{
for(*j = 0; *j < *rows_number; ++*j) // rows_number is inverse proportional
// to the number of threads
{
for(*k = 0; *k < *row_length; ++*k) // row_length is always the same
{
(*a + *j * *row_length)[*k] -= 2. * *k;
}
}
}
return NULL;
}
所以我的问题是:为什么我在这里使用两个线程时只能获得 15% 的加速而不是 50%?任何帮助或建议将不胜感激。我正在运行 64 位 Ubuntu Linux,内核 3.19.0-39-generic,CPU Intel Core i5 4200M(具有多线程的两个物理内核),但我也在另外两台机器上进行了测试,结果相同。
编辑:如果我用 a[0] -= 2.;
替换 a[j*row_length + k] -= 2.;
,我得到预期的加速:
./threads_test1
Duration: 1823.689481 msec.
./threads_test2
Duration: 949.745232 msec.
编辑 2:现在,当我用 a[k] -= 2.;
替换它时,我得到以下信息:
./threads_test1
Duration: 1039.666979 msec.
./threads_test2
Duration: 1323.460080 msec.
这个我根本买不到。
最佳答案
这是一个经典问题,交换 i 和 j for 循环。
您首先遍历列,然后在内部循环中处理行,这意味着您有比必要更多的缓存未命中。
我的原始代码结果(第一个没有填充的版本):
$ ./matrix_test1
Duration: 4620.799763 msec.
$ ./matrix_test2
Duration: 2800.486895 msec.
(实际上比你的改进更好)
在为 i 和 j 切换 for 循环之后:
$ ./matrix_test1
Duration: 1450.037651 msec.
$ ./matrix_test2
Duration: 728.690853 msec.
这里是 2 倍的加速。
编辑:事实上,原来的并不坏,因为 k 索引仍然遍历行迭代列,但在外循环中迭代行仍然要好得多。当 i 上升时,您在最内层循环中处理的项目越来越少,所以它仍然很重要。
EDIT2:(删除了 block 解决方案,因为它实际上产生了不同的结果)——但仍然应该可以利用 block 来提高缓存性能。
关于c++ - 使用多线程时性能几乎没有提高,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/34137578/
我将 Bootstrap 与 css 和 java 脚本结合使用。在不影响前端代码的情况下,我真的很难在css中绘制这个背景。在许多问题中,人们将宽度和高度设置为 0%。但是由于我的导航栏,我不能使用
我正在用 c 编写一个程序来读取文件的内容。代码如下: #include void main() { char line[90]; while(scanf("%79[^\
我想使用 javascript 获取矩阵数组的所有对 Angular 线。假设输入输出如下: input = [ [1,2,3], [4,5,6], [7,8,9], ] output =
可以用pdfmake绘制lines,circles和other shapes吗?如果是,是否有documentation或样本?我想用jsPDF替换pdfmake。 最佳答案 是的,有可能。 pdfm
我有一个小svg小部件,其目的是显示角度列表(参见图片)。 现在,角度是线元素,仅具有笔触,没有填充。但是现在我想使用一种“内部填充”颜色和一种“笔触/边框”颜色。我猜想line元素不能解决这个问题,
我正在为带有三角对象的 3D 场景编写一个非常基本的光线转换器,一切都工作正常,直到我决定尝试从场景原点 (0/0/0) 以外的点转换光线。 但是,当我将光线原点更改为 (0/1/0) 时,相交测试突
这个问题已经有答案了: Why do people write "#!/usr/bin/env python" on the first line of a Python script? (22 个回
如何使用大约 50 个星号 * 并使用 for 循环绘制一条水平线?当我尝试这样做时,结果是垂直(而不是水平)列出 50 个星号。 public void drawAstline() { f
这是一个让球以对角线方式下降的 UI,但球保持静止;线程似乎无法正常工作。你能告诉我如何让球移动吗? 请下载一个球并更改目录,以便程序可以找到您的球的分配位置。没有必要下载足球场,但如果您愿意,也可以
我在我的一个项目中使用 Jmeter 和 Ant,当我们生成报告时,它会在报告中显示 URL、#Samples、失败、成功率、平均时间、最短时间、最长时间。 我也想在报告中包含 90% 的时间线。 现
我有一个不寻常的问题,希望有人能帮助我。我想用 Canvas (android) 画一条 Swing 或波浪线,但我不知道该怎么做。它将成为蝌蚪的尾部,所以理想情况下我希望它的形状更像三角形,一端更大
这个问题已经有答案了: Checking Collision of Shapes with JavaFX (1 个回答) 已关闭 8 年前。 我正在使用 JavaFx 8 库。 我的任务很简单:我想检
如何按编号的百分比拆分文件。行数? 假设我想将我的文件分成 3 个部分(60%/20%/20% 部分),我可以手动执行此操作,-_-: $ wc -l brown.txt 57339 brown.tx
我正在努力实现这样的目标: 但这就是我设法做到的。 你能帮我实现预期的结果吗? 更新: 如果我删除 bootstrap.css 依赖项,问题就会消失。我怎样才能让它与 Bootstrap 一起工作?
我目前正在构建一个网站,但遇到了 transform: scale 的问题。我有一个按钮,当用户将鼠标悬停在它上面时,会发生两件事: 背景以对 Angular 线“扫过” 按钮标签颜色改变 按钮稍微变
我需要使用直线和仿射变换绘制大量数据点的图形(缩放图形以适合 View )。 目前,我正在使用 NSBezierPath,但我认为它效率很低(因为点在绘制之前被复制到贝塞尔路径)。通过将我的数据切割成
我正在使用基于 SVM 分类的 HOG 特征检测器。我可以成功提取车牌,但提取的车牌除了车牌号外还有一些不必要的像素/线。我的图像处理流程如下: 在灰度图像上应用 HOG 检测器 裁剪检测到的区域 调
我有以下图片: 我想填充它的轮廓(即我想在这张图片中填充线条)。 我尝试了形态学闭合,但使用大小为 3x3 的矩形内核和 10 迭代并没有填满整个边界。我还尝试了一个 21x21 内核和 1 迭代,但
我必须找到一种算法,可以找到两组数组之间的交集总数,而其中一个数组已排序。 举个例子,我们有这两个数组,我们向相应的数字画直线。 这两个数组为我们提供了总共 7 个交集。 有什么样的算法可以帮助我解决
简单地说 - 我想使用透视投影从近裁剪平面绘制一条射线/线到远裁剪平面。我有我认为是使用各种 OpenGL/图形编程指南中描述的方法通过单击鼠标生成的正确标准化的世界坐标。 我遇到的问题是我的光线似乎
我是一名优秀的程序员,十分优秀!