- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我正在计算稠密非对称矩阵 A 的特征值。为此,我使用 xGEHRD和 xHSEQR Lapack 例程首先计算 A 的上 Hessenberg 形式,然后仅计算所获得矩阵的特征值。
这两个例程都需要参数 LWORK,并且都提供了一种计算其最佳值的机制。我相信这个参数与缓冲技术的内部阻塞有关,但我不知道它是如何确定的。
使用查询机制获取LWORK的最优值的工作流应该是这样的:
int LWORK = -1;
float* OPT_LWORK = (float*) malloc( sizeof(float));
sgehrd_ (..., OPT_LWORK ,&LWORK, ...) // query optimal value for sgehrd
LWORK = (int) OPT_WORK
float* WORK = (float*) malloc( (int) sizeof(float) * LWORK );
sgehrd_ (..., WORK ,&LWORK, ...) // calculate Hessenberg
int LWORK = -1;
shseqr_ (..., OPT_LWORK ,&LWORK, ...) // query optimal value for shseqr
LWORK = (int) OPT_WORK
float* WORK = // possibly realloc with the new LWORK value
shseqr_ (..., WORK ,&LWORK, ...) // finally obtain eigenvalues
我做了一些测试,始终获得 WORK 数组维度的相同最佳值。如果值相同,我可以大大简化我的代码(不需要重新分配,只需一次调用即可确定 LWORK 的值,减少错误检查...)。
我的问题是,对于相同的矩阵和相同的 ILO 和 IHI 值,我可以假设两个例程的值相等吗?
最佳答案
查看sgehrd.f , 例程 sgehrd()
的 WORK
的最佳大小似乎是 N*NB
,其中 NB
是
NB = MIN( NBMAX, ILAENV( 1, 'SGEHRD', ' ', N, ILO, IHI, -1 ) )
其中 NBMAX=64
。因此,最优的LWORK
取决于N
、ILO
和IHI
。
查看shseqr.f ,最佳长度的计算 LWORK
更复杂:例程 slaqr0()
被调用...但是文件中的文档指出:
If LWORK = -1, then SHSEQR does a workspace query. In this case, SHSEQR checks the input parameters and estimates the optimal workspace size for the given values of N, ILO and IHI. The estimate is returned in WORK(1). No error message related to LWORK is issued by XERBLA. Neither H nor Z are accessed.
对于 sgehrd()
和 shseqr()
,WORK
的最佳长度可能不同。这是一个例子:
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
extern void sgehrd_(int *n,int* ilo, int* ihi, float* a,int* lda,float* tau,float* work, int* lwork,int* info);
extern void shseqr_(char* job,char* compz,int *n,int* ilo, int* ihi,float* h,int* ldh,float* wr,float* wi,float* z,int* ldz,float* work, int* lwork,int* info);
int main()
{
int n=10;
int ilo=n;
int ihi=n;
float* a=malloc(sizeof(float)*n*n);
int lda=n;
float* tau=malloc(sizeof(float)*(n-1));
int info;
char job='S';
char compz='I';
float* h=malloc(sizeof(float)*n*n);
int ldh=n;
float* wr=malloc(sizeof(float)*(n));
float* wi=malloc(sizeof(float)*(n));
float* z=malloc(sizeof(float)*n*n);
int ldz=n;
int LWORK = -1;
float* OPT_LWORK =(float*) malloc( sizeof(float));
sgehrd_ (&n,&ilo, &ihi, a,&lda,tau,OPT_LWORK ,&LWORK,&info); // query optimal value for sgehrd
if(info!=0){printf("sgehrd lwork=-1 : failed\n");}
LWORK = (int) OPT_LWORK[0];
printf("sgehrd,length of optimal work : %d \n",LWORK);
float* WORK = (float*) malloc( (int) sizeof(float) * LWORK );
sgehrd_ (&n,&ilo, &ihi, a,&lda,tau,WORK ,&LWORK,&info);// calculate Hessenberg
if(info!=0){printf("sgehrd execution : failed\n");}
LWORK = -1;
shseqr_ (&job,&compz,&n,&ilo, &ihi, h,&ldh, wr, wi,z,&ldz, OPT_LWORK ,&LWORK, &info); // query optimal value for shseqr
if(info!=0){printf("shgeqr lwork=-1 : failed\n");}
LWORK = (int) OPT_LWORK[0];
printf("shseqr,length of optimal work : %d \n",LWORK);
WORK = realloc(WORK,(int) sizeof(float) * LWORK );// possibly realloc with the new LWORK value
shseqr_ (&job,&compz,&n,&ilo, &ihi, h,&ldh, wr, wi,z,&ldz, WORK ,&LWORK, &info); // finally obtain eigenvalues
if(info!=0){printf("shgeqr execution : failed\n");}
free(OPT_LWORK);
free(WORK);
free(a);
free(tau);
free(h);
free(wr);
free(wi);
free(z);
}
通过gcc main.c -o main -llapack -lblas
编译
我的输出是:
sgehrd,length of optimal work : 320
shseqr,length of optimal work : 10
关于c - xGEHRD 和 xHSEQR 例程中 WORK 数组的维数在特征值计算中是否应该相等?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/26957900/
我有一个 if 语句,如下所示 if (not(fullpath.lower().endswith(".pdf")) or not (fullpath.lower().endswith(tup
然而,在 PHP 中,可以: only appears if $foo is true. only appears if $foo is false. 在 Javascript 中,能否在一个脚
XML有很多好处。它既是机器可读的,也是人类可读的,它具有标准化的格式,并且用途广泛。 它也有一些缺点。它是冗长的,不是传输大量数据的非常有效的方法。 XML最有用的方面之一是模式语言。使用模式,您可
由于长期使用 SQL2000,我并没有真正深入了解公用表表达式。 我给出的答案here (#4025380)和 here (#4018793)违背了潮流,因为他们没有使用 CTE。 我很欣赏它们对于递
我有一个应用程序: void deleteObj(id){ MyObj obj = getObjById(id); if (obj == null) { throw n
我的代码如下。可能我以类似的方式多次使用它,即简单地说,我正在以这种方式管理 session 和事务: List users= null; try{ sess
在开发J2EE Web应用程序时,我通常会按以下方式组织我的包结构 com.jameselsey.. 控制器-控制器/操作转到此处 服务-事务服务类,由控制器调用 域-应用程序使用的我的域类/对象 D
这更多是出于好奇而不是任何重要问题,但我只是想知道 memmove 中的以下片段文档: Copying takes place as if an intermediate buffer were us
路径压缩涉及将根指定为路径上每个节点的新父节点——这可能会降低根的等级,并可能降低路径上所有节点的等级。有办法解决这个问题吗?有必要处理这个吗?或者,也许可以将等级视为树高的上限而不是确切的高度? 谢
我有两个类,A 和 B。A 是 B 的父类,我有一个函数接收指向 A 类型类的指针,检查它是否也是 B 类型,如果是将调用另一个函数,该函数接受一个指向类型 B 的类的指针。当函数调用另一个函数时,我
有没有办法让 valgrind 使用多个处理器? 我正在使用 valgrind 的 callgrind 进行一些瓶颈分析,并注意到我的应用程序中的资源使用行为与在 valgrind/callgrind
假设我们要使用 ReaderT [(a,b)]超过 Maybe monad,然后我们想在列表中进行查找。 现在,一个简单且不常见的方法是: 第一种可能性 find a = ReaderT (looku
我的代码似乎有问题。我需要说的是: if ( $('html').attr('lang').val() == 'fr-FR' ) { // do this } else { // do
根据this文章(2018 年 4 月)AKS 在可用性集中运行时能够跨故障域智能放置 Pod,但尚不考虑更新域。很快就会使用更新域将 Pod 放入 AKS 中吗? 最佳答案 当您设置集群时,它已经自
course | section | type comart2 : bsit201 : lec comart2 :
我正在开发自己的 SDK,而这又依赖于某些第 3 方 SDK。例如 - OkHttp。 我应该将 OkHttp 添加到我的 build.gradle 中,还是让我的 SDK 用户包含它?在这种情况下,
随着 Rust 越来越充实,我对它的兴趣开始激起。我喜欢它支持代数数据类型,尤其是那些匹配的事实,但是对其他功能习语有什么想法吗? 例如标准库中是否有标准过滤器/映射/归约函数的集合,更重要的是,您能
关闭。这个问题不符合Stack Overflow guidelines .它目前不接受答案。 这个问题似乎与 help center 中定义的范围内的编程无关。 . 关闭 9 年前。 Improve
我一直在研究 PHP 中的对象。我见过的所有示例甚至在它们自己的对象上都使用了对象构造函数。 PHP 会强制您这样做吗?如果是,为什么? 例如: firstname = $firstname;
...比关联数组? 关联数组会占用更多内存吗? $arr = array(1, 1, 1); $arr[10] = 1; $arr[] = 1; // <- index is 11; does the
我是一名优秀的程序员,十分优秀!