Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F88594144
system_factory.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
Sat, Oct 19, 15:51
Size
2 KB
Mime Type
text/x-c
Expires
Mon, Oct 21, 15:51 (2 d)
Engine
blob
Format
Raw Data
Handle
21795044
Attached To
rMUSPECTRE µSpectre
system_factory.hh
View Options
/**
* file system_factory.hh
*
* @author Till Junge <till.junge@epfl.ch>
*
* @date 15 Dec 2017
*
* @brief System factories to help create systems with ease
*
* @section LICENSE
*
* Copyright © 2017 Till Junge
*
* µSpectre is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License as
* published by the Free Software Foundation, either version 3, or (at
* your option) any later version.
*
* µSpectre 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
* General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with GNU Emacs; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* Boston, MA 02111-1307, USA.
*/
#ifndef SYSTEM_FACTORY_H
#define SYSTEM_FACTORY_H
#include "common/common.hh"
#include "common/ccoord_operations.hh"
#include "system/system_base.hh"
#include "fft/projection_finite_strain_fast.hh"
#include "fft/projection_small_strain.hh"
#include "fft/fftw_engine.hh"
#include <memory>
namespace
muSpectre
{
/**
* Create a unique ptr to a Projection operator (with appropriate
* FFT_engine) to be used in a system constructor
*/
template
<
Dim_t
DimS
,
Dim_t
DimM
,
typename
FFT_Engine
=
FFTW_Engine
<
DimS
,
DimM
>>
inline
std
::
unique_ptr
<
ProjectionBase
<
DimS
,
DimM
>>
system_input
(
Ccoord_t
<
DimS
>
resolutions
,
Rcoord_t
<
DimS
>
lengths
,
Formulation
form
)
{
auto
fft_ptr
{
std
::
make_unique
<
FFT_Engine
>
(
resolutions
,
lengths
)};
switch
(
form
)
{
case
Formulation
::
finite_strain:
{
using
Projection
=
ProjectionFiniteStrainFast
<
DimS
,
DimM
>
;
return
std
::
make_unique
<
Projection
>
(
std
::
move
(
fft_ptr
));
break
;
}
case
Formulation
::
small_strain:
{
using
Projection
=
ProjectionSmallStrain
<
DimS
,
DimM
>
;
return
std
::
make_unique
<
Projection
>
(
std
::
move
(
fft_ptr
));
break
;
}
default
:
{
throw
std
::
runtime_error
(
"unknow formulation"
);
break
;
}
}
}
/**
* convenience function to create a system (avoids having to build
* and move the chain of unique_ptrs
*/
template
<
size_t
DimS
,
size_t
DimM
=
DimS
,
typename
System
=
SystemBase
<
DimS
,
DimM
>
,
typename
FFT_Engine
=
FFTW_Engine
<
DimS
,
DimM
>>
inline
System
make_system
(
Ccoord_t
<
DimS
>
resolutions
,
Rcoord_t
<
DimS
>
lengths
,
Formulation
form
)
{
auto
&&
input
=
system_input
<
DimS
,
DimM
,
FFT_Engine
>
(
resolutions
,
lengths
,
form
);
auto
system
{
System
{
std
::
move
(
input
)}};
return
system
;
}
}
// muSpectre
#endif
/* SYSTEM_FACTORY_H */
Event Timeline
Log In to Comment