wangdingxi大哥，duns并行你调通了吗？

shyang5

我编译通过了，但他带的算例parallel/cylinder却做不下去，你算通了不？

wangdingxi

shyang5老兄，已经调通了。开始我调用了算例，后来发现并行部分源代码有一个小小的错误，不改正是不能用的。但是他的算例却凑巧能够使用（不改正源代码）。

shyang5

能否把错误贴出了？奇怪了，我计算的时候报出“non-physical temperature"的错误！大概计算50步后。

wangdingxi

in the following files:
parupdate.f
  real_precision  qdon(ninflw,nsmall2),qrcv(ninflw,nsmall2)
parprepar.f
    real_precision coeftmp(8,nsmall2),grdtmp(3*nbig)
      integer pimx,pjmx,pkmx,ptotmx,pzon
      integer inds(ibig*ibig),ind_on(nbig),ind_on_tmp(nbig),
     >        iptstmp(0:8,nsmall2),ibuff(2+10*nblk)
we should change nsmall2 to nsmall5,if you want to know the reason,you can
refer to the DUNS document and get the meaning of nsmall2 and nsmall5,then you can understand them,which cost me several days to find them!Good luck!
In addition,I wanna to know whether you use DUNS under Linux platform?

shyang5

dingxi兄，非常感谢你能把错误贴出来，我想对于我们duns使用者都是很好的帮助，顺便说一下，我是用linux platform。

wangdingxi

It is my pleasure!At the same time,I want to thank you for your help!
I am using DUNS under linux platform,too!

shyang5

我在调试并行时，发现如下错误：
1.makesystem_linuxpar中，当我们加上 －dp参数时，MPIPREC＝MPI_REAL应当改为MPI_DOUBLE_PRECISION.最好是在dunssource中加上。
2.在parupdate中，更新时二维和三维没有分开，二维时只用到4个点，后4个点不要。
3.在update中dq1没有定义。
好在他带的算例能算了，但有一个奇怪的现象，运行并行时，输出信息：
the number of doner cells:319
the number of recv cells:400
the number of doner cells:0
the number of recv cells:319
怎么可能呢，在actrual loop时是irecv＝400，怎么在后面的loop中变为0了。
加上 print*,'irecv=',irecv后显示：
the number of doner cells:319
the number of recv cells:400
the number of doner cells:400
the number of recv cells:319
莫名奇妙！
至今没有找到原因。

wangdingxi

shyang5老兄说得对，parupdate只用到前四个参数，再者MPI_REAL也是（你可看得仔细），但是update.f里面那个dq1你可以不管，因为预编译后的代码中没有。那段用到dq1的代码，与编译后被去掉了