iostats.h

对积分过程中的核函数和积分值进行记录

Author

Zhu Dengda (zhudengda@mail.iggcas.ac.cn)

Date

2024-07-24

Functions

void grt_write_stats(FILE *f0, real_t k, const cplx_t QWV[GRT_SRC_M_NUM][GRT_QWV_NUM])

将积分过程中计算的核函数写入文件

备注

文件记录的值均为波数积分的中间结果,与最终的结果还差一系列的系数, 记录其值主要用于参考其变化趋势。

参数:
  • f0[out] 二进制文件指针

  • k[in] 波数

  • QWV[in] 不同震源,不同阶数的核函数 \( q_m, w_m, v_m \)

int grt_extract_stats(FILE *bf0, FILE *af0)

从二进制核函数文件读出一个数据块,写入到文本文件中

参数:
  • bf0[inout] 二进制文件指针,如果为NULL则打印标题

  • af0[out] 文本文件指针

返回:

0表示读取成功,-1表示读取结果/失败

void grt_write_stats_ptam(FILE *f0, real_t Kpt[GRT_SRC_M_NUM][GRT_INTEG_NUM][GRT_PTAM_PT_MAX], cplx_t Fpt[GRT_SRC_M_NUM][GRT_INTEG_NUM][GRT_PTAM_PT_MAX])

记录峰谷平均法的峰谷位置

备注

文件记录的积分值与最终的结果还差一系列的系数, 记录其值主要用于参考其变化趋势。

参数:
  • f0[out] 二进制文件指针

  • Kpt[in] 最终收敛积分值使用的波峰波谷位置

  • Fpt[in] 最终收敛积分值使用的波峰波谷幅值

int grt_extract_stats_ptam(FILE *bf0, FILE *af0)

从二进制峰谷位置文件读出一个数据块,写入到文本文件中

参数:
  • bf0[inout] 二进制文件指针,如果为NULL则打印标题

  • af0[out] 文本文件指针

返回:

0表示读取成功,-1表示读取结果/失败