read_from_file() [2/2]

template<typename T >

void cntr::function< T >::read_from_file	(	int	nt1,
		const char *	file
	)

Read the function data from a txt file

Purpose

Read the function data from a txt file. The file is probably written by print_to_file method, or is assumed to have the same format.

Parameters

nt1	the number of time steps to be read
file	the file name

Definition at line 790 of file cntr_function_impl.hpp.

                                                          {
    int i, n, l, size1, sg;
    double real, imag;
    std::string s;
    std::ifstream out;
    out.open(file, std::ios::in);
    if (!(out >> s >> n >> size1)) {
        std::cerr << "read function from file " << file << " error in file"
                  << std::endl;
        abort();
    }

    assert(nt1 >= -1 && nt1 <= nt_);
    assert(nt1 <= n);
    assert(size1 == size1_);
    
    sg = size1 * size1;
    for (i = -1; i <= n; i++) {
        for (l = 0; l < sg; l++) {
            if (!(out >> real >> imag)) {
                std::cerr << "read function from file " << file
                          << " error at  t= " << i << std::endl;
                abort();
            }
            if (i <= nt1)
                ptr(i)[l] = std::complex<T>(real, imag);
        }
    }
    out.close();
}