Mirrors/VILLASnode
1
0
Fork 0
mirror of https://git.rwth-aachen.de/acs/public/villas/node/ synced 2025-03-30 00:00:11 +01:00
Code Issues Releases Wiki Activity

cbuilder: coding style

Browse source
This commit is contained in:
Steffen Vogel 2017-02-12 14:44:29 -03:00
parent dea0f1d61c
commit ac7a87b8fc
1 changed files with 82 additions and 89 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

171
plugins/models/simple_circuit.c
View file

@ -1,114 +1,113 @@
// This is c-code for CBuilder component for Subsystem 2
// Solver used as in RTDS: Resistive companion (Dommel's algo)
// Subsystem 1 is modelled in RSCAD
//
//% Circuit topology
// %
//% *** Subsystem 1 (SS1) *** % *** Subsystem 2 (SS2) ***
//% %
//% |---------| |---------| %
//% |---------| R1 |-------| L1 |----%-------|---------|
//% | |---------| |---------| % | |
//% | % | |
//% ----- % ----- -----
//% | + | % | | | |
//% | E | % |C2 | | R2|
//% | - | % | | | |
//% ----- % ----- -----
//% | % | |
//% |------------------------------------------%------------------
// %
// %
/** This is c-code for CBuilder component for Subsystem 2
* Solver used as in RTDS: Resistive companion (Dommel's algo)
* Subsystem 1 is modelled in RSCAD
*
* % Circuit topology
* %
* % *** Subsystem 1 (SS1) *** % *** Subsystem 2 (SS2) ***
* % %
* % |---------| |---------| %
* % |---------| R1 |-------| L1 |----%-------|---------|
* % | |---------| |---------| % | |
* % | % | |
* % ----- % ----- -----
* % | + | % | | | |
* % | E | % |C2 | | R2|
* % | - | % | | | |
* % ----- % ----- -----
* % | % | |
* % |------------------------------------------%------------------
* %
* %
*/
// -----------------------------------------------
// Variables declared here may be used as parameters
// inputs or outputs
// The have to match with whats in Subsystem.h
// -----------------------------------------------
/* Variables declared here may be used as parameters inputs or outputs
* The have to match with whats in Subsystem.h */
#if defined(VILLAS) || SECTION == INPUTS
double IntfIn;
#endif
#endif /* defined(VILLAS) || SECTION == INPUTS */
#if defined(VILLAS) || SECTION == OUTPUTS
double IntfOut;
#endif
#endif /* defined(VILLAS) || SECTION == OUTPUTS */
#if defined(VILLAS) || SECTION == PARAMS
double R2; // Resistor [Ohm] in SS2
double C2; // Capacitance [F] in SS2
#endif
double R2; /**< Resistor [Ohm] in SS2 */
double C2; /**< Capacitance [F] in SS2 */
#endif /* defined(VILLAS) || SECTION == PARAMS */
// -----------------------------------------------
// Variables declared here may be used in both the
// RAM: and CODE: sections below.
// -----------------------------------------------
/* Variables declared here may be used in both the RAM: and CODE: sections below. */
#if defined(VILLAS) || SECTION == STATIC
double dt;
double GR2, GC2; //Inductances of components
double GnInv; //Inversion of conductance matrix (here only scalar)
double vC2Hist, iC2Hist, AC2; // history meas. and current of dynamic elements
double Jn; //source vector in equation Gn*e=Jn
double eSS2; //node voltage solution
#endif
double dt;
double GR2, GC2; /**< Inductances of components */
double GnInv; /**< Inversion of conductance matrix (here only scalar) */
double vC2Hist, iC2Hist, AC2; /**< History meas. and current of dynamic elements */
double Jn; /**< Source vector in equation Gn*e=Jn */
double eSS2; /**< Node voltage solution */
#endif /* defined(VILLAS) || SECTION == STATIC */
// -----------------------------------------------
// This section should contain any 'c' functions
// to be called from the RAM section (either
// RAM_PASS1 or RAM_PASS2). Example:
//
// static double myFunction(double v1, double v2)
// {
// return(v1*v2);
// }
// -----------------------------------------------
/* This section should contain any 'c' functions
* to be called from the RAM section (either
* RAM_PASS1 or RAM_PASS2). Example:
*
* static double myFunction(double v1, double v2)
* {
* return(v1*v2);
* }
*/
#if defined(VILLAS) || SECTION == RAM_FUNCTIONS
/* Nothing here */
#endif
// -----------------------------------------------
// Place C code here which computes constants
// required for the CODE: section below. The C
// code here is executed once, prior to the start
// of the simulation case.
// -----------------------------------------------
/* Nothing here */
#endif /* defined(VILLAS) || SECTION == RAM_FUNCTIONS */
/* Place C code here which computes constants
* required for the CODE: section below. The C
* code here is executed once, prior to the start
* of the simulation case.
*/
#if defined(VILLAS) || SECTION == RAM
void simple_circuit_ram() {
GR2 = 1/R2;
GC2 = 2*C2/dt; //trapezoidal rule
GnInv = 1/(GR2+GC2); //eq. conductance (inverted)
vC2Hist = 0.0; //Voltage over C2 in previous time step
iC2Hist = 0.0; //Current through C2 in previous time step
GR2 = 1/R2;
GC2 = 2*C2/dt; /**< Trapezoidal rule */
GnInv = 1/(GR2+GC2); /**< eq. conductance (inverted) */
vC2Hist = 0.0; /**< Voltage over C2 in previous time step */
iC2Hist = 0.0; /**< Current through C2 in previous time step */
}
#endif
#endif /* defined(VILLAS) || SECTION == RAM */
// -----------------------------------------------
// Place C code here which runs on the RTDS. The
// code below is entered once each simulation
// step.
// -----------------------------------------------
// Place C code here which runs on the RTDS. The
// code below is entered once each simulation
// step.
// -----------------------------------------------
#if defined(VILLAS) || SECTION == CODE
void simple_circuit_code() {
//Update source vector
/* Update source vector */
AC2 = iC2Hist+vC2Hist*GC2;
Jn = IntfIn+AC2;
//Solution of the equation Gn*e=Jn;
eSS2 = GnInv*Jn;
//Post step -> calculate the voltage and current for C2 for next step and set interface output
/* Solution of the equation Gn*e=Jn; */
eSS2 = GnInv*Jn;
/* Post step -> calculate the voltage and current for C2 for next step and set interface output */
vC2Hist= eSS2;
iC2Hist = vC2Hist*GC2-AC2;
IntfOut = eSS2;
}
#endif
#endif /* defined(VILLAS) || SECTION == CODE */
// -----------------------------------------------
// Interface to VILLASnode
// -----------------------------------------------
/* Interface to VILLASnode */
#if defined(VILLAS)
#include <villas/nodes/cbuilder.h>
@ -119,9 +118,7 @@ double getTimeStep()
return dt;
}
// -----------------------------------------------
// Place C code here which intializes parameters
// -----------------------------------------------
/** Place C code here which intializes parameters */
int simple_circuit_init(struct cbuilder *cb)
{
if (cb->paramlen < 2)
@ -136,9 +133,7 @@ int simple_circuit_init(struct cbuilder *cb)
return 0; /* success */
}
// -----------------------------------------------
// Place C code here reads model outputs
// -----------------------------------------------
/** Place C code here reads model outputs */
int simple_circuit_read(float outputs[], int len)
{
if (len < 1)
@ -149,9 +144,7 @@ int simple_circuit_read(float outputs[], int len)
return 1; /* 1 value per sample */
}
// -----------------------------------------------
// Place C code here which updates model inputs
// -----------------------------------------------
/** Place C code here which updates model inputs */
int simple_circuit_write(float inputs[], int len)
{
if (len < 1)
@ -177,4 +170,4 @@ static struct plugin p = {
REGISTER_PLUGIN(&p)
#endif
#endif /* defined(VILLAS) */