Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F74856920
cast.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, Jul 30, 02:09
Size
6 KB
Mime Type
text/x-c++
Expires
Thu, Aug 1, 02:09 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
19347218
Attached To
rTAMAAS tamaas
cast.hh
View Options
/**
* @file
*
* @author Lucas Frérot <lucas.frerot@epfl.ch>
*
* @section LICENSE
*
* Copyright (©) 2017 EPFL (Ecole Polytechnique Fédérale de
* Lausanne) Laboratory (LSMS - Laboratoire de Simulation en Mécanique des
* Solides)
*
* Tamaas is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Tamaas is distributed in the hope that it will be useful, but WITHOUT ANY
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
* A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Tamaas. If not, see <http://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#ifndef __CAST_HH__
#define __CAST_HH__
/* -------------------------------------------------------------------------- */
#include "grid_base.hh"
#include "numpy.hh"
#include "surface.hh"
/* -------------------------------------------------------------------------- */
#include <boost/preprocessor/seq.hpp>
#include <pybind11/cast.h>
/* -------------------------------------------------------------------------- */
namespace
pybind11
{
// Format descriptor necessary for correct wrap of tamaas complex type
template
<
typename
T
>
struct
format_descriptor
<
tamaas
::
complex
<
T
>
,
detail
::
enable_if_t
<
std
::
is_floating_point
<
T
>::
value
>>
{
static
constexpr
const
char
c
=
format_descriptor
<
T
>::
c
;
static
constexpr
const
char
value
[
3
]
=
{
'Z'
,
c
,
'\0'
};
static
std
::
string
format
()
{
return
std
::
string
(
value
);
}
};
#ifndef PYBIND11_CPP17
template
<
typename
T
>
constexpr
const
char
format_descriptor
<
tamaas
::
complex
<
T
>
,
detail
::
enable_if_t
<
std
::
is_floating_point
<
T
>::
value
>>::
value
[
3
];
#endif
namespace
detail
{
// declare tamaas complex as a complex type for pybind11
template
<
typename
T
>
struct
is_complex
<
tamaas
::
complex
<
T
>>
:
std
::
true_type
{};
template
<
typename
T
>
struct
is_fmt_numeric
<
tamaas
::
complex
<
T
>
,
detail
::
enable_if_t
<
std
::
is_floating_point
<
T
>::
value
>>
:
std
::
true_type
{
static
constexpr
int
index
=
is_fmt_numeric
<
T
>::
index
+
3
;
};
static
inline
handle
policy_switch
(
return_value_policy
policy
,
handle
parent
)
{
switch
(
policy
)
{
case
return_value_policy
::
copy:
case
return_value_policy
::
move:
return
handle
();
case
return_value_policy
::
automatic_reference:
// happens in python-derived
// classes
case
return_value_policy
::
reference:
return
none
();
case
return_value_policy
::
reference_internal:
return
parent
;
default
:
TAMAAS_EXCEPTION
(
"Policy is not handled"
);
}
}
template
<
class
array
,
typename
T
,
unsigned
int
dim
>
handle
grid_to_python
(
const
tamaas
::
Grid
<
T
,
dim
>&
grid
,
return_value_policy
policy
,
handle
parent
)
{
parent
=
policy_switch
(
policy
,
parent
);
// reusing variable
std
::
vector
<
unsigned
int
>
sizes
(
dim
);
std
::
copy
(
grid
.
sizes
().
begin
(),
grid
.
sizes
().
end
(),
sizes
.
begin
());
if
(
grid
.
getNbComponents
()
!=
1
)
sizes
.
push_back
(
grid
.
getNbComponents
());
return
array
(
sizes
,
grid
.
getInternalData
(),
parent
).
release
();
}
/**
* Type caster for grid classes
* inspired by https://tinyurl.com/y8m47qh3 from T. De Geus
* and pybind11/eigen.h
*/
template
<
template
<
typename
,
unsigned
int
>
class
G
,
typename
T
,
unsigned
int
dim
>
struct
type_caster
<
G
<
T
,
dim
>>
{
using
type
=
G
<
T
,
dim
>
;
using
array_type
=
array_t
<
typename
type
::
value_type
,
array
::
c_style
|
array
::
forcecast
>
;
public
:
PYBIND11_TYPE_CASTER
(
type
,
_
(
"GridWrap<T, dim>"
));
/**
* Conversion part 1 (Python->C++): convert a PyObject into a grid
* instance or return false upon failure. The second argument
* indicates whether implicit conversions should be applied.
*/
bool
load
(
handle
src
,
bool
convert
)
{
if
(
!
convert
&&
!
array_type
::
check_
(
src
))
return
false
;
auto
buf
=
array_type
::
ensure
(
src
);
value
.
move
(
tamaas
::
wrap
::
GridNumpy
<
G
<
T
,
dim
>>
(
buf
));
return
true
;
}
/**
* Conversion part 2 (C++ -> Python): convert a grid instance into
* a Python object. The second and third arguments are used to
* indicate the return value policy and parent object (for
* ``return_value_policy::reference_internal``) and are generally
* ignored by implicit casters.
*/
static
handle
cast
(
const
type
&
src
,
return_value_policy
policy
,
handle
parent
)
{
return
grid_to_python
<
array_type
>
(
src
,
policy
,
parent
);
}
};
/**
* Type caster for GridBase classes
*/
template
<
typename
T
>
struct
type_caster
<
tamaas
::
GridBase
<
T
>>
{
using
type
=
tamaas
::
GridBase
<
T
>
;
using
array_type
=
array_t
<
typename
type
::
value_type
,
array
::
c_style
|
array
::
forcecast
>
;
public
:
PYBIND11_TYPE_CASTER
(
type
,
_
(
"GridBaseWrap<T>"
));
bool
load
(
handle
src
,
bool
convert
)
{
if
(
!
convert
&&
!
array_type
::
check_
(
src
))
return
false
;
auto
buf
=
array_type
::
ensure
(
src
);
value
.
move
(
tamaas
::
wrap
::
GridBaseNumpy
<
T
>
(
buf
));
return
true
;
}
static
handle
cast
(
const
type
&
src
,
return_value_policy
policy
,
handle
parent
)
{
#define GRID_BASE_CASE(unused1, unused2, dim) \
case dim: { \
const tamaas::Grid<T, dim>& conv = dynamic_cast<decltype(conv)>(src); \
return grid_to_python<array_type>(conv, policy, parent); \
}
switch
(
src
.
getDimension
())
{
BOOST_PP_SEQ_FOR_EACH
(
GRID_BASE_CASE
,
~
,
(
1
)(
2
)(
3
));
default
:
return
nullptr
;
}
#undef GRID_BASE_CASE
}
};
/**
* Type caster for surface class
*/
template
<
typename
T
>
struct
type_caster
<
tamaas
::
Surface
<
T
>>
{
using
type
=
tamaas
::
Surface
<
T
>
;
using
array_type
=
array_t
<
typename
type
::
value_type
,
array
::
c_style
|
array
::
forcecast
>
;
public
:
PYBIND11_TYPE_CASTER
(
type
,
_
(
"SurfaceWrap<T, dim>"
));
bool
load
(
handle
src
,
bool
convert
)
{
if
(
!
convert
&&
!
array_type
::
check_
(
src
))
return
false
;
auto
buf
=
array_type
::
ensure
(
src
);
value
.
move
(
tamaas
::
wrap
::
SurfaceNumpy
<
type
>
(
buf
));
return
true
;
}
static
handle
cast
(
const
type
&
src
,
return_value_policy
policy
,
handle
parent
)
{
return
grid_to_python
<
array_type
>
(
src
,
policy
,
parent
);
}
};
}
// namespace detail
}
// namespace pybind11
#endif
// __CAST_HH__
Event Timeline
Log In to Comment