Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F87909248
new_solvers.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, Oct 15, 16:56
Size
3 KB
Mime Type
text/x-c
Expires
Thu, Oct 17, 16:56 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
21679034
Attached To
rMUSPECTRE µSpectre
new_solvers.hh
View Options
/**
* file new_solvers.hh
*
* @author Till Junge <till.junge@epfl.ch>
*
* @date 24 Apr 2018
*
* @brief Free functions for solving rve problems
*
* @section LICENSE
*
* Copyright © 2018 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 NEW_SOLVERS_H
#define NEW_SOLVERS_H
#include "solvers.hh"
#include "solver/new_solver_base.hh"
#include <Eigen/Dense>
#include <vector>
#include <string>
namespace
muSpectre
{
using
LoadSteps_t
=
std
::
vector
<
Eigen
::
MatrixXd
>
;
/**
* Uses the Newton-conjugate Gradient method to find the static
* equilibrium of a cell given a series of mean applied strains
*/
std
::
vector
<
OptimizeResult
>
newton_cg_dyn
(
Cell
&
cell
,
const
LoadSteps_t
&
load_steps
,
SolverBaseDyn
&
solver
,
Real
newton_tol
,
Real
equil_tol
,
Dim_t
verbose
=
0
);
/**
* Uses the Newton-conjugate Gradient method to find the static
* equilibrium of a cell given a mean applied strain
*/
OptimizeResult
newton_cg_dyn
(
Cell
&
cell
,
const
Eigen
::
Ref
<
Eigen
::
MatrixXd
>
load_step
,
SolverBaseDyn
&
solver
,
Real
newton_tol
,
Real
equil_tol
,
Dim_t
verbose
=
0
)
{
LoadSteps_t
load_steps
{
load_step
};
return
newton_cg_dyn
(
cell
,
load_steps
,
solver
,
newton_tol
,
equil_tol
,
verbose
).
front
();
}
/* ---------------------------------------------------------------------- */
/**
* Uses the method proposed by de Geus method to find the static
* equilibrium of a cell given a series of mean applied strains
*/
std
::
vector
<
OptimizeResult
>
de_geus_dyn
(
Cell
&
cell
,
const
LoadSteps_t
&
load_steps
,
SolverBaseDyn
&
solver
,
Real
newton_tol
,
Real
equil_tol
,
Dim_t
verbose
=
0
);
/* ---------------------------------------------------------------------- */
/**
* Uses the method proposed by de Geus method to find the static
* equilibrium of a cell given a mean applied strain
*/
OptimizeResult
de_geus_dyn
(
Cell
&
cell
,
const
Eigen
::
Ref
<
Eigen
::
MatrixXd
>
load_step
,
SolverBaseDyn
&
solver
,
Real
newton_tol
,
Real
equil_tol
,
Dim_t
verbose
=
0
){
return
de_geus_dyn
(
cell
,
LoadSteps_t
{
load_step
},
solver
,
newton_tol
,
equil_tol
,
verbose
)[
0
];
}
}
// muSpectre
#endif
/* NEW_SOLVERS_H */
Event Timeline
Log In to Comment