1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
#include "mtt_Solver.hh"
class Hybrd_Solver : public Solver {
// http://www.netlib.org/minpack/hybrd.f
// used by Octave's fsolve
public:
Hybrd_Solver (sys_ae ae,
const int npar,
const int nu,
const int nx,
const int ny,
const int nyz)
: Solver (ae,npar,nu,nx,ny,nyz)
{
static_ptr = this;
}
static ColumnVector
f_hybrd (const ColumnVector &tryUi);
~Hybrd_Solver (void) {};
protected:
void
Solve (void);
public:
static Hybrd_Solver *static_ptr;
};
|
>
>
|
>
>
>
>
>
|
|
|
|
|
|
|
<
|
|
|
|
|
<
<
<
|
|
|
|
|
|
|
|
|
|
>
|
|
>
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
|
#ifndef MTT_HYBRDSOLVER
#define MTT_HYBRDSOLVER
#include <octave/NLEqn.h>
#include "mtt_AlgebraicSolver.hh"
namespace MTT
{
class Hybrd_Solver : public MTT::AlgebraicSolver
{
// http://www.netlib.org/minpack/hybrd.f
// used by Octave's fsolve
public:
Hybrd_Solver (const int npar,
const int nu,
const int nx,
const int ny,
const int nyz)
: MTT::AlgebraicSolver (npar,nu,nx,ny,nyz);
static ColumnVector
f_hybrd (const ColumnVector &tryUi);
~Hybrd_Solver (void) {};
protected:
void
Solve (void);
public:
static Hybrd_Solver *static_ptr;
};
}
#endif // MTT_HYBRDSOLVER
|