powerflow/solver_nr.h

/* $Id
 * Newton-Raphson solver
 */

#ifndef _SOLVER_NR
#define _SOLVER_NR

#include "complex.h"

typedef struct  {
    int type;               
    unsigned char phases;   
    unsigned char origphases;   
    complex *V;             
    complex *S;             
    complex *Y;             
    complex *I;             
    complex *extra_var;     
    complex *house_var;     
    int *Link_Table;        
    unsigned int Link_Table_Size;   
    double PL[3];           
    double QL[3];           
    double PG[3];           
    double QG[3];           
    double kv_base;         
    double mva_base;        
    double Jacob_A[3];      // Element a in equation (37), which is used to update the Jacobian matrix at each iteration
    double Jacob_B[3];      // Element b in equation (38), which is used to update the Jacobian matrix at each iteration
    double Jacob_C[3];      // Element c in equation (39), which is used to update the Jacobian matrix at each iteration
    double Jacob_D[3];      // Element d in equation (40), which is used to update the Jacobian matrix at each iteration
    unsigned int Matrix_Loc;// Starting index of this object's place in all matrices/equations
    double max_volt_error;  
    char *name;             
    int Parent_Node;                    
    int *Child_Nodes;                   
    unsigned int Child_Node_idx;        
} BUSDATA;

typedef struct {
    complex *Yfrom;         
    complex *Yto;           
    complex *YSfrom;        
    complex *YSto;          
    unsigned char phases;   
    unsigned char origphases;   
    int from;               
    int to;                 
    bool *status;           
    double v_ratio;         
    char *name;             
} BRANCHDATA;

typedef struct Y_NR{
    int row_ind;  
    int col_ind;  
    double Y_value; 
} Y_NR;

typedef struct {
    int row_ind;  
    int col_ind;  
    complex Y[3][3]; 
    char size;      
} Bus_admit;


int64 solver_nr(unsigned int bus_count, BUSDATA *bus, unsigned int branch_count, BRANCHDATA *branch, bool *bad_computations);

#endif

---

GridLAB-D^TM Version 2.0

An open-source project initiated by the US Department of Energy