首页 文章
  • 1 votes
     answers
     views

    具有16字节整数索引的ARPACK特征值

    我的代码可以很好地计算我的测试用例中的特征值为ARPACK无耻地取自here并适应快速的4x4矩阵 . (为简洁起见,我在示例代码中删除了顶部的注释) . 好的,我的问题 . 我有非常大的矩阵,或者至少,我会为我的实际问题 . 但是,当我使整数类型为16时,ARPACK会出错 . 有没有一种简单的方法来转换ARPACK函数以允许我的16字节索引的东西?或者,是否可以改变它使库允许的方式?我用gf...
  • 0 votes
     answers
     views

    如何计算Fortran中的第一个特征值和特征向量

    我尝试使用ARPACK使用一些例子(here),但我甚至无法弄清楚如何输入我的矩阵 . 从这个question,似乎Python和Matlab中的实现是避免ARPACK复杂性的唯一解决方案 . 是否有一些(基于Fortran)的方法来做到这一点,还避免计算所有的特征值/特征向量?
  • 2 votes
     answers
     views

    ARPACK页面排名错误

    我在R中使用igraph包来运行PageRank . 我为我的图运行它并得到了这个错误: 在arpack.c:1130:ARPACK错误,由LAPACK例程dlahqr计算的Schur形式无法由LAPACK例程dtrsen重新排序 . 我在网上看到,这可能是由于你的图表有多个不相交的组件造成的 . 所以我跑了 gList = decompose.graph(g) is.connected(gL...
  • 0 votes
     answers
     views

    kknn中的specClust() - arpack迭代限制增加

    我正在将谱聚类应用于具有4200行和2列的数据集 . spec <- specClust(df1, centers=7, nn = 7, method = "symmetric") 我有以下错误 . n .Call("R_igraph_arpack", func, extra, options, env, sym, PACKAGE = "i...
  • 5 votes
     answers
     views

    ARPACK线程安全吗?

    从用C编写的程序中同时使用不同线程的ARPACK本征解析器是否安全?或者,如果ARPACK本身不是线程安全的,那么是否存在与API兼容的线程安全实现?一个快速的谷歌搜索没有发现它是第一个需要一个线程安全的稀疏eigensolver非常令人惊讶 . 我对Fortran不太熟悉,所以我使用 f2c 将ARPACK源代码翻译成C,似乎有很多静态变量 . 基本上,转换例程中的所有局部变量似乎都是静态的,这...
  • 2 votes
     answers
     views

    在ARPACK中将值分配给数组ipntr

    我试图弄清楚如何将值分配给名为ipntr的数组,该数组存储指针,在ARPACK中 . 我对fortran相对较新,我无法弄清楚这个数组的值是如何初始设置的 . 我正在看的代码是一个不对称的驱动程序示例dndrv1.f与ARPACK一起分发,我不理解的部分如下 program dndrv1 integer maxn, maxnev, maxncv, ldv ...
  • 1 votes
     answers
     views

    链接到VS2010的ARPACK

    我想让 ARPACK 库在 VS2010 上运行 . 我想使用 ARPACK 提供的C包装器(一些背景 - 我需要得到巨大矩阵的特征值和特征向量) . 有关该主题的非常好的教程here . 在本教程之后,我已经设法使用 g77 在 mingw 上编译了fortran代码,我成功生成了所描述的dll和lib . 尝试将我的visual studio项目链接到库时会出现问题 . 我试图链接的方式如下:...
  • 3 votes
     answers
     views

    Scipy arpack eigs与eigsh特征值的数量

    在scipy的ARPACK绑定中,无法计算矩阵的所有特征值 . 但是,我发现eigsh能够计算n - 1个特征值,而eigs只能计算n - 2个特征值 . 任何人都可以验证这实际上是ARPACK的基本限制而不是scipy中的错误吗? 这是示例代码: import scipy.sparse, scipy.sparse.linalg t = scipy.sparse.eye(3,3).tocs...
  • 2 votes
     answers
     views

    Windows上的Arpack()

    我试图让Arpack和Arpack在Windows 7中使用Visual Studio 2010,并使用swig创建一个可以与C#一起使用的界面 . 这基本上是为了创建一个基准,我将用它来衡量用C#解决特征向量问题的所有进一步尝试,但首先我需要这个工作 . 我已经设法使用this guide编译Arpack,虽然稍有改动,因为g77不再是MinGW附带的Fortran编译器,而是我们得到了gfor...
  • 1 votes
     answers
     views

    使用ARPACK查找稀疏矩阵的特征向量和特征值(称为PYTHON,MATLAB或FORTRAN子例程)

    几天前我问了一个如何找到大稀疏矩阵的特征值的问题 . 我没有答案,所以我决定描述一个潜在的解决方案 . One question remains: Can I use the python implementation of ARPACK to compute the eigenvalues of a asymmetric sparse matrix. 首先,我想说,根本不需要使用FOTRA...
  • 0 votes
     answers
     views

    用Arpack计算特征值

    我正在使用Arpack来计算问题的特征值\ lambda Mx = Ax,采用具有复位移的反向移位方法 . 然后,我使用znaupd e zneupd . 一切似乎都没问题 . 我写了一个函数,它采用Arpack计算的特征值和特征向量,并将它们放在原始问题中:\ lambda Mx -Ax =残差 . 残差的范数应该几乎为零或非常小 . 无论如何,我得到了一个非常高的10 ^ -3的错误!你知道问...
  • 0 votes
     answers
     views

    ARPACK对于Hermitian矩阵的奇怪行为

    我想在数值上获得一些厄米特矩阵的基态能量(参见下面代码中该矩阵的定义),并根据矩阵 - 参数“相位”绘制它 . import scipy.sparse as sparse import scipy import numpy import numpy as np import math from scipy.special import binom import cmath import sympy...
  • 0 votes
     answers
     views

    试图构建ARPACK,但许多源文件没有被编译

    首先我要说的是,我不是超级用户,有点像菜鸟,但我对我正在做的事情有一个基本的了解 . 我按照here的指示安装了MSYS2和MinGW-w64 . 我正在使用64位Windows 10,我想最终在我的Fortran 90代码中使用ARPACK库 . 我目前正在使用Microsoft Visual Studio 2017和Intel Fortran ifort编译器和MKL库 . 所以, 1. 我应...

热门问题