#pragma once

// -----------------------------------------------------------------------
// Linear algebra

/* 
 * Vector (double precision)
 * len : dimension of space
 * inc : stride to access in buffer
 * buf : pointer to start of data
 * NOTE: minimum buffer size MUST be len*inc
 */
typedef struct math·Vecd
{
    int     len;
    int     inc;
    double *buf;
} math·Vecd;

double linalg·dotd(math·Vecd a, math·Vecd b);