首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >在ARPACK中为数组ipntr赋值

在ARPACK中为数组ipntr赋值
EN

Stack Overflow用户
提问于 2011-09-02 12:01:23
回答 1查看 245关注 0票数 2

我正在尝试弄清楚如何在ARPACK中将值赋给名为ipntr的数组,该数组用于存储指针。我对fortran比较陌生,我不知道这个数组的值最初是如何设置的。我正在查看的代码是随ARPACK一起分发的非对称驱动程序示例dndrv1.f,我不理解的部分如下

代码语言:javascript
运行
复制
      program dndrv1 

      integer           maxn, maxnev, maxncv, ldv
      parameter         (maxn=256, maxnev=12, maxncv=30, ldv=maxn)
c
c     %--------------%
c     | Local Arrays |
c     %--------------%
c
      integer           iparam(11), ipntr(14)
      logical           select(maxncv)
      Double precision
     &                  ax(maxn), d(maxncv,3), resid(maxn), 
     &                  v(ldv,maxncv), workd(3*maxn), 
     &                  workev(3*maxncv), 
     &                  workl(3*maxncv*maxncv+6*maxncv)
c
c     %---------------%
c     | Local Scalars |
c     %---------------%
c
      character         bmat*1, which*2
      integer           ido, n, nx, nev, ncv, lworkl, info, j,
     &                  ierr, nconv, maxitr, ishfts, mode
      Double precision
     &                  tol, sigmar, sigmai
      logical           first, rvec
c
c     %------------%
c     | Parameters |
c     %------------%
c
      Double precision
     &                  zero
      parameter         (zero = 0.0D+0)
c
c     %-----------------------------%
c     | BLAS & LAPACK routines used |
c     %-----------------------------%
c
      Double precision
     &                  dlapy2, dnrm2
      external          dlapy2, dnrm2, daxpy
c
c     %--------------------%
c     | Intrinsic function |
c     %--------------------%
c
      intrinsic         abs

      write(*,*) ipntr

write语句的结果是:

1606679396 32767 1606696480 32767 1606918048 32767 0 0 0

显然不存在类似于ipntr(1) = 1606679396的值,那么这些值是如何分配的呢?

谢谢

EN

回答 1

Stack Overflow用户

发布于 2011-09-02 14:32:34

如果变量没有赋值,那么它们的值将由内存中恰好驻留的任意随机字节序列决定。您应该看到,重新运行该程序将产生不同的输出序列。

票数 3
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/7279096

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档