Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F79596153
matrix.hh
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Tue, Aug 27, 06:57
Size
2 KB
Mime Type
text/x-c++
Expires
Thu, Aug 29, 06:57 (2 d)
Engine
blob
Format
Raw Data
Handle
20108297
Attached To
R7561 SP4E_HW1
matrix.hh
View Options
#ifndef MATRIX_HH
#define MATRIX_HH
/* ------------------------------------------------------ */
#include "my_types.hh"
#include <algorithm>
#include <array>
#include <tuple>
#include <vector>
/* ------------------------------------------------------ */
template
<
typename
T
>
struct
MatrixIterator
{
MatrixIterator
(
UInt
index
,
UInt
size
,
T
*
ptr
)
:
index
(
index
),
size
(
size
),
ptr
(
ptr
){};
T
&
operator
*
()
{
return
ptr
[
index
];
};
void
operator
++
()
{
index
++
;
};
bool
operator
!=
(
MatrixIterator
&
other
)
{
return
this
->
index
!=
other
.
index
;
};
int
index
,
size
;
T
*
ptr
;
};
/* ------------------------------------------------------ */
template
<
typename
T
>
struct
Matrix
{
Matrix
(){};
Matrix
(
UInt
size
)
{
resize
(
size
);
};
UInt
rows
()
{
return
this
->
size
();
};
UInt
cols
()
{
return
this
->
size
();
};
UInt
size
()
{
return
std
::
sqrt
(
storage
.
size
());
};
void
resize
(
UInt
size
)
{
storage
.
resize
(
size
*
size
);
}
T
&
operator
()(
UInt
i
,
UInt
j
)
{
return
storage
[
j
+
i
*
this
->
size
()];
}
Matrix
&
operator
/=
(
const
T
&
c
)
{
std
::
for_each
(
storage
.
begin
(),
storage
.
end
(),
[
&
c
](
auto
&
v
)
{
v
/=
c
;
});
return
*
this
;
};
T
*
data
()
{
return
&
storage
[
0
];
};
std
::
vector
<
T
>
storage
;
MatrixIterator
<
T
>
begin
()
{
return
MatrixIterator
<
T
>
(
0
,
this
->
size
(),
this
->
data
());
};
MatrixIterator
<
T
>
end
()
{
return
MatrixIterator
<
T
>
(
storage
.
size
(),
this
->
size
(),
this
->
data
());
};
};
/* ------------------------------------------------------ */
template
<
typename
T
>
struct
MatrixIndexIterator
:
public
MatrixIterator
<
T
>
{
MatrixIndexIterator
(
UInt
index
,
UInt
size
,
T
*
ptr
)
:
MatrixIterator
<
T
>
(
index
,
size
,
ptr
){};
std
::
tuple
<
UInt
,
UInt
,
T
&>
operator
*
()
{
int
i
=
this
->
index
/
this
->
size
;
int
j
=
this
->
index
%
this
->
size
;
return
std
::
tuple
<
UInt
,
UInt
,
T
&>
(
i
,
j
,
this
->
ptr
[
this
->
index
]);
};
};
/* ------------------------------------------------------ */
template
<
typename
T
>
struct
IndexedMatrix
{
IndexedMatrix
(
Matrix
<
T
>&
mat
)
:
mat
(
mat
){};
MatrixIndexIterator
<
T
>
begin
()
{
return
MatrixIndexIterator
<
T
>
(
0
,
mat
.
size
(),
mat
.
data
());
};
MatrixIndexIterator
<
T
>
end
()
{
return
MatrixIndexIterator
<
T
>
(
mat
.
storage
.
size
(),
mat
.
size
(),
mat
.
data
());
};
private
:
Matrix
<
T
>&
mat
;
};
/* ------------------------------------------------------ */
template
<
typename
T
>
IndexedMatrix
<
T
>
index
(
Matrix
<
T
>&
mat
)
{
return
IndexedMatrix
<
T
>
(
mat
);
}
/* ------------------------------------------------------ */
template
<
typename
T
>
struct
std
::
iterator_traits
<
MatrixIterator
<
T
>>
{
using
value_type
=
T
;
};
/* ------------------------------------------------------ */
#endif
//MATRIX
Event Timeline
Log In to Comment