architecture course work

这个因为 specify 了学校目录所以都放那个机子上跑，所以需要使用 SILO=1 和 OMP_NUM_THREADS=6，加到 bashrc 里面就可

然后要在本地安装 visit

https://visit-sphinx-github-user-manual.readthedocs.io/en/latest/Intro/Installing_VisIt.html

这个还是比较麻烦的，下载相应版本的 targz 包，然后再下载 install script

https://wci.llnl.gov/simulation/computer-codes/visit/executables

添加权限 chmod +x visit-install，使用 visit-install version platform directory 放在和 targz 同目录下，给的范例是

visit-install 3.0.0 linux-x86_64 /usr/local/visit 不修改安装地点，看这个地方要 sudo 了

然后添加到 path，在 ~/.bashrc 最后一行加入 export PATH="/usr/local/visit/bin:$PATH"

然后直接输入 visit 就可以运行，出现了找不到 py 的问题，我是 ubuntu 20 所以安装 sudo apt install python-is-python3 就可以启动了

使用 visit 的话，左边 open file 打开 silo 的那个文件夹，记得 visit 要在 ACACGS 的目录下打开这样好找。然后下面的 plot, add - volume x_nodal or p_nodal.

先使用 AVX 搞 ddot waxpby sparsemv

补充一下，不能自动从 deployment 那里下载新增文件，所以直接在 tool - deployment - browse remote host 打开下载文件

决定使用 AVX，反正这个也快一些。先处理 ddot，256 位寄存 64 的 double 精度，一共可以存四个

1 2	__m256d localResult = _mm256_setzero_pd(); __m256d xi = _mm256_load_pd(x + i);

加个 d 是代表 double

草，记住 jb 的 deployment 的第二个是相对路经啊！

尝试 remote 开发？https://www.jetbrains.com/help/clion/remote-projects-support.html#remote-toolchain

file - build/execute- toolchain choose remote

算了不要弄了，一会儿还得用 remote host的呢

在 cppflags 里面加入 -fopenmp -mavx

这次一个 vector 包含 4 个 double，将 4 个加起来

1
2
3

__m256d x;
__m256d s = _mm256_hadd_pd(x,x);
return ((double*)&s)[0] + ((double*)&s)[2];

因为总是出现异常访问问题所以用远程 deployment 试试看

https://youtrack.jetbrains.com/issue/CPP-24324

在 ssh 的 keepalive 要取消，或者设置成 300？也不知道为什么这样就能连的上……

然后参考这个

https://www.jetbrains.com/help/clion/remote-projects-support.html?keymap=secondary_intellij_idea_classic#Makefile-toolchain

说句实在话真的超级麻烦啊啊啊啊！牙白！

在 toolchain 设置好 remote 之后，在 makefile 选中（都在 setting 里的 build 里面搞）,然后 makefile 指定本地的文件，修改好 deployment 的 mapping，会自动导航到 remote 的 executable 去？是的，自动的。有时候刚刚打开会连接不上远程服务器，重新设置一下 setting 里的 deployment 就可以了

使用 debug 可以找到问题了，如果不行就重启罢（无慈悲

如果看 log 的话在 help 里面 show log in file 那里。

vmovapd (%rax),%ymm0 这个时候 segmentation fualt，访问 x ptr 的时候……出现了问题，访问方式出现了错误？

回到 sse 进行访问，求和的话 https://stackoverflow.com/questions/6996764/fastest-way-to-do-horizontal-sse-vector-sum-or-other-reduction

使用 v = _mm_hadd_pd(v, v);

然后 store 低位的 64 位 double re; _mm_storel_pd(&re, sum);

这样加起来就可以了

Basics of SIMD Programming

看这里，其实 array 和 vector 是可以互相转换的，直接改写类型就好了，像这里的 simd-sample/simd_float.c at master · sonots/simd-sample · GitHub

还有 smaple code

int a[8] __attribute__((aligned(16))) = { 1, 2, 3, 4, 5, 6, 7, 8 };

__vector signed int *va = (__vector signed int *) a;
/* va[0] = { 1, 2, 3, 4}, va[1] = { 5, 6, 7, 8 } */

vb = vec_add(va[0], va[1]);

DEBUG !

在 flag 里面加入 -g 获得更多

p var 打印变量

horizontal add in vector

__m256d x1, x2;
// calculate 4 two-element horizontal sums:
// lower 64 bits contain x1[0] + x1[1]
// next 64 bits contain x2[0] + x2[1]
// next 64 bits contain x1[2] + x1[3]
// next 64 bits contain x2[2] + x2[3]
__m256d sum = _mm256_hadd_pd(x1, x2);
// extract upper 128 bits of result
__m128d sum_high = _mm256_extractf128_pd(sum1, 1);
// add upper 128 bits of sum to its lower 128 bits
__m128d result = _mm_add_pd(sum_high, _mm256_castpd256_pd128(sum));
// lower 64 bits of result contain the sum of x1[0], x1[1], x1[2], x1[3]
// upper 64 bits of result contain the sum of x2[0], x2[1], x2[2], x2[3]

_mm256_extractf128_pd 是提取低位，如果是 1 那就高位

_mm256_castpd256_pd128 是转换到 128，因为直接切掉高位所以没有延迟

fmadd 是 512 的指令，残念，请用 add 吧

不要使用 clion 的功能了！命令行不香吗！

openmp

atmoic 不是什么都可以用的

c - Invalid form of pragma atomic - Stack Overflow

只有这些可以

x++;
x--;
++x;
--x;
x binop= expr;
x = x binop expr;
x = expr binop x;

嘛反正就这样了，我摆！