Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F82824325
virtualmatrix.cpp
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 13, 15:52
Size
4 KB
Mime Type
text/x-c
Expires
Sun, Sep 15, 15:52 (2 d)
Engine
blob
Format
Raw Data
Handle
19559124
Attached To
rLAMMPS lammps
virtualmatrix.cpp
View Options
/*
*_________________________________________________________________________*
* POEMS: PARALLELIZABLE OPEN SOURCE EFFICIENT MULTIBODY SOFTWARE *
* DESCRIPTION: SEE READ-ME *
* FILE NAME: virtualmatrix.cpp *
* AUTHORS: See Author List *
* GRANTS: See Grants List *
* COPYRIGHT: (C) 2005 by Authors as listed in Author's List *
* LICENSE: Please see License Agreement *
* DOWNLOAD: Free at www.rpi.edu/~anderk5 *
* ADMINISTRATOR: Prof. Kurt Anderson *
* Computational Dynamics Lab *
* Rensselaer Polytechnic Institute *
* 110 8th St. Troy NY 12180 *
* CONTACT: anderk5@rpi.edu *
*_________________________________________________________________________*/
#include "virtualmatrix.h"
#include "matrixfun.h"
#include <cstdlib>
using
namespace
std
;
VirtualMatrix
::
VirtualMatrix
(){
numrows
=
numcols
=
0
;
}
VirtualMatrix
::~
VirtualMatrix
(){
}
int
VirtualMatrix
::
GetNumRows
()
const
{
return
numrows
;
}
int
VirtualMatrix
::
GetNumCols
()
const
{
return
numcols
;
}
double
&
VirtualMatrix
::
operator
()
(
int
i
,
int
j
){
// array access
return
operator_2int
(
i
,
j
);
}
double
VirtualMatrix
::
Get
(
int
i
,
int
j
)
const
{
return
Get_2int
(
i
,
j
);
}
void
VirtualMatrix
::
Set
(
int
i
,
int
j
,
double
value
){
Set_2int
(
i
,
j
,
value
);
}
double
VirtualMatrix
::
BasicGet
(
int
i
,
int
j
)
const
{
return
BasicGet_2int
(
i
,
j
);
}
void
VirtualMatrix
::
BasicSet
(
int
i
,
int
j
,
double
value
){
BasicSet_2int
(
i
,
j
,
value
);
}
void
VirtualMatrix
::
BasicIncrement
(
int
i
,
int
j
,
double
value
){
BasicIncrement_2int
(
i
,
j
,
value
);
}
double
&
VirtualMatrix
::
operator
()
(
int
i
){
return
operator_1int
(
i
);
}
double
VirtualMatrix
::
Get
(
int
i
)
const
{
return
Get_1int
(
i
);
}
void
VirtualMatrix
::
Set
(
int
i
,
double
value
){
Set_1int
(
i
,
value
);
}
double
VirtualMatrix
::
BasicGet
(
int
i
)
const
{
return
BasicGet_1int
(
i
);
}
void
VirtualMatrix
::
BasicSet
(
int
i
,
double
value
){
BasicSet_1int
(
i
,
value
);
}
void
VirtualMatrix
::
BasicIncrement
(
int
i
,
double
value
){
BasicIncrement_1int
(
i
,
value
);
}
double
&
VirtualMatrix
::
operator_1int
(
int
i
)
{
cerr
<<
"Error: single dimensional access is not defined for matrices of type "
<<
GetType
()
<<
endl
;
exit
(
0
);
return
*
(
new
double
);
}
double
VirtualMatrix
::
Get_1int
(
int
i
)
const
{
cerr
<<
"Error: single dimensional access is not defined for matrices of type "
<<
GetType
()
<<
endl
;
exit
(
0
);
return
0.0
;
}
void
VirtualMatrix
::
Set_1int
(
int
i
,
double
value
){
cerr
<<
"Error: single dimensional access is not defined for matrices of type "
<<
GetType
()
<<
endl
;
exit
(
0
);
}
double
VirtualMatrix
::
BasicGet_1int
(
int
i
)
const
{
cerr
<<
"Error: single dimensional access is not defined for matrices of type "
<<
GetType
()
<<
endl
;
exit
(
0
);
return
0.0
;
}
void
VirtualMatrix
::
BasicSet_1int
(
int
i
,
double
value
)
{
cerr
<<
"Error: single dimensional access is not defined for matrices of type "
<<
GetType
()
<<
endl
;
exit
(
0
);
}
void
VirtualMatrix
::
BasicIncrement_1int
(
int
i
,
double
value
){
cerr
<<
"Error: single dimensional access is not defined for matrices of type "
<<
GetType
()
<<
endl
;
exit
(
0
);
}
void
VirtualMatrix
::
Zeros
(){
Const
(
0.0
);
}
void
VirtualMatrix
::
Ones
(){
Const
(
1.0
);
}
ostream
&
VirtualMatrix
::
WriteData
(
ostream
&
c
)
const
{
cerr
<<
"Error: no output definition for matrices of type "
<<
GetType
()
<<
endl
;
exit
(
0
);
}
istream
&
VirtualMatrix
::
ReadData
(
istream
&
c
){
cerr
<<
"Error: no input definition for matrices of type "
<<
GetType
()
<<
endl
;
exit
(
0
);
}
//
// operators and functions
//
ostream
&
operator
<<
(
ostream
&
c
,
const
VirtualMatrix
&
A
){
//output
c
<<
A
.
GetType
()
<<
' '
;
A
.
WriteData
(
c
);
c
<<
endl
;
return
c
;
}
istream
&
operator
>>
(
istream
&
c
,
VirtualMatrix
&
A
){
//input
VirtualMatrix
*
vm
;
int
matrixtype
;
c
>>
matrixtype
;
if
(
MatrixType
(
matrixtype
)
==
A
.
GetType
()
)
A
.
ReadData
(
c
);
else
{
// issue a warning?
cerr
<<
"Warning: During matrix read expected type "
<<
A
.
GetType
()
<<
" and got type "
<<
matrixtype
<<
endl
;
vm
=
NewMatrix
(
matrixtype
);
if
(
!
vm
){
cerr
<<
"Error: unable to instantiate matrix of type "
<<
matrixtype
<<
endl
;
exit
(
0
);
}
vm
->
ReadData
(
c
);
A
.
AssignVM
(
*
vm
);
delete
vm
;
}
return
c
;
}
Event Timeline
Log In to Comment