- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
我必须执行 Cannon's Algorithm这是一种并行算法,用于乘以维度为正方形的矩阵,维度可以被处理器数量的平方根整除。我写的这段代码运行得很好,但在生活中运行而不会崩溃只是故事的一半。它不会正确地将 A x B 乘以新矩阵 C。这是代码,请帮助我并引导我帮助我解决我做错的事情。很明显,这是家庭作业。
void shift_left(datatype** mat, int s, int row, int n, int amount) {
datatype* temp_buffer = malloc(sizeof(datatype) * n);
for(int col = 0; col < n; col++) {
datatype temp = mat[row][(col+amount)%s];
temp_buffer[(col+amount)%s] = mat[row][col];
temp_buffer[col] = temp;
}
memcpy(mat[row], temp_buffer, n);
free(temp_buffer);
}
void shift_up(datatype** mat, int s, int col, int n, int amount) {
datatype* temp_buffer = malloc(sizeof(datatype) * n);
for(int row = 0; row < n; row++) {
datatype temp = mat[(row+amount)%s][col];
temp_buffer[(row+amount)%s] = mat[row][col];
temp_buffer[row] = temp;
}
memcpy(&mat[0][col], temp_buffer, n);
free(temp_buffer);
}
void cannon_mul(int p_sqrt,datatype** a, datatype** b, datatype** c, int n) {
/* 2D matrices and n^2 sized only!*/
int i = 0, j = 0, k = 0;
int s = p_sqrt;
for(i = 0; i < (s-1); i++) {
shift_left(a, s, i, s-1, i); // Skew matrix a
}
for (i = 0; i < (s-1); i++) {
shift_up(b, s, i, s-1, i); // Skew matrix b
}
for(k = 0; k < (s-1); k++) {
for(i = 0; i < (s-1); i++) {
for(j = 0; j < (s-1); j++) {
c[i][j] += a[i][j]*b[i][j];
shift_left(a, s, i, s-1, 1);
shift_up(b, s, i, s-1, 1);
}
}
}
}
我的预感是移位错误或者我遗漏了算法的重要部分。我原来的 shift 函数没有使用临时缓冲区,所以这次我想使用临时缓冲区,但这并没有什么区别。如果有帮助,我可以显示一些示例输出,但结果不甚至与期望的结果相差甚远。好消息它运行得很快:)
1.48 0.14 9.47 8.99 8.06 0.06 6.68 1.04 4.44 7.50
7.26 8.87 2.21 6.27 2.12 7.91 0.65 5.24 0.45 4.94
0.47 4.13 1.87 2.25 6.83 1.52 6.41 9.14 9.22 8.91
7.34 2.70 6.78 2.78 3.51 4.95 5.27 0.85 9.51 6.82
0.28 6.73 0.70 8.88 7.14 9.09 2.36 5.38 6.43 9.00
7.13 6.71 6.92 9.81 5.13 9.35 7.50 5.16 4.68 3.62
1.30 6.26 4.55 4.27 0.51 2.23 3.19 8.75 6.57 9.07
7.49 6.41 1.04 7.78 7.16 2.78 2.25 6.23 9.42 0.32
3.21 3.60 2.04 2.93 4.29 3.88 2.78 8.01 4.57 6.47
7.52 3.77 0.63 5.97 7.32 4.90 9.63 4.90 8.46 1.90
用我的代码将上面的矩阵与自身相乘产生:
2.20 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
50.81 0.00 0.00 0.00 0.00 87.51 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
顺序程序产生这个结果:
163.41 212.17 144.32 227.10 251.03 205.60 245.63 277.33 368.99 334.11
257.85 230.82 203.60 314.08 246.02 240.12 228.37 197.90 264.38 228.24
234.13 272.10 110.75 294.84 263.16 242.07 209.54 316.13 339.23 260.51
185.33 215.59 192.26 283.31 270.80 208.38 265.08 291.49 312.24 319.73
313.23 301.95 182.04 348.11 283.20 337.49 266.54 284.57 355.28 281.07
293.25 323.29 281.35 393.92 325.24 313.62 313.48 342.95 418.37 401.91
255.88 238.25 122.17 254.52 243.58 204.49 217.69 273.03 314.89 214.45
219.26 239.07 200.18 309.98 262.21 242.68 190.02 245.85 297.96 308.56
209.03 213.11 126.24 266.48 233.88 199.33 193.28 228.92 277.50 202.27
210.31 264.67 227.59 337.79 261.40 250.35 225.77 295.00 331.92 352.17
重要的是要注意我正在展示我的程序的相关部分,如果您觉得我需要展示更多请这样做,我将提供更多代码。最后,为什么作业标签不见了?
一个人指出缓冲区太小并且缺少 sizeof 是一个愚蠢的错误,现在已更正。我试过了,我得到了相同的结果,所以显然问题与此无关。希望在 2 天内我可以开赏金来吸引一些人至少给我一个问题所在的线索。这是一个我似乎无法调试的错误,我必须承认这个算法的理解是零到零。我依赖的网络资源对增加我的理解作用甚微。
尝试使用 calloc 来零分配缓冲区,但它不会改变结果。很奇怪,但谢谢你的建议;我忘了内存没有分配归零。
我试过这个:
void shift_left(datatype** mat, int s, int row, int n, int amount) {
datatype* temp_buffer = calloc(n, sizeof(datatype) * n);
for(int col = 0; col < n; col++) {
/* temp_buffer[(col+amount)%s] = mat[row][col];
temp_buffer[col] = mat[row][(col+amount)%s]; */
temp_buffer[(col+amount)%s] = 0;
temp_buffer[col] = 0;
}
memcpy(mat[row], temp_buffer, sizeof(datatype) * n);
//free(temp_buffer);
}
void shift_up(datatype** mat, int s, int col, int n, int amount) {
datatype* temp_buffer = calloc(n, sizeof(datatype) * n);
for(int row = 0; row < n; row++) {
/* temp_buffer[(row+amount)%s] = mat[row][col];
temp_buffer[row] = mat[(row+amount)%s][col]; */
temp_buffer[(row+amount)%s] = 0;
temp_buffer[row] = 0;
}
memcpy(&mat[0][col], temp_buffer, sizeof(datatype) * n);
free(temp_buffer);
}
令人惊讶的是,结果是一样的。自从我注释了代码并将其替换为零后,它什么时候应该打印全零。我猜 memcpy 不工作。
我已经确认 memcpy 是罪魁祸首。但是为什么我不知道我被难住了,如果它有助于 datatype is just an alias for double 教授出于某种奇怪的原因写道,因为它确实无助于使代码更具可读性。
但如果我确实自己解决了它,我会很乐意展示我的问题的解决方案。
最佳答案
副本太小。
// memcpy(mat[row], temp_buffer, n);
memcpy(mat[row], temp_buffer, n * sizeof(datatype) );
同样适用于 memcpy(&mat[0][col], temp_buffer, n);
关于c - 实现 Cannon 的算法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/29637513/
滑动窗口限流 滑动窗口限流是一种常用的限流算法,通过维护一个固定大小的窗口,在单位时间内允许通过的请求次数不超过设定的阈值。具体来说,滑动窗口限流算法通常包括以下几个步骤: 初始化:设置窗口
表达式求值:一个只有+,-,*,/的表达式,没有括号 一种神奇的做法:使用数组存储数字和运算符,先把优先级别高的乘法和除法计算出来,再计算加法和减法 int GetVal(string s){
【算法】前缀和 题目 先来看一道题目:(前缀和模板题) 已知一个数组A[],现在想要求出其中一些数字的和。 输入格式: 先是整数N,M,表示一共有N个数字,有M组询问 接下来有N个数,表示A[1]..
1.前序遍历 根-左-右的顺序遍历,可以使用递归 void preOrder(Node *u){ if(u==NULL)return; printf("%d ",u->val);
先看题目 物品不能分隔,必须全部取走或者留下,因此称为01背包 (只有不取和取两种状态) 看第一个样例 我们需要把4个物品装入一个容量为10的背包 我们可以简化问题,从小到大入手分析 weightva
我最近在一次采访中遇到了这个问题: 给出以下矩阵: [[ R R R R R R], [ R B B B R R], [ B R R R B B], [ R B R R R R]] 找出是否有任
我正在尝试通过 C++ 算法从我的 outlook 帐户发送一封电子邮件,该帐户已经打开并记录,但真的不知道从哪里开始(对于 outlook-c++ 集成),谷歌也没有帮我这么多。任何提示将不胜感激。
我发现自己像这样编写了一个手工制作的 while 循环: std::list foo; // In my case, map, but list is simpler auto currentPoin
我有用于检测正方形的 opencv 代码。现在我想在检测正方形后,代码运行另一个命令。 代码如下: #include "cv.h" #include "cxcore.h" #include "high
我正在尝试模拟一个 matlab 函数“imfill”来填充二进制图像(1 和 0 的二维矩阵)。 我想在矩阵中指定一个起点,并像 imfill 的 4 连接版本那样进行洪水填充。 这是否已经存在于
我正在阅读 Robert Sedgewick 的《C++ 算法》。 Basic recurrences section it was mentioned as 这种循环出现在循环输入以消除一个项目的递
我正在思考如何在我的日历中生成代表任务的数据结构(仅供我个人使用)。我有来自 DBMS 的按日期排序的任务记录,如下所示: 买牛奶(18.1.2013) 任务日期 (2013-01-15) 任务标签(
输入一个未排序的整数数组A[1..n]只有 O(d) :(d int) 计算每个元素在单次迭代中出现在列表中的次数。 map 是balanced Binary Search Tree基于确保 O(nl
我遇到了一个问题,但我仍然不知道如何解决。我想出了如何用蛮力的方式来做到这一点,但是当有成千上万的元素时它就不起作用了。 Problem: Say you are given the followin
我有一个列表列表。 L1= [[...][...][.......].......]如果我在展平列表后获取所有元素并从中提取唯一值,那么我会得到一个列表 L2。我有另一个列表 L3,它是 L2 的某个
我们得到二维矩阵数组(假设长度为 i 和宽度为 j)和整数 k我们必须找到包含这个或更大总和的最小矩形的大小F.e k=7 4 1 1 1 1 1 4 4 Anwser是2,因为4+4=8 >= 7,
我实行 3 类倒制,每周换类。顺序为早类 (m)、晚类 (n) 和下午类 (a)。我固定的订单,即它永远不会改变,即使那个星期不工作也是如此。 我创建了一个函数来获取 ISO 周数。当我给它一个日期时
假设我们有一个输入,它是一个元素列表: {a, b, c, d, e, f} 还有不同的集合,可能包含这些元素的任意组合,也可能包含不在输入列表中的其他元素: A:{e,f} B:{d,f,a} C:
我有一个子集算法,可以找到给定集合的所有子集。原始集合的问题在于它是一个不断增长的集合,如果向其中添加元素,我需要再次重新计算它的子集。 有没有一种方法可以优化子集算法,该算法可以从最后一个计算点重新
我有一个包含 100 万个符号及其预期频率的表格。 我想通过为每个符号分配一个唯一(且前缀唯一)的可变长度位串来压缩这些符号的序列,然后将它们连接在一起以表示序列。 我想分配这些位串,以使编码序列的预
我是一名优秀的程序员,十分优秀!