check_aif_flx_format.hpp

Go to the documentation of this file.

#ifndef CHECK_AIF_FLX_FORMAT_HPP
#define CHECK_AIF_FLX_FORMAT_HPP
 
#include <iostream>
#include <fstream>
#include <sstream>
#include <string>
#include "globals.hpp"
 
/* This subroutine checks if the flx.aif file is in the old format
 * It assumes nnode>1
 */
inline bool is_old_flx_format(const std::string& filename){
    // Open file
    std::ifstream file(filename);
 
    // Check if file could be opened
    if(!file.is_open()){
        printf("\nError: %s not found", filename.c_str());
        exit_XGC("\nThe flx.aif file is required.\n");
    }
 
    // Read file line by line
    std::string line;
    int line_count = 0;
    int n_entries_on_line_5 = 0;
    int n_entries_on_line_6 = 0;
    while (std::getline(file, line)){
        line_count++;
 
        // To distinguish between the two formats, we need only count the number of
        // entries on Lines 5 and 6.
        if(line_count==5) {
            // Count how many entries are on this line
            std::istringstream iss(line);
            int entry;
            while (iss >> entry){
                n_entries_on_line_5++;
            }
        }
 
        if(line_count==6) {
            // Count how many entries are on this line
            std::istringstream iss(line);
            int entry;
            while (iss >> entry){
                n_entries_on_line_6++;
            }
 
            // Exit loop after Line 6
            break;
        }
    }
    file.close();
 
    // In the new format, the only way that both Line 5 and Line 6 could be size 1 is if nnode==1, which we assume is not the case.
    // Also, if there were 0 entries due to the file being too short, that could only be the old format. Therefore we can distinguish them:
    if(n_entries_on_line_5<=1 && n_entries_on_line_6<=1){
        // Old format
        printf("  Flux surface file %s identified as using the OLD format. Consider converting to the new format.\n", filename.c_str());
        return true;
    }else{
        // New format
        printf("  Flux surface file %s identified as using the NEW format.\n", filename.c_str());
        return false;
    }
}
 
#endif