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
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
/*1:*/
#line 5 "./twod_matrix.cweb"
#include "twod_matrix.h"
#include "tl_exception.h"
/*2:*/
#line 17 "./twod_matrix.cweb"
ConstTwoDMatrix::ConstTwoDMatrix(const TwoDMatrix&m)
:ConstGeneralMatrix(m){}
ConstTwoDMatrix::ConstTwoDMatrix(const TwoDMatrix&m,int first_col,int num)
:ConstGeneralMatrix(m,0,first_col,m.nrows(),num){}
ConstTwoDMatrix::ConstTwoDMatrix(const ConstTwoDMatrix&m,int first_col,int num)
:ConstGeneralMatrix(m,0,first_col,m.nrows(),num){}
ConstTwoDMatrix::ConstTwoDMatrix(int first_row,int num,const TwoDMatrix&m)
:ConstGeneralMatrix(m,first_row,0,num,m.ncols()){}
ConstTwoDMatrix::ConstTwoDMatrix(int first_row,int num,const ConstTwoDMatrix&m)
:ConstGeneralMatrix(m,first_row,0,num,m.ncols()){}
/*:2*/
#line 10 "./twod_matrix.cweb"
;
/*3:*/
#line 34 "./twod_matrix.cweb"
void ConstTwoDMatrix::writeMat(mat_t*fd,const char*vname)const
{
#if MATIO_MAJOR_VERSION > 1 || (MATIO_MAJOR_VERSION == 1 && MATIO_MINOR_VERSION >= 5)
size_t dims[2];
const matio_compression compression= MAT_COMPRESSION_NONE;
#else
int dims[2];
const int compression= COMPRESSION_NONE;
#endif
dims[0]= nrows();
dims[1]= ncols();
double*data= new double[nrows()*ncols()];
for(int j= 0;j<ncols();j++)
for(int i= 0;i<nrows();i++)
data[j*nrows()+i]= get(i,j);
matvar_t*v= Mat_VarCreate(vname,MAT_C_DOUBLE,MAT_T_DOUBLE,2,dims,data,0);
Mat_VarWrite(fd,v,compression);
Mat_VarFree(v);
delete[]data;
}
/*:3*/
#line 11 "./twod_matrix.cweb"
;
/*4:*/
#line 61 "./twod_matrix.cweb"
void TwoDMatrix::copyRow(int from,int to)
{
if(from!=to)
copyRow(ConstTwoDMatrix(*this),from,to);
}
void TwoDMatrix::copyRow(const ConstTwoDMatrix&m,int from,int to)
{
ConstVector fr_row(from,m);
Vector to_row(to,*this);
to_row= fr_row;
}
void TwoDMatrix::addRow(double d,const ConstTwoDMatrix&m,int from,int to)
{
ConstVector fr_row(from,m);
Vector to_row(to,*this);
to_row.add(d,fr_row);
}
/*:4*/
#line 12 "./twod_matrix.cweb"
;
/*5:*/
#line 84 "./twod_matrix.cweb"
void TwoDMatrix::copyColumn(int from,int to)
{
if(from!=to)
copyColumn(ConstTwoDMatrix(*this),from,to);
}
void TwoDMatrix::copyColumn(const ConstTwoDMatrix&m,int from,int to)
{
ConstVector fr_col(m,from);
Vector to_col(*this,to);
to_col= fr_col;
}
void TwoDMatrix::addColumn(double d,const ConstTwoDMatrix&m,int from,int to)
{
ConstVector fr_col(m,from);
Vector to_col(*this,to);
to_col.add(d,fr_col);
}
/*:5*/
#line 13 "./twod_matrix.cweb"
;
/*6:*/
#line 106 "./twod_matrix.cweb"
void TwoDMatrix::save(const char*fname)const
{
FILE*fd;
if(NULL==(fd= fopen(fname,"w"))){
TL_RAISE("Cannot open file for writing in TwoDMatrix::save");
}
for(int row= 0;row<nrows();row++){
for(int col= 0;col<ncols();col++)
fprintf(fd," %20.10g",get(row,col));
fprintf(fd,"\n");
}
fclose(fd);
}
/*:6*/
#line 14 "./twod_matrix.cweb"
;
/*:1*/
|
