- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
是否有更好的解决方案来快速初始化 C 数组(在堆上创建)?就像我们使用大括号一样
double** matrix_multiply(const double **l_matrix, const double **r_matrix);
foo() {
double DCT_matrix[8][8] = {
{ 0.3536, 0.3536, 0.3536, 0.3536, 0.3536, 0.3536, 0.3536, 0.3536 },
{ 0.4904, 0.4157, 0.2778, 0.0975, -0.0975, -0.2778, -0.4157, -0.4904 },
{ 0.4619, 0.1913, -0.1913, -0.4619, -0.4619, -0.1913, 0.1913, 0.4619 },
{ 0.4157, -0.0975, -0.4904, -0.2778, 0.2778, 0.4904, 0.0975, -0.4157 },
{ 0.3536, -0.3536, -0.3536, 0.3536, 0.3536, -0.3536, -0.3536, 0.3536 },
{ 0.2778, -0.4904, 0.0975, 0.4157, -0.4157, -0.0975, 0.4904, -0.2778 },
{ 0.1913, -0.4619, 0.4619, -0.1913, -0.1913, 0.4619, -0.4619, 0.1913 },
{ 0.0975, -0.2778, 0.4157, -0.4904, 0.4904, -0.4157, 0.2778, -0.0975 }
};
const double other_matrix[8][8] = {
{ 26, -5, -5, -5, -5, -5, -5, 8 },
{ 64, 52, 8, 26, 26, 26, 8, -18 },
{ 126, 70, 26, 26, 52, 26, -5, -5 },
{ 111, 52, 8, 52, 52, 38, -5, -5 },
{ 52, 26, 8, 39, 38, 21, 8, 8 },
{ 0, 8, -5, 8, 26, 52, 70, 26 },
{ -5, -23, -18, 21, 8, 8, 52, 38 },
{ -18, 8, -5, -5, -5, 8, 26, 8 }
};
matrix_multiply(DCT_matrix, other_matrix); // Segfault
}
最佳答案
在评论部分,很明显主要缺陷是matrix_multiply
的设计。没有理由将 double**
传递给它。
首先,double**
的布局表明矩阵是这样创建的:
double **mat = malloc(8 * sizeof *mat);
for(int i=0; i<8; i++)
mat[i] = malloc(8 * sizeof *mat[0]);
这就是入门类(class)中学习如何处理指针时的教学方式。但这会导致许多不必要的 malloc 调用,这也会使释放所有内存变得笨拙。最重要的是,它很慢。在分配/释放和使用内存时都如此,因为矩阵不适合缓存。 Here is an answer I wrote that is about the cache
我还有其他几个here
所以,在你的函数中,我假设你会有很多像l_matrix[x][y]
这样的东西。从 double**
到 double*
时,将其替换为 l_matrix[x + 8*y]
。
最重要的是,由于您没有输出参数,我假设您在函数内部分配输出。不要那样做。在很多情况下,您想要重用为矩阵分配的内存。相反,提供一个输出参数并在函数外部分配内存。这给了我们签名:
void matrix_multiply(const double *A, const double *B, double *output)
此外,由于在您的情况下矩阵大小是固定的,我建议以某种反射(reflect)这一点的方式重命名该函数。也许是 DCT_matrix_multiply
但这取决于你。
如果需要,您可以将返回类型更改为 double*
并以 return output;
结束函数。这不是强制性的,但它提供了执行此操作的选项:
matrix_multiply(matrix_multiply(A, B, C),
matrix_multiply(A, D, E),
F);
但是,如果这样做,请务必小心,因为未指定参数的求值顺序。
如果您确实想要一个执行分配的函数,请编写一个包装器:
double *wrapper(const double *A, const double *B) {
double *output = malloc(8*8 * sizeof *output);
if(!output) return NULL;
return matrix_multiply(A, B, output);
}
如果您想使用大小不固定的矩阵,则有点棘手,但总的来说,我仍然建议传递带有有关大小的附加信息的 double*
。像这样的东西:
double *mul(const double *A, const double *B, double *C, size_t x, size_t y);
如果您想支持矩阵和 vector 的乘法,签名将变得更加复杂,并且完整的解决方案超出了本答案的范围。但在这种情况下我会编写一个单独的函数。
您可能想要研究的另一个选项是使用此签名。
typedef dm double[8][8]; // Short for dct matrix
dm *matrix_multiply(const dm *A, const *dm B, dm *output);
这里有一定的意义,因为您使用的是固定大小,但一般来说,这种方法并不常见。原因很简单,它的好处并没有真正超过必须处理指向二维数组的指针的缺点。
如果您想推迟重写但仍然能够按照您想要的方式初始化,请编写一个转换器。
double **convert (double mat[8][8])
{
double **ret = malloc (8 * sizeof *ret);
if (!ret) return NULL;
for (int i = 0; i < 8; i++)
ret[i] = &mat[i][0];
return ret;
}
然后你就可以了
double **A = convert(DCT_matrix);
double **B = convert(other_matrix);
// Use A and B
free(A);
free(B);
关于arrays - 快速初始化C数组(写入),我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/72798455/
如何更改循环中变量的名称?比如 number1 、 number2 、 number3 、 number4 ? var array = [2,4,6,8] func ap ( number1: Int
我想设置 View 的背景颜色并在一定延迟后将其更改为另一种颜色。这是我的尝试方式: print("setting color 1") self.view.backgroundColor = UICo
我在使用 express-session 时遇到问题。 session 数据不会在请求之间持续存在。 正如您在下面的代码中看到的那样,/join 路由设置了一些 session 属性,但是当 /sur
我试图从叶渲染器获得一个非常简单的结果,用于快速 Steam 的 for 循环。 我正在上传叶文件 HTML,因为它不接受此处格式正确的代码 - 下面的pizza.swift代码- import
你们中有人有什么好的链接可以与我分享吗?我正在寻找一个 FAST 程序员编辑器,它可以非常快速地打开包含超过 100, 000 行代码的文件?我目前正在使用记事本自动取款机,打开一个 29000 行长
我现在正在处理眼动追踪数据,因此拥有一个巨大的数据集(想想数百万行),因此希望有一种快速的方法来完成此任务。这是它的简化版本。 数据告诉您眼睛在每个时间点正在查看的位置以及我们正在查看的每个文件。 X
我是新手,想为计时器或其他设备选择提示音。 如何打开此列表,以选择其中一种声音? Alert sound list 最佳答案 您将无法在应用中使用系统声音。 但是,您可以包括自己的声音文件,并将其显示
我编写了以下代码来构建具有顺序字符串的数组。 它的工作方式与我预期的一样,但我希望它能更快地运行。有没有更有效的方法在PowerShell中产生我想要的结果? 我是PowerShell的新手,非常感谢
我有一个包含一些非唯一行的矩阵,例如: x 尝试 y <- rle(apply(x, 1, paste, collapse = " ")) # y$lengths is the vector con
我的函数“keyboardWillShown”有问题。所以我想要的是菜单打开时,菜单正好出现在键盘上方。它可以在Iphone 8 plus,8、7、6上完美运行。但是,当我在模拟器上运行Iphone
我正在尝试通过Swift 5中的HTTP get方法从API提取数据。它在启动时成功加载了数据,但是当我刷新页面时,它说“索引超出范围”,这是因为数据是不再会在我的日志中读取,因此索引中没有任何内容。
我想做什么: 从我的数据库中获取时间戳并将其转换为用户的时区。 我的代码: let tryItNow = "\(model.timestampName)" let format = D
给定字体名称和字体大小,如何查找字符串的宽度(CGFloat)? (目标是将UIView的宽度设置为足以容纳字符串的宽度。) 我有两个字符串:一个重复“1”,重复36次,另一个重复“M”,重复36次。
我正在尝试解析此JSON ["Items": ( { AccountBalance = 0; AlphabetType = 3; Description = "\U0631\U
我在UINavigationBar内放置了一个UILabel。 我想根据navigationBar的高度增加该标签的字体大小。当navigationBar很大时,我希望字体大小更大;当滚动并缩小nav
我想将用户输入限制为仅有效数字并使用以下内容: func textView(_ textView: UITextView, shouldChangeTextIn range: NSRange, rep
目前我有一个包含超过 100.000 张图像的数据库,它们大小不一或类似,但我想为我的公司制作以下内容: 我插入/上传一张图片,系统返回最有可能相同的图片。我不知道使用什么算法,但它需要快速。我可以预
在我的 swift 项目中,我有一个按钮,我想在标签上打印按下该按钮的时间。 如何解决这个问题? 最佳答案 添加到DHEERAJ的答案中,您只需在func press(sender: UIButton
我必须发表评论,尝试在解析中导入数组。然而,有一个问题。 当我尝试从 Parse 加载数组时,我的输出是 ("Blah","Blah","Blah")这是一个元组...而不是一个数组 TT... 如何
我的应用程序有一个名为 MyDevice 的类,我用它来与硬件通信。该硬件是可选的,实例变量也是可选的: var theDevice:MyDevice = nil 然后,在应用程序中,我必须初始化设备
我是一名优秀的程序员,十分优秀!