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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
// 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 "LinearPrediction.h"
#include "Io/IoSimple.h"
#include <cstdlib>
#include <fstream>
#include <iostream>
#include <string>
typedef double FieldType;
typedef PsimagLite::LinearPrediction<FieldType> LinearPredictionType;
void usage(const char* progName)
{
std::cerr << "Usage: " << progName << " -f file -l label -p n -q q[>1]\n";
}
int main(int argc, char* argv[])
{
int opt = 0;
PsimagLite::String file = "";
PsimagLite::String label = "";
SizeType p = 0;
SizeType q = 2;
while ((opt = getopt(argc, argv, "f:p:l:q:")) != -1) {
switch (opt) {
case 'f':
file = optarg;
break;
case 'l':
label = optarg;
break;
case 'p':
p = atoi(optarg);
break;
case 'q':
q = atoi(optarg);
break;
default:
usage(argv[0]);
return 1;
}
}
// sanity checks:
if (file == "" || label == "" || p == 0 || q < 2) {
usage(argv[0]);
return 1;
}
PsimagLite::IoSimple::In io(file);
PsimagLite::Vector<FieldType>::Type y;
io.read(y, label);
SizeType n = y.size();
std::cout << "#Found " << n << " points in file " << file << "\n";
LinearPredictionType linearPrediction(y, q);
linearPrediction.predict(q);
for (SizeType i = 0; i < p; i++) {
linearPrediction.linearPredictionfunction(y, q);
linearPrediction.predict(q);
}
for (SizeType i = 0; i < p + n; i++) {
std::cout << linearPrediction(i) << "\n";
}
}
// Test Function in series.txt is: f[i] =
// 0.5*(cos((i-1)*pi/100)+cos((i-1)*pi/20))*exp(-(i-1)/100)
|
