1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
// BEGIN LICENSE BLOCK
/*
Copyright (c) 2009 , UT-Battelle, LLC
All rights reserved
[PsimagLite, Version 1.0.0]
*********************************************************
THE SOFTWARE IS SUPPLIED BY THE COPYRIGHT HOLDERS AND
CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED
WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
PARTICULAR PURPOSE ARE DISCLAIMED.
Please see full open source license included in file LICENSE.
*********************************************************
*/
// END LICENSE BLOCK
#include "ContinuedFraction.h"
#include "ContinuedFractionCollection.h"
#include "Io/IoSimple.h"
#include "TridiagonalMatrix.h"
#include <cstdlib>
#include <fstream>
#include <iostream>
#include <string>
using namespace PsimagLite;
typedef double RealType;
typedef TridiagonalMatrix<RealType> TridiagonalMatrixType;
typedef ContinuedFraction<TridiagonalMatrixType> ContinuedFractionType;
typedef ContinuedFractionCollection<ContinuedFractionType>
ContinuedFractionCollectionType;
void usage(const char* progName)
{
std::cerr << "Usage: " << progName << " file1 file2\n";
}
int main(int argc, char* argv[])
{
if (argc < 2) {
usage(argv[0]);
return 1;
}
ContinuedFractionCollectionType cfCollection(FREQ_REAL);
String s = "#Avector";
for (int x = 1; x < argc; x++) {
IoSimple::In io(argv[x]);
io.advance(s, IoSimple::In::LAST_INSTANCE);
ContinuedFractionType cf(io);
cfCollection.push(cf);
}
IoSimple::Out ioOut(std::cout);
ioOut.setPrecision(12);
cfCollection.write(ioOut);
}
|
