powerflow/line.cpp

Go to the documentation of this file.

#include <stdlib.h>
#include <stdio.h>
#include <errno.h>
#include <math.h>
#include <iostream>
using namespace std;

#include "node.h"
#include "line.h"

CLASS* line::oclass = NULL;
CLASS* line::pclass = NULL;
CLASS *line_class = NULL;

line::line(MODULE *mod) : link(mod) {
    if(oclass == NULL)
    {
        pclass = link::oclass;
        
        line_class = oclass = gl_register_class(mod,"line",sizeof(line),PC_PRETOPDOWN|PC_BOTTOMUP|PC_POSTTOPDOWN|PC_UNSAFE_OVERRIDE_OMIT);
        if(oclass == NULL)
            GL_THROW("unable to register object class implemented by %s",__FILE__);
        if(gl_publish_variable(oclass,
            PT_INHERIT, "link",
            PT_object, "configuration",PADDR(configuration),
            PT_double, "length[ft]",PADDR(length),
            NULL) < 1) GL_THROW("unable to publish line properties in %s",__FILE__);
    }
}

int line::create()
{
    int result = link::create();

    configuration = NULL;
    length = 0;

    return result;
}

int line::init(OBJECT *parent)
{
    int result = link::init(parent);
    

    node *pFrom = OBJECTDATA(from,node);
    node *pTo = OBJECTDATA(to,node);

    /* check for node nominal voltage mismatch */
    if ((pFrom->nominal_voltage - pTo->nominal_voltage) > fabs(0.001*pFrom->nominal_voltage))
        throw "from and to node nominal voltage mismatch of greater than 0.1%";

    if (solver_method == SM_NR && length == 0.0)
        throw "Newton-Raphson method does not support zero length lines at this time";

    return result;
}

int line::isa(char *classname)
{
    return strcmp(classname,"line")==0 || link::isa(classname);
}


// IMPLEMENTATION OF CORE LINKAGE: line

/* This can be added back in after tape has been moved to commit
EXPORT int commit_line(OBJECT *obj)
{
    if (solver_method==SM_FBS)
    {
        line *plink = OBJECTDATA(obj,line);
        plink->calculate_power();
    }
    return 1;
}
*/
EXPORT int create_line(OBJECT **obj, OBJECT *parent)
{
    try
    {
        *obj = gl_create_object(line::oclass);
        if (*obj!=NULL)
        {
            line *my = OBJECTDATA(*obj,line);
            gl_set_parent(*obj,parent);
            return my->create();
        }
    }
    catch (const char *msg)
    {
        gl_error("create_line: %s", msg);
    }
    return 0;
}

EXPORT TIMESTAMP sync_line(OBJECT *obj, TIMESTAMP t0, PASSCONFIG pass)
{
    line *pObj = OBJECTDATA(obj,line);
    try {
        TIMESTAMP t1 = TS_NEVER;
        switch (pass) {
        case PC_PRETOPDOWN:
            return pObj->presync(t0);
        case PC_BOTTOMUP:
            return pObj->sync(t0);
        case PC_POSTTOPDOWN:
            t1 = pObj->postsync(t0);
            obj->clock = t0;
            return t1;
        default:
            throw "invalid pass request";
        }
    } catch (const char *error) {
        GL_THROW("%s (line:%d): %s", pObj->get_name(), pObj->get_id(), error);
        return 0; 
    } catch (...) {
        GL_THROW("%s (line:%d): %s", pObj->get_name(), pObj->get_id(), "unknown exception");
        return 0;
    }
}

EXPORT int init_line(OBJECT *obj)
{
    line *my = OBJECTDATA(obj,line);
    try {
        return my->init(obj->parent);
    }
    catch (const char *msg)
    {
        GL_THROW("%s (line:%d): %s", my->get_name(), my->get_id(), msg);
        return 0; 
    }
}

EXPORT int isa_line(OBJECT *obj, char *classname)
{
    return OBJECTDATA(obj,line)->isa(classname);
}

---

GridLAB-D^TM Version 2.0

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