[PIPE2D-460] Evaluate the interpolations between model templates Created: 27/Sep/19  Updated: 12/Nov/19  Resolved: 08/Nov/19

Status: Done
Project: DRP 2-D Pipeline
Component/s: None
Affects Version/s: None
Fix Version/s: None

Type: Task Priority: Normal
Reporter: Takuji Yamashita Assignee: Takuji Yamashita
Resolution: Done Votes: 0
Labels: flux-calibration, model-templates
Remaining Estimate: Not Specified
Time Spent: Not Specified
Original Estimate: Not Specified

Attachments: PNG File numspec_time.png     PNG File pmodel_rms_4param.png     PNG File pmodel_rms_alpha.png     PNG File pmodel_rms_logg.png     PNG File pmodel_rms_numspec.png     PNG File pmodel_rms_Teff.png     PNG File pmodel_rms_Z.png     PNG File rms_table.png     PNG File smodel_rms_4param.png     PNG File smodel_rms_alpha.png     PNG File smodel_rms_logg.png     PNG File smodel_rms_numspec.png     PNG File smodel_rms_Teff.png     PNG File smodel_rms_Z.png     PNG File time_min.png    
Issue Links:
Blocks
blocks PIPE2D-365 Interpolate between model templates Done
Story Points: 4
Sprint: 2DDRP-2019 H, 2DDRP-2019 I
Reviewers: hassan

 Description   

Evaluate the accuracy of interpolations of the AMBRE model spectra. The interpolations were provided in PIPE2D-365, where 4D radial basis functions were used for interpolations.



 Comments   
Comment by Takuji Yamashita [ 25/Oct/19 ]

To evaluate the accuracies, the interpolations were made at each parameter point that we have. 
An RBF fitting was performed for every spectrum. A document is in preparation. Target value of the accuracy is <1%.
 
Model templates
All the AMBRE templates with 6000 K < Teff < 8000 K

  • 2590 models of the plan-parallel model (p-model)
  • 2125 models of the spherical model (s-model)

Computer
3.8GHz CPU, 62GB memory, CentOS7, 7 cores are used.
 
Processing time
~10 days for 2590 p-model spectra; ~5 days for 2125 s-model spectral. On average, 5 min for p-model, 3 min for s-model. Most of the interpolations takes only ~10 sec. Interpolations which use many neighbor spectra take much time. 

 
RMS of residuals
For each interpolation spectrum, a residual is calculated as a function of wavelength.
Residual (%) == (Interpolation - Reference) / Reference * 100, 
where Reference is an original AMBRE spectrum.
Overall, the accuracies is good enough. On average, they are 0.03 (/-0.03)% and 0.02 (/-0.04)% in p-mode and s-model, respectively. I found that accuracies become worse a bit at lower Teff and high metallicity in both p/s-models.
 
P-model






 
S-model






 

P and S-models
Both models reproduce with high accuracies. The results show that the accuracy of p-model is systematically higher than that of s-model. This is because the parameter grid of the p-model templates is dense like a regular grid while one of s-model is sparse. 

Comment by Takuji Yamashita [ 08/Nov/19 ]

The results of the evaluations was provided and a code for this work was uploaded at github. 

Generated at Sat Feb 10 15:53:36 JST 2024 using Jira 8.3.4#803005-sha1:1f96e09b3c60279a408a2ae47be3c745f571388b.