#ifndef __PTVECTEUR_4F_C_H__
#define __PTVECTEUR_4F_C_H__

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

#ifdef __cplusplus
extern "C" {
#endif

typedef float PtVect4fc[4];

void setPtVect4f(PtVect4fc vect4fOut, float x, float y, float z, float t);

void setXPtVect4f(PtVect4fc vect4fOut, float x);
void setYPtVect4f(PtVect4fc vect4fOut, float y);
void setZPtVect4f(PtVect4fc vect4fOut, float z);
void setTPtVect4f(PtVect4fc vect4fOut, float t);

void setValueNPtVect4f(PtVect4fc vect4fOut, unsigned int n, float value);
void setValueNTestPtVect4f(PtVect4fc vect4fOut, unsigned int n, float value);

void addPtVect4f(PtVect4fc vect4fOut, const PtVect4fc vect4fIn1, const PtVect4fc vect4fIn2);

void subPtVect4f(PtVect4fc vect4fOut, const PtVect4fc vect4fIn1, const PtVect4fc vect4fIn2);

void oppPtVect4f(PtVect4fc vect4fOut, const PtVect4fc vect4fIn);
void oppSelfPtVect4f(PtVect4fc vect4fInOut);

void multPtVect4f(PtVect4fc vect4fOut, PtVect4fc vect4fIn, float scal);
void multSelfPtVect4f(PtVect4fc vect4fInOut, float scal);

void divPtVect4f(PtVect4fc vect4fOut, const PtVect4fc vect4fIn, float scal);
void divSelfPtVect4f(PtVect4fc vect4fInOut, float scal);

float dotProductPtVect4f(const PtVect4fc vect4fIn1, const PtVect4fc vect4fIn2);

float normePtVect4f(const PtVect4fc vect4fIn);
float normeSQRPtVect4f(const PtVect4fc vect4fIn);

void normaliseSelfPtVect4f(PtVect4fc vect4fInOut);
void normalisePtVect4f(PtVect4fc vect4fOut, const PtVect4fc vect4fIn);
void normaliseVect4f0(PtVect4fc vect4fOut, const PtVect4fc vect4fIn);

void printPtVect4f(const PtVect4fc vect4f);

#ifdef __cplusplus
} /* closing brace for extern "C" */
#endif

#endif