Science/cloudy-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/cloudy-agn.git synced 2024-12-05 02:25:08 +00:00
Code Issues Packages Projects Releases Wiki Activity

added bulk interpolation fix script

Browse Source
This commit is contained in:
caes 2016-06-10 01:05:29 -04:00
parent 83dcad1348
commit 84c621f35e
3 changed files with 43 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
scripts/meta/bulk_interpolation_fix.sh Executable file
View File

@ -0,0 +1,25 @@
#!/bin/bash
script_dir=$( cd $(dirname $0) ; pwd -P |sed 's@^\(.*\)/scripts.*@\1/scripts@')
basedir=`pwd`
if [[ -x $script_dir/../bin/interpolate_fix]];
then
echo "Binary available. Proceeding."
else
echo "Binary unavailable. Exiting."
exit 9
fi
#for fortfile in $(echo fort.*);
for fortfile in $@
do
if [[ -e $(basename $fortfile) ]]
then
mkdir -p uninterpolated
mv $fortfile uninterpolated
echo $script_dir/../bin/interpolate_fix uninterpolated/$(basename $fortfile) $fortfile
else
echo $script_dir/../bin/interpolate_fix $fortfile $(basename $fortfile)
fi
exit 0

28
src/interpolation_fix.cpp
View File

@ -34,7 +34,7 @@ int main(int argc, char const *argv[]) {
<< outlier_it->second
<< "\n"
<< std::scientific
<< "Old value: "
<< "Old value: "
<< jagged_table.value[outlier_it->first][outlier_it->second]
<< "\t"
<< "New value: "
@ -43,7 +43,7 @@ int main(int argc, char const *argv[]) {
outlier_it++;
}
std::ofstream smoothed_table_file;
smoothed_table_file.open("smoothed_eqwidth_table");
smoothed_table_file.open(argv[2]);
smoothed_table_file << agn::format_eqwidth_table(smooth_table);
smoothed_table_file.close();
@ -118,13 +118,13 @@ agn::eqwidth_table agn::smooth(agn::eqwidth_table & jagged_table,agn::gridcoordl
if ( agn::debug ) {
std::cout
<< std::scientific
<< "Smoothed: "
<< "Smoothed: "
<< newx
<< ", "
<< oldy
<< oldy
<< " -> "
<< newy
<< "\n";
<< "\n";
}
coord_it++;
@ -150,23 +150,23 @@ agn::gridcoordlist agn::known_outliers() {
}
agn::gridcoordlist agn::find_outliers(agn::eqwidth_table jagged_table) {
if (agn::debug) std::cout
if (agn::debug) std::cout
<< "Scanning for outliers.\n";
agn::gridcoordlist outliers;
// The distribution is
// parameterised by x=log(jagged_table.hden) and
// The distribution is
// parameterised by x=log(jagged_table.hden) and
// y=log(jagged_table.eqwidth).
agn::iterator2d hden_it = jagged_table.value.begin();
while ( hden_it != jagged_table.value.end() ) {
agn::table1d slice = hden_it->second;
if( agn::debug ) std::cout
if( agn::debug ) std::cout
<< "\nhden= "
<< (hden_it->first)
<< (hden_it->first)
<< ": ";
// Crawl, checking slope, until a positive slope is found.
// Crawl, checking slope, until a positive slope is found.
agn::table1d::iterator phi_it = slice.begin();
double x1,x2,y1,y2,slope;
std::vector<double> ref_curve_x,ref_curve_y,curve_back_x,curve_back_y;
@ -220,7 +220,7 @@ agn::gridcoordlist agn::find_outliers(agn::eqwidth_table jagged_table) {
phi_it++;
}
if ( phi_it == slice.end() ) {
if ( phi_it == slice.end() ) {
if (agn::debug) std::cout
<< "Clean";
hden_it++;
@ -282,7 +282,7 @@ agn::gridcoordlist agn::find_outliers(agn::eqwidth_table jagged_table) {
std::vector<double>::iterator curve_back_x_it = curve_back_x.begin();
std::vector<double>::iterator curve_back_y_it = curve_back_y.begin();
while ( curve_back_x_it != curve_back_x.end() &&
while ( curve_back_x_it != curve_back_x.end() &&
curve_back_y_it != curve_back_y.end() ) {
ref_curve_x.push_back(*curve_back_x_it);
ref_curve_y.push_back(*curve_back_y_it);
@ -329,7 +329,7 @@ agn::gridcoordlist agn::find_outliers(agn::eqwidth_table jagged_table) {
// << ", ";
phi_it = slice.find(anomalies_start_x);
while (phi_it != (++(slice.find(anomalies_end_x)))) {
y1 = log10(phi_it->second);
y1 = log10(phi_it->second);
x1 = phi_it->first;
y2 = testing_spline[x1];
//threshold = RATIO_THRESHOLD_MULTIPLIER * (y2 - EQWIDTH_MIN_VAL_LOG);

8
src/spectral_lines.hpp
View File

@ -143,7 +143,7 @@ std::string agn::format_eqwidth_table_slice(agn::eqwidth_table table,agn::iterat
std::ostream& agn::operator<< (std::ostream& outstream, agn::eqwidth_table table) {
outstream
outstream
<< table.header[0]
<< std::endl
<< table.header[1]
@ -211,9 +211,9 @@ std::list<agn::eqwidth_table> agn::compile_eqwidth_tables(agn::cloudy_grid grid,
double x,y;
agn::cloudy_line_data data;
agn::eqwidth_table new_table;
if(agn::debug) std::cout
<< "Processing label "
<< *line_label_it
if(agn::debug) std::cout
<< "Processing label "
<< *line_label_it
<< std::endl;
new_table.header[0] = *line_label_it;
new_table.header[0].append(" relative to Inci 1215.00A scaled to 1215.00A");