input_output.h

Go to the documentation of this file.

/*
 * This file is part of Quantum++.
 *
 * MIT License
 *
 * Copyright (c) 2013 - 2020 Vlad Gheorghiu.
 *
 * Permission is hereby granted, free of charge, to any person obtaining a copy
 * of this software and associated documentation files (the "Software"), to deal
 * in the Software without restriction, including without limitation the rights
 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
 * copies of the Software, and to permit persons to whom the Software is
 * furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in
 * all copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
 * SOFTWARE.
 */

#ifndef INPUT_OUTPUT_H_
#define INPUT_OUTPUT_H_

namespace qpp {
template <typename Derived>
internal::IOManipEigen disp(const Eigen::MatrixBase<Derived>& A,
                            double chop = qpp::chop) {
    return internal::IOManipEigen(A, chop);
}

inline internal::IOManipEigen disp(cplx z, double chop = qpp::chop) {
    return internal::IOManipEigen(z, chop);
}

template <typename InputIterator>
internal::IOManipRange<InputIterator>
disp(InputIterator first, InputIterator last, const std::string& separator,
     const std::string& start = "[", const std::string& end = "]",
     double chop = qpp::chop) {
    return internal::IOManipRange<InputIterator>(first, last, separator, start,
                                                 end, chop);
}

template <typename Container>
internal::IOManipRange<typename Container::const_iterator>
disp(const Container& c, const std::string& separator,
     const std::string& start = "[", const std::string& end = "]",
     double chop = qpp::chop,
     typename std::enable_if<is_iterable<Container>::value>::type* = nullptr) {
    return internal::IOManipRange<typename Container::const_iterator>(
        std::begin(c), std::end(c), separator, start, end, chop);
}

template <typename PointerType>
internal::IOManipPointer<PointerType>
disp(const PointerType* p, idx N, const std::string& separator,
     const std::string& start = "[", const std::string& end = "]",
     double chop = qpp::chop) {
    return internal::IOManipPointer<PointerType>(p, N, separator, start, end,
                                                 chop);
}

template <typename Derived>
void save(const Eigen::MatrixBase<Derived>& A, const std::string& fname) {
    const dyn_mat<typename Derived::Scalar>& rA = A.derived();

    // EXCEPTION CHECKS

    // check zero-size
    if (!internal::check_nonzero_size(rA))
        throw exception::ZeroSize("qpp::save()");

    std::fstream fout;
    fout.open(fname, std::ios::out | std::ios::binary);

    if (fout.fail()) {
        throw std::runtime_error("qpp::save(): Error writing output file \"" +
                                 std::string(fname) + "\"!");
    }
    // END EXCEPTION CHECKS

    // write the header to file
    const std::string header_ = "TYPE::Eigen::Matrix";
    fout.write(header_.c_str(), header_.length());

    idx rows = static_cast<idx>(rA.rows());
    idx cols = static_cast<idx>(rA.cols());
    fout.write(reinterpret_cast<const char*>(&rows), sizeof(rows));
    fout.write(reinterpret_cast<const char*>(&cols), sizeof(cols));

    fout.write(reinterpret_cast<const char*>(rA.data()),
               sizeof(typename Derived::Scalar) * rows * cols);

    fout.close();
}

template <typename Derived>
dyn_mat<typename Derived::Scalar> load(const std::string& fname) {
    std::fstream fin;
    fin.open(fname, std::ios::in | std::ios::binary);

    // EXCEPTION CHECKS

    if (fin.fail()) {
        throw std::runtime_error("qpp::load(): Error opening input file \"" +
                                 std::string(fname) + "\"!");
    }

    const std::string header_ = "TYPE::Eigen::Matrix";
    std::unique_ptr<char[]> fheader_{new char[header_.length()]};

    // read the header from file
    fin.read(fheader_.get(), header_.length());
    if (std::string(fheader_.get(), header_.length()) != header_) {
        throw std::runtime_error("qpp::load(): Input file \"" +
                                 std::string(fname) + "\" is corrupted!");
    }
    // END EXCEPTION CHECKS

    idx rows, cols;
    fin.read(reinterpret_cast<char*>(&rows), sizeof(rows));
    fin.read(reinterpret_cast<char*>(&cols), sizeof(cols));

    dyn_mat<typename Derived::Scalar> A(rows, cols);

    fin.read(reinterpret_cast<char*>(A.data()),
             sizeof(typename Derived::Scalar) * rows * cols);

    fin.close();

    return A;
}

} /* namespace qpp */

#endif /* INPUT_OUTPUT_H_ */