- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我是数值线性代数的新手,我刚刚开始使用 LAPACK 和 BLAS。
是否有可以在打包存储和完整存储之间复制/转换对称矩阵的例程?
我找到了 dtrttp
,我可以用它来将 double 全对称矩阵转换为压缩存储。但是,这些例程适用于三角矩阵,因此对应的 dtpttr
只填充完整矩阵的三角形。我怎样才能填满另一半?
最佳答案
显而易见的解决方案是通过“自制/DIY”代码对矩阵进行对称化,风险是重新发明轮子。在 dtpttr
之后编写对称矩阵所需的 for 循环非常容易。 .
for(i=0;i<n;i++){
for(j=i+1;j<n;j++){
a[i*n+j]=a[j*n+i];
}
}
dtpttr
持续 0.24 秒。
gcc main.c -o main -llapack -lblas -lm
编译它:
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>
void dtrttp_(char* UPLO,int* N,double* A,int* LDA,double* AP,int* INFO);
void dtpttr_(char* UPLO,int* N,double* AP,double* A,int* LDA,int* INFO);
void daxpy_(int* N,double* DA,double* DX,int* INCX,double* DY,int* INCY);
void dtpttf_(char* TRANSR,char* UPLO,int* N,double* AP,double* ARF,int* INFO);
int main(int argc, char **argv)
{
int n=10;
int info;
double *a=malloc(n*n*sizeof(double));
double *packed=malloc((n*(n+1))/2*sizeof(double));
int i,j;
for(i=0;i<n;i++){
for(j=0;j<n;j++){
a[i*n+j]=i+j;
}
}
printf("matrix before pack\n");
for(i=0;i<n;i++){
for(j=0;j<n;j++){
printf("%g ",a[i*n+j]);
}
printf("\n");
}
printf("\n");
//pack
dtrttp_("U",&n,a,&n,packed,&info);
//unpack
memset(a,0,n*n*sizeof(double));
dtpttr_("U",&n,packed,a,&n,&info);
for(i=0;i<n;i++){
for(j=i+1;j<n;j++){
a[i*n+j]=a[j*n+i];
}
}
printf("matrix after unpack\n");
for(i=0;i<n;i++){
for(j=0;j<n;j++){
printf("%g ",a[i*n+j]);
}
printf("\n");
}
free(a);
free(packed);
printf("timing...\n");
n=10000;
a=malloc(n*n*sizeof(double));
packed=malloc((n*(n+1))/2*sizeof(double));
for(i=0;i<n;i++){
for(j=0;j<n;j++){
a[i*n+j]=i+j;
}
}
//pack
dtrttp_("U",&n,a,&n,packed,&info);
//unpack
memset(a,0,n*n*sizeof(double));
clock_t t;
t = clock();
dtpttr_("U",&n,packed,a,&n,&info);
t = clock() - t;
printf ("dtpttr took %f seconds.\n",((double)t)/CLOCKS_PER_SEC);
t = clock();
for(i=0;i<n;i++){
for(j=i+1;j<n;j++){
a[i*n+j]=a[j*n+i];
}
}
t = clock() - t;
printf ("symmetrize took %f seconds.\n",((double)t)/CLOCKS_PER_SEC);
free(a);
free(packed);
return 0;
}
关于lapack - 在打包存储和完整存储之间转换对称矩阵?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/28445648/
具体详细介绍请看下文: 在使用文件进行交互数据的应用来说,使用FTP服务器是一个很好的选择。本文使用Apache Jakarta Commons Net(commons-net-3.3.jar)
我在日志文件中收到这些警告: WARN 2013-01-15 00:08:15,550 org.eclipse.jetty.http.HttpParser- HttpParser Full for
我在使用特定网页时遇到问题。当我按下链接时,我收到应用程序错误(不是 http 错误等,而是应用程序级别错误)。 但是我打开了开发人员工具和网络控制台,我看到没有请求发送到服务器。 所以我双击并选择查
我没有组装经验,但这是我一直在做的。如果在通过程序集中的指针传递参数和调用函数时缺少任何基本方面,我希望输入。 例如,我想知道是否应该还原ecx,edx,esi,edi,。我读到它们是通用寄存器,但我
我没有组装经验,但这是我一直在做的。如果在通过程序集中的指针传递参数和调用函数时缺少任何基本方面,我希望输入。 例如,我想知道是否应该还原ecx,edx,esi,edi,。我读到它们是通用寄存器,但我
我正在尝试创建完整 uiscrollview 的快照,所有内容大小,我已经搜索了很多,并且我在 SO 上找到了一些东西,如下所示: Getting a screenshot of a UIScroll
我想复制一个包含以下结构的Vector,对我来说重要的是在修改复制的 vector 时保持原始Vector完整: public class objet_poid_n { public int
给定一个示例字符串 s = '嗨,我的名字是 Humpty-Dumpty,来自“爱丽丝,爱丽丝镜中奇遇记”',我想将其分成以下 block : # To Do: something like {l =
已关闭。此问题旨在寻求有关书籍、工具、软件库等的建议。不符合Stack Overflow guidelines .它目前不接受答案。 我们不允许提问寻求书籍、工具、软件库等的推荐。您可以编辑问题,以
我正在尝试创建一个正则表达式来查找文本中的 Linux 文件路径,但是正则表达式对我来说非常陌生。我有下面的代码片段,它将识别下面文件结构的开头。 .*(/bin/|/home/).* 完成正则表达式
我正在寻找远程托管的 JPG 的尺寸、宽度和高度。我已经了解了如何通过下载完整图像来执行此操作。 但是,如果我可以通过仅下载足以获取此信息的方式来做到这一点,那将是理想的。 典型的图像大小为 200K
有没有办法让下面的代码: import traceback def log(message): print "%s: %s" %(traceback.extract_stack()[0:-1]
关闭。这个问题需要更多focused .它目前不接受答案。 想改进这个问题吗? 更新问题,使其只关注一个问题 editing this post . 关闭 3 年前。 Improve this qu
git show 显示了修订版中所做的所有更改。但是,它会打印出所有更改——而不仅仅是文件名。 git show --stat 只显示文件名,但它把它们截断了!有没有办法获得已更改文件名的完整列表?
Closed. This question does not meet Stack Overflow guidelines。它当前不接受答案。 想要改善这个问题吗?更新问题,以便将其作为on-topi
当我在模板中调用我的模型 get_absolute_url 方法时,我想要一个绝对/完整的 url。在我的入门模型中,我有以下内容: def get_absolute_url(self): r
我正在使用 jQuery 1.5.1 这是我的代码: $('.cellcontent').animate({ left: '-=190'}, { easing: alert('start
我正在使用下面的方法删除条形图并使用新数据更新条形图,但这样做时出现了一个小故障/完整的图表消失 1 秒,直到加载新数据。但是是否可以通过仅增加/减少柱形而不实际消失图表来实现相同的目的。 d3.se
基于 this question 中的讨论,任何人都可以提供代码或代码链接,显示 NumericLiteralX 模块的完整实现(例如 this one )?我对 NumericLiteralX 模块
我的目标是检索网站的 html,并将其转换为可读的String。我下面的代码可以工作,但我遇到了一个技术问题:当我尝试检索 http://time.gov/HTML5 的 html 时,我在 andr
我是一名优秀的程序员,十分优秀!