vector_operations.h

/* vim: ft=cpp
 * */
#pragma once
#include <vector_types.h> // float4
#include "config.h"
//__device__ float4 operator*(float4 a, float4 b);
//__device__ float4 operator-(float4 a, float4 b);
//__device__ float4 operator-(float4 a);
//__device__ float4 operator+(float4 a, float4 b);
//__device__ float4 operator/(float4 a, float k);
static __device__ float4 operator*(float4 a, float4 b) {
    return make_float4(a.x*b.x,
                       a.y*b.y,
                       a.z*b.z,
                       0);
}
static __device__ float4 operator-(float4 a, float4 b) {
    return make_float4(a.x-b.x,
                       a.y-b.y,
                       a.z-b.z,
                       0);
}
static __device__ float4 operator-(float4 a) {
    return make_float4(-a.x,
                       -a.y,
                       -a.z,
                       0);
}
static __device__ float4 operator+(float4 a, float4 b) {
    return make_float4(a.x+b.x,
                       a.y+b.y,
                       a.z+b.z,
                       0);
}
static __device__ float4 operator/(float4 a, float k) {
    return make_float4(a.x/k,
                       a.y/k,
                       a.z/k,
                       0);
}
static __device__ float l1_float4 (float4 a) {
    return ( fabsf (a.x) +
             fabsf (a.y) +
             fabsf (a.z))*0.3333333f;

}
static __device__ float l2_float4 (float4 a) {
    return sqrtf( pow2 (a.x) +
             pow2 (a.y) +
             pow2 (a.z));

}