main.c

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <time.h>

#include "gebt.h"


int main(int argc, char* argv[])
{
    int i;
    double load;

    // Initializing the Beam
    printf("Initializing the beam...\n");
    Gebt *bm = malloc(sizeof(Gebt));
    initGebt(bm);
    printf("...Done!\n");

    gebt_setTwistFromFile(bm,"Twist.dat");

    // Creating an array with the uniform load for each member
    printf("Creating the loads...\n");
    double *loads = malloc(sizeof(double)*bm->nmemb);
    for(i=0;i<bm->nmemb;i++)
    {
        loads[i] = 5e3;
    }
    gebt_setLoads(bm,loads,1);   // 1 -> Fy (flapwise)
    for(i=0;i<bm->nmemb;i++)
    {
        loads[i] = 10e3;
    }
    gebt_setLoads(bm,loads,2);   // 2 -> Fz (edgewise)
    printf("...Done!\n");

    char fname[32] = "output";
    double meanTime = 0.0;
    clock_t start;
    
    printf("Performing the analysis...\n");
    int ntot=  100;
    for (i=0;i<ntot;i++)
    {
        sprintf(fname,"output_%d_",i);
        printf("Global Step = %d...\n",i);
        start = clock();
        gebt_analysis(bm);
        meanTime += ((double) (clock() - start))/CLOCKS_PER_SEC;
        gebt_writeSolToFile(bm,fname);
        printf("... Ok! \n");
    }

    gebt_writeGebtToFile(bm,"beamFile.out");
    
    free(loads);
    freeGebt(bm);
    free(bm);

    printf("Mean time: %e",meanTime/((double)ntot));
    
    return EXIT_SUCCESS;
}