summaryrefslogtreecommitdiffstatshomepage
path: root/src/lib/netlist/solver/nld_solver.h
blob: 986d14f401a39b23dff3a25cf957c6f3e4e86ed0 (plain) (blame)
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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
// license:GPL-2.0+
// copyright-holders:Couriersud
/*
 * nld_solver.h
 *
 */

#ifndef NLD_SOLVER_H_
#define NLD_SOLVER_H_

#include <map>

#include "../nl_base.h"
#include "../plib/pstream.h"
#include "nld_matrix_solver.h"

//#define ATTR_ALIGNED(N) __attribute__((aligned(N)))
#define ATTR_ALIGNED(N) ATTR_ALIGN

// ----------------------------------------------------------------------------------------
// solver
// ----------------------------------------------------------------------------------------

namespace netlist
{
	namespace devices
	{
class NETLIB_NAME(solver);


class matrix_solver_t;

NETLIB_OBJECT(solver)
{
	NETLIB_CONSTRUCTOR(solver)
	, m_fb_step(*this, "FB_step")
	, m_Q_step(*this, "Q_step")
	, m_freq(*this, "FREQ", 48000.0)

	/* iteration parameters */
	, m_gs_sor(*this, "SOR_FACTOR", 1.059)
	, m_method(*this, "METHOD", "MAT_CR")
	, m_accuracy(*this, "ACCURACY", 1e-7)
	, m_gs_loops(*this, "GS_LOOPS",9)              // Gauss-Seidel loops

	/* general parameters */
	, m_gmin(*this, "GMIN", NETLIST_GMIN_DEFAULT)
	, m_pivot(*this, "PIVOT", 0)                    // use pivoting - on supported solvers
	, m_nr_loops(*this, "NR_LOOPS", 250)            // Newton-Raphson loops
	, m_nr_recalc_delay(*this, "NR_RECALC_DELAY", NLTIME_FROM_NS(10).as_double()) // Delay to next solve attempt if nr loops exceeded
	, m_parallel(*this, "PARALLEL", 0)

	/* automatic time step */
	, m_dynamic_ts(*this, "DYNAMIC_TS", 0)
	, m_dynamic_lte(*this, "DYNAMIC_LTE", 1e-5)                     // diff/timestep
	, m_dynamic_min_ts(*this, "DYNAMIC_MIN_TIMESTEP", 1e-6)   // nl_double timestep resolution

	, m_log_stats(*this, "LOG_STATS", 0)   // log statistics on shutdown
	, m_params()
	{
		// internal staff

		connect(m_fb_step, m_Q_step);
	}

	virtual ~NETLIB_NAME(solver)() override;

	void post_start();
	void stop();

	inline nl_double gmin() { return m_gmin(); }

	void create_solver_code(std::map<pstring, pstring> &mp);

	NETLIB_UPDATEI();
	NETLIB_RESETI();
	// NETLIB_UPDATE_PARAMI();

protected:
	logic_input_t m_fb_step;
	logic_output_t m_Q_step;

	param_double_t m_freq;
	param_double_t m_gs_sor;
	param_str_t m_method;
	param_double_t m_accuracy;
	param_int_t m_gs_loops;
	param_double_t m_gmin;
	param_logic_t  m_pivot;
	param_int_t m_nr_loops;
	param_double_t m_nr_recalc_delay;
	param_int_t m_parallel;
	param_logic_t  m_dynamic_ts;
	param_double_t m_dynamic_lte;
	param_double_t m_dynamic_min_ts;

	param_logic_t  m_log_stats;

	std::vector<std::unique_ptr<matrix_solver_t>> m_mat_solvers;
private:

	solver_parameters_t m_params;

	template <std::size_t m_N, std::size_t storage_N>
	std::unique_ptr<matrix_solver_t> create_solver(std::size_t size, const pstring &solvername);
};

	} //namespace devices
} // namespace netlist

#endif /* NLD_SOLVER_H_ */