doublevector_utilities.c

Go to the documentation of this file.

#include "turtle.h"
#include "t_datamanipulation.h"
#include "t_utilities.h"
#include "linear_span.h"

#include "doublevector_utilities.h"

#define MAX_REITERATON 1

#define DELTA_MIN epsilon

int linear_comb_doublevector(DOUBLEVECTOR *result,DOUBLEVECTOR *a, DOUBLEVECTOR *b, double ca, double cb) {
        long i;

        if((result->nh!=a->nh) || (result->nh!=b->nh)  || (b->nh!=a->nh) ) {
                printf("Error: in linear_comb vectors result and a and b do not have the same number of elements!!  /n");
                return -1;
        }

        for (i=a->nl;i<=a->nh;i++){
                result->element[i]=ca*a->element[i]+cb*b->element[i];

        }

        return 0;

}




double max_doublevector(DOUBLEVECTOR *v) {
         long j;
         double MK=fabs(v->element[v->nl]);

         for (j=v->nl+1;j<=v->nh;j++){
                 MK=fmax(MK,fabs(v->element[j]));
         }

         return MK;
}

double min_doublevector(DOUBLEVECTOR *v) {
         long j;
         char *function_name="min_doublevector";
         double MK=fabs(v->element[v->nl]);
         for (j=v->nl;j<=v->nh;j++) {
                 if (v->element[j]<0.0) {
                         printf("Warning in function %s, the at %ld position is %le (negative) \n",function_name,j,v->element[j]);
                 }
         }

         for (j=v->nl+1;j<=v->nh;j++){
                 MK=fmin(MK,fabs(v->element[j]));
         }

         return MK;
}