-
Notifications
You must be signed in to change notification settings - Fork 1
/
computeRMSE.cpp
79 lines (58 loc) · 1.7 KB
/
computeRMSE.cpp
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
#include <iostream>
#include "matmethods.hpp"
#include "h5objcpp.hpp"
using namespace std;
int main(int argc, char* argv[]){
char* h5file;
char* outputfile;
H5File file;
int chunkstart, chunknum,Achunksize,Bchunksize;
int totalAchunks,totalBchunks,totalchunks;
int kfold,bsi;
mkl_set_num_threads(1);
int rownum;
if(argc!=9){
cerr<<"Usage:RMSE h5file chunkstart chunknum Achunksize Bchunksize bsi kfold outputfile"<<endl;
cerr<<"argc: "<<argc<<endl;
return(2);
}
h5file = argv[1];
chunkstart = atoi(argv[2]);
chunknum = atoi(argv[3]);
Achunksize = atoi(argv[4]);
Bchunksize = atoi(argv[5]);
bsi = atoi(argv[6]);
kfold = atoi(argv[7]);
outputfile = argv[8];
H5std_string Amatfield = "matA";
H5std_string Bmatfield= "matB";
try{
file = H5File(h5file,H5F_ACC_RDONLY);
}
catch(...)
{
cerr<<"Can't open file: "<<h5file<<endl;
throw 3;
}
vector <hsize_t> sizevec(2);
sizevec = getdims(file,Amatfield);
rownum = sizevec[1];
totalAchunks = ceil((double) sizevec[0]/(double) Achunksize);
sizevec = getdims(file,Bmatfield);
totalBchunks = ceil((double) sizevec[0]/(double) Bchunksize);
totalchunks = totalAchunks*totalBchunks;
cout<<"Total Chunks="<<totalchunks<<endl;
int chunkstop = min(chunkstart+chunknum,totalchunks)+1;
for( int i=chunkstart; i<chunkstop; i++){
if(i>=totalchunks){
return(0);
}
int Achunk = i/totalBchunks;
int Bchunk = i%totalBchunks;
mat A=readh5mat(file,Amatfield,Achunk,Achunksize,0,rownum);
mat B=readh5mat(file,Bmatfield,Bchunk,Bchunksize,0,rownum);
cout<<"Starting chunk: "<<i<<endl;
KfoldCV(A,B,kfold,i,bsi,outputfile,10);
}
return(0);
}