powerflow/triplexline_test.h

#ifndef _TRIPLEXLINE_TEST_H
#define _TRIPLEXLINE_TEST_H

#include "line.h"

#ifndef _NO_CPPUNIT

class triplex_line_tests : public test_helper
{
    
    
public:
    triplex_line_tests(){}
    void setUp(){
    }

    void tearDown(){
        local_callbacks->remove_objects();
    }

    /*void test_create(){
        CLASS *cl = get_class_by_name("sample");
        CPPUNIT_ASSERT(cl != NULL);
        obj = gl_create_object(cl,sizeof(sample_obj)); // Sample_obj is the name of the cpp class or struct.
        CPPUNIT_ASSERT(obj != NULL);
        
    }*/

    void test_triplex_line(){
        OBJECT *node1,*node2; // Triplex Node objects
        OBJECT *tpl1; // Triplex line object
        // Test values
        complex* node1_phA_I_test = new complex(159.91,-79.788);
        complex* node1_phB_I_test = new complex(-199.4,106.36);
        complex* node1_phC_I_test = new complex(39.494,-26.569);

        complex* node2_phA_V_test = new complex(115.44,-.26192);
        complex* node2_phB_V_test = new complex(122.45,-2.542);
        complex* node2_phC_V_test = new complex(0,0);

        CLASS *cl = get_class_by_name("triplex_node");
        node1 = gl_create_object(cl,sizeof(triplex_node));
        OBJECTDATA(node1,triplex_node)->create();
        cl = get_class_by_name("triplex_node");
        node2 = gl_create_object(cl,sizeof(triplex_node));
        OBJECTDATA(node2,triplex_node)->create();
        cl = get_class_by_name("triplex_line");
        tpl1 = gl_create_object(cl,sizeof(triplex_line));
        OBJECTDATA(tpl1,triplex_line)->create();

        OBJECTDATA(tpl1,triplex_line)->from = node1;
        OBJECTDATA(tpl1,triplex_line)->to = node2;
        OBJECTDATA(tpl1,triplex_line)->length = 100;

        // Next we set up the parameters needed for the test
        OBJECTDATA(tpl1,triplex_line)->phases = PHASE_ABCN;
        OBJECTDATA(node1,triplex_node)->phases = PHASE_ABCN;
        OBJECTDATA(node2,triplex_node)->phases = PHASE_ABCN;
        OBJECTDATA(node1,triplex_node)->phaseA_V = OBJECTDATA(node2,triplex_node)->phaseA_V = complex(115.44,-0.26192);
        OBJECTDATA(node2,triplex_node)->phaseA_I = complex(159.91,-79.788);
        OBJECTDATA(node1,triplex_node)->phaseB_V = OBJECTDATA(node2,triplex_node)->phaseB_V = complex(122.45,-2.5421);
        OBJECTDATA(node2,triplex_node)->phaseB_I = complex(-199.4,106.36);
        OBJECTDATA(node1,triplex_node)->phaseC_V = OBJECTDATA(node2,triplex_node)->phaseC_V = complex(0,0);
        OBJECTDATA(node2,triplex_node)->phaseC_I = complex(39.494,-26.569);
        
        cl = get_class_by_name("triplex_line_conductor");
        OBJECT *olc_a = gl_create_object(cl,sizeof(triplex_line_conductor)); // L1 conductor
        OBJECTDATA(olc_a,triplex_line_conductor)->geometric_mean_radius = 0.0111;
        OBJECTDATA(olc_a,triplex_line_conductor)->resistance = 0.97;

        OBJECT *olc_b = gl_create_object(cl,sizeof(triplex_line_conductor)); // L2 conductor
        OBJECTDATA(olc_b,triplex_line_conductor)->geometric_mean_radius = 0.0111;
        OBJECTDATA(olc_b,triplex_line_conductor)->resistance = 0.97;

        OBJECT *olc_c = gl_create_object(cl,sizeof(triplex_line_conductor)); // LN conductor
        OBJECTDATA(olc_c,triplex_line_conductor)->geometric_mean_radius = 0.0111;
        OBJECTDATA(olc_c,triplex_line_conductor)->resistance = 0.97;

        // Do I need to set up the other conductors and set them to INF?
        cl = get_class_by_name("triplex_line_configuration");
        OBJECT *line_config = gl_create_object(cl,sizeof(triplex_line_configuration));
        triplex_line_configuration *lc = OBJECTDATA(line_config,triplex_line_configuration);
        //lc->line_spacing = ls;
        lc->phaseA_conductor = olc_a;
        lc->phaseB_conductor = olc_b;
        lc->phaseC_conductor = olc_c;
        lc->phaseN_conductor = NULL;
        lc->ins_thickness = 0.08; // replace with test value
        lc->diameter = 0.368; // replace with test value
        
        OBJECTDATA(tpl1,triplex_line)->configuration = line_config;
        
        // Now we call init on the objects
        CPPUNIT_ASSERT(local_callbacks->init_objects() != FAILED);

        // Rank?
        CPPUNIT_ASSERT(local_callbacks->setup_test_ranks() != FAILED);

        // Run the test
        local_callbacks->sync_all(PC_PRETOPDOWN);
        local_callbacks->sync_all(PC_BOTTOMUP);

        // Check the bottom-up pass values
        double rI_1 = OBJECTDATA(node1,triplex_node)->get_phaseA_I_in().Re(); // Real part of the current at node 1
        double iI_1 = OBJECTDATA(node1,triplex_node)->get_phaseA_I_in().Im(); // Imaginary part of the current at node 1

        // Get the percentage of the difference between the calculated value and the test value
        double validate_re = 1 - (rI_1 <= node1_phA_I_test->Re() ? rI_1 / node1_phA_I_test->Re() : node1_phA_I_test->Re() /rI_1);
        double validate_im = 1 - (iI_1 <= node1_phA_I_test->Im() ? iI_1 / node1_phA_I_test->Im() : node1_phA_I_test->Im() /iI_1);

        CPPUNIT_ASSERT(validate_re < VALIDATE_THRESHOLD);
        CPPUNIT_ASSERT(validate_im < VALIDATE_THRESHOLD);

        
        OBJECTDATA(node1,triplex_node)->phaseA_V = -1024.9;
        OBJECTDATA(node1,triplex_node)->phaseA_V.SetImag(2108.2);
        local_callbacks->sync_all(PC_POSTTOPDOWN);
        // check the top down values

        double rV_2 = OBJECTDATA(node2,triplex_node)->phaseA_V.Re();
        double iV_2 = OBJECTDATA(node2,triplex_node)->phaseA_V.Im();

        validate_re = 1 - (rV_2 <= node2_phA_V_test->Re() ? rV_2 / node2_phA_V_test->Re() : node2_phA_V_test->Re() /rV_2);
        validate_im = 1 - (iV_2 <= node2_phA_V_test->Im() ? iV_2 / node2_phA_V_test->Im() : node2_phA_V_test->Im() /iV_2);
        CPPUNIT_ASSERT(validate_re < VALIDATE_THRESHOLD);
        CPPUNIT_ASSERT(validate_im < VALIDATE_THRESHOLD);

    }

    /*
     * This section creates the suite() method that will be used by the
     * CPPUnit testrunner to execute the tests that we have registered.
     * This section needs to be in the .h file
     */
    CPPUNIT_TEST_SUITE(triplex_line_tests);
    /*
     * For each test method defined above, we should have a separate
     * CPPUNIT_TEST() line.
     */
    CPPUNIT_TEST(test_triplex_line);
    CPPUNIT_TEST_SUITE_END();

    
};

#endif
#endif

---

GridLAB-D^TM Version 1.0

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