目前问题:有不轴对称并按约5180及接近2000位置对称的强点,显示奇怪的干涉
下图是电子密度早期分布
下图是场强@t500
禁掉傅里叶之后这个问题就没有了?
图片上看不到了
考虑是傅里叶自带的周期性导致的类似于干涉的效果?
在修订E0函数时,意外发现,else配对if可能需要一个分号?
至少我的if……if……可能不对
所以应该是if……elseif…………else?
面试官 | 写 if 时不带 else,你的代码会更好!_Java精选的博客-CSDN博客
暂时不清楚情况,先把所有的if……if……if串改成if……elseif……else串。
跑起来感觉速度没有提升,但是可以考虑等出结果看看结果有没有变化(带线性傅里叶部分)
只让光场占入射区域一半空间
运算结果没有区别……
发现一个bug:乘的过程怎么出现k1没有被kk取代的情况?
修正。再算一次
场强@t=500
神tm诡异
出现了六大分隔
应该是某个地方kk的赋值出现了问题
去掉傅里叶部分
CUDA学习笔记3:CUFFT(CUDA提供了封装好的CUFFT库)的使用例子 - 爱国呐 - 博客园
现在考虑是原程序cufftplan部分行列关系搞反了,把数颠倒一下试试
好像正确了!
这个长度和直径看起来比较正确
拉低z轴分辨率,提高尺度到总长150μm,整体在半径50μm的入射下基本符合常识
附录:数据传送回主机用277ms,写hdf文件183ms