pardiso.parser.hpp

Go to the documentation of this file.

/*******************************************************************************
 * Copyright (C) 2025-2026 Theodore Chang
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 ******************************************************************************/
#ifndef PARDISO_PARSER_HPP
#define PARDISO_PARSER_HPP
 
#ifdef EZP_MKL
 
#include "pardiso.hpp"
#include "parser.hpp"
 
namespace ezp {
    template<data_t DT, index_t IT> auto pardiso_set(const std::string& command, pardiso<DT, IT>& solver) {
        for(const auto& [token, value] : detail::parse(command)) {
            if(token == "--default-value") solver.iparm_default_value(value);
            else if(token == "--reducing-ordering") solver.iparm_reducing_ordering(value);
            else if(token == "--user-permutation") solver.iparm_user_permutation(value);
            else if(token == "--iterative-refinement") solver.iparm_iterative_refinement(value);
            else if(token == "--pivoting-perturbation") solver.iparm_pivoting_perturbation(value);
            else if(token == "--scaling") solver.iparm_scaling(value);
            else if(token == "--transpose-matrix") solver.iparm_transpose_matrix(value);
            else if(token == "--weighted-matching") solver.iparm_weighted_matching(value);
            else if(token == "--nnz-factor") solver.iparm_nnz_factor(value);
            else if(token == "--pivoting-type") solver.iparm_pivoting_type(value);
            else if(token == "--matrix-checker") solver.iparm_matrix_checker(value);
            else if(token == "--partial-solve") solver.iparm_partial_solve(value);
            else if(token == "--zero-based-indexing") solver.iparm_zero_based_indexing(value);
            else if(token == "--schur-complement") solver.iparm_schur_complement(value);
            else if(token == "--out-of-core") solver.iparm_out_of_core(value);
        }
    }
} // namespace ezp
 
#endif
 
#endif