etranger的博客

0.不带参数与rawlinson版本完全相同,不同点在于不需要另外计算ak135走时表，内置iasp91模型并根据以前的算法实时计算远震。 1.带一个正整数参数或负整数参数，表示并行多线程的个数，正数表示存储全部波场按检波点输出，负数不存储波场按震源输出 注意原FMTOMO无法在no_pp_mode下计算Pn震相，这个bug在这个参数下已经修复，输出arrivals.dat和frechet.dat。 2.带时间走时参数和阻尼参数，表示根据当前走时otimes.dat数据对比arrivals.dat和frechet.dat导数数据进行一次反演。阻尼参数为正表示进行 subspace子空间阻尼反演，而负阻尼表示进行全稀疏矩阵的阻尼最小二乘反演。 3.带时间走时参数和阻尼参数和多线程的个数，表示根据当前模型(包括界面)进行多线程计算走时和偏导数(不输出)，多线程参数同1， 然后进行迭代反演，阻尼参数如2， 直至拟合otimes.dat残差变化小于10ms. 4.根据参数1读取走时，生成震源文件sources.in和接收点文件receivers.in以及走时文件otimes.dat. 根据参数2，3，4生成界面interfaces.in， 根据crust2.0模型 和iasp91模型生成速度模型vgrids.in。 2，3，4分别为地表,Moho和底界面的grd文件 5.根据lat1 lon1 lat2 lon2 dep对模型进行纵切片(前4个参数)和深度切片(dep)，并生成grd文件和ps文件。 6.根据存在的模型文件，结合给定的参数进行分辨率测试反演。参数1是速度分辨率，参数2是界面分辨率，参数3是震源定位分辨率，参数4给予 走时扰动，5，6参数分别位阻尼和线程数参数。

在命令行的基础上增加了界面，左键震源，右键接收点，中间滚动放大缩小，按住中键拖动移动位置 状态条显示坐标和震中距，射线上显示走时

Usage:iasp<phase/Psdep> <delta> [depth=0] [initrayp] [recdep] [rayperr] [deltaerr] [Nsamp]

if(atof(argv[1])>0){code[0]=2;code[1]=-1;z3[0]=R-atof(argv[1]);nm=2;} ("SKS")) ("SKiKS") ("SKIKS") ("PKiKP") ("PKIKP") ("P")) ("p")) ("S")) ("s")) ("Pm") ("Sm") ("Pms" ("Smp" ("Ps") ("Sp") ("PP") ("SS") ("PS") ("SP") ("PcP" ("ScS" ("PcS" ("ScP" ("Pn") ("Sn") ("Pdiff")) ("Sdiff"))

采用mpreal任意精度数学库，可以计算pi到2亿位

算法数独C++界面加高级建议

大地坐标经纬度互相转换，有北京坐标，西安坐标和3度分带和6度分带

c 匹配认字的字母c .(句点) 匹配任意的单个字符 ^ 匹配输入字符串的开头 $ 匹配输入字符串的结尾 * 匹配前一个字符的零个或者多个出现

能够解较难的题，用法:sudoku_console <谜题文本文件> 谜题例子 0 7 0 0 0 0 0 0 8 0 0 0 6 0 3 0 7 0 0 0 6 0 0 1 0 9 2 0 4 1 0 0 7 9 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 4 8 2 1 7 0 0 1 0 0 5 0 0 5 0 0 0 0 0 0 0 3 0 3 0 0 7 0 0 0 0 结果在ans.txt中，屏幕上显示解题过程