Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F83977642
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
Fri, Sep 20, 02:15
Size
2 KB
Mime Type
text/x-c++
Expires
Sun, Sep 22, 02:15 (2 d)
Engine
blob
Format
Raw Data
Handle
20788959
Attached To
R9484 sp4e-homework-lars-bertil
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
*
this
->
size
()
+
i
];
}
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