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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
|
// license:GPL-2.0+
// copyright-holders:Couriersud
/*
* nld_ms_direct1.h
*
*/
#ifndef NLD_MS_DIRECT2_H_
#define NLD_MS_DIRECT2_H_
#include "nld_ms_direct.h"
#include "nld_solver.h"
namespace netlist
{
namespace devices
{
// ----------------------------------------------------------------------------------------
// matrix_solver - Direct2
// ----------------------------------------------------------------------------------------
template <typename FT>
class matrix_solver_direct2_t: public matrix_solver_direct_t<FT, 2>
{
public:
using float_type = FT;
matrix_solver_direct2_t(netlist_state_t &anetlist, const pstring &name, const solver_parameters_t *params)
: matrix_solver_direct_t<double, 2>(anetlist, name, params, 2)
{}
unsigned vsolve_non_dynamic(const bool newton_raphson) override
{
this->build_LE_A(*this);
this->build_LE_RHS(*this);
const float_type a = this->A(0,0);
const float_type b = this->A(0,1);
const float_type c = this->A(1,0);
const float_type d = this->A(1,1);
const float_type v1 = (a * this->RHS(1) - c * this->RHS(0)) / (a * d - b * c);
const float_type v0 = (this->RHS(0) - b * v1) / a;
std::array<float_type, 2> new_V = {v0, v1};
this->m_stat_calculations++;
const float_type err = (newton_raphson ? this->delta(new_V) : 0.0);
this->store(new_V);
return (err > this->m_params.m_accuracy) ? 2 : 1;
}
};
} //namespace devices
} // namespace netlist
#endif /* NLD_MS_DIRECT2_H_ */
|