1381 out.open(file, std::ios::in);
1382 if (!(out >> s >> n >> m >>
size1 >>
sig)) {
1383 std::cerr <<
"read G from file " << file <<
" error in file" 1387 if (n > nt_ || m != ntau_ ||
size1 != size1_)
1391 for (j = 0; j <= ntau_; j++) {
1393 for (l = 0; l < sg; l++) {
1394 if (!(out >> real >> imag)) {
1395 std::cerr <<
"read G from file " << file <<
" error at mat (" 1396 << j <<
")" << std::endl;
1399 matptr(j)[l] = std::complex<T>(real, imag);
1403 for (i = 0; i <= n; i++) {
1404 for (j = 0; j <= i; j++) {
1406 for (l = 0; l < sg; l++) {
1407 if (!(out >> real >> imag)) {
1408 std::cerr <<
"read G from file " << file
1409 <<
" error at ret (" << i <<
"," << j <<
")" 1413 retptr(i, j)[l] = std::complex<T>(real, imag);
1417 for (i = 0; i <= n; i++) {
1418 for (j = 0; j <= ntau_; j++) {
1420 for (l = 0; l < sg; l++) {
1421 if (!(out >> real >> imag)) {
1422 std::cerr <<
"read G from file " << file
1423 <<
" error at tv (" << i <<
"," << j <<
")" 1427 tvptr(i, j)[l] = std::complex<T>(real, imag);
1431 for (j = 0; j <= n; j++) {
1432 for (i = 0; i <= j; i++) {
1434 for (l = 0; l < sg; l++) {
1435 if (!(out >> real >> imag)) {
1436 std::cerr <<
"read G from file " << file
1437 <<
" error at les (" << i <<
"," << j <<
")" 1441 lesptr(i, j)[l] = std::complex<T>(real, imag);
void resize(int nt, int ntau, int size1)
Resizes herm_matrix object with respect to the number of time points nt, points on the Matsubara bra...