Skip to content

FIX: fillup the elastic params before returning the reading model function #89

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 2 commits into from
Sep 14, 2025
Merged

FIX: fillup the elastic params before returning the reading model function #89

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
06f93ec
FIX: fillup the elastic params before returning the reading model fun…
Dengda98 Sep 14, 2025
6e1be3a
fix typo
Dengda98 Sep 14, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions pygrt/C_extension/include/grt/common/model.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -77,12 +77,12 @@ GRT_MODEL1D * grt_init_mod1d(MYINT n);
GRT_MODEL1D * grt_copy_mod1d(const GRT_MODEL1D *mod1d1);

/**
* 根据不同的omega， 更新模型的 \f$ k_a^2, k_b^2 \f$ 参数
* 根据不同的 omega， 计算衰减系数，更新弹性模量
*
* @param[in,out] mod1d `MODEL1D` 结构体指针
* @param[in] omega 复数频率
*/
void grt_update_mod1d_omega(GRT_MODEL1D *mod1d, MYCOMPLEX omega);
void grt_attenuate_mod1d(GRT_MODEL1D *mod1d, MYCOMPLEX omega);

/**
* 根据相速度和层位，计算 iy 层的 (c/vp)^2, (c/vs)^2 以及归一化垂直波数
Expand Down
12 changes: 8 additions & 4 deletions pygrt/C_extension/src/common/model.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -145,16 +145,17 @@ GRT_MODEL1D * grt_copy_mod1d(const GRT_MODEL1D *mod1d1){
}


void grt_update_mod1d_omega(GRT_MODEL1D *mod1d, MYCOMPLEX omega){
void grt_attenuate_mod1d(GRT_MODEL1D *mod1d, MYCOMPLEX omega){
MYREAL Va0, Vb0;
// MYCOMPLEX ka0, kb0;
MYCOMPLEX atna, atnb;
for(MYINT i=0; i<mod1d->n; ++i){
Va0 = mod1d->Va[i];
Vb0 = mod1d->Vb[i];

atna = (mod1d->Qainv[i] > 0.0)? grt_attenuation_law(mod1d->Qainv[i], omega) : 1.0;
atnb = (mod1d->Qbinv[i] > 0.0)? grt_attenuation_law(mod1d->Qbinv[i], omega) : 1.0;
// 圆频率实部为负数表明不考虑模型的 Q 值属性
// 在读入模型后需要需要运行一次本函数以填充弹性模量，见 grt_read_mod1d_from_file 函数
atna = (creal(omega) >= 0.0 && mod1d->Qainv[i] > 0.0)? grt_attenuation_law(mod1d->Qainv[i], omega) : 1.0;
atnb = (creal(omega) >= 0.0 && mod1d->Qbinv[i] > 0.0)? grt_attenuation_law(mod1d->Qbinv[i], omega) : 1.0;

mod1d->atna[i] = atna;
mod1d->atnb[i] = atnb;
Expand Down Expand Up @@ -391,6 +392,9 @@ GRT_MODEL1D * grt_read_mod1d_from_file(const char *command, const char *modelpat
GRT_SAFE_FREE_PTR(modarr);
GRT_SAFE_FREE_PTR(line);

// 填充弹性模量
grt_attenuate_mod1d(mod1d, -1);

return mod1d;
}

Expand Down
2 changes: 1 addition & 1 deletion pygrt/C_extension/src/dynamic/grn.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -159,7 +159,7 @@ void grt_integ_grn_spec(
#else
local_mod1d = mod1d;
#endif
grt_update_mod1d_omega(local_mod1d, omega);
grt_attenuate_mod1d(local_mod1d, omega);

// 是否要输出积分过程文件
bool needfstats = (statsstr!=NULL && ((grt_findElement_MYINT(statsidxs, nstatsidxs, iw) >= 0) || (grt_findElement_MYINT(statsidxs, nstatsidxs, -1) >= 0)));
Expand Down
Loading