/* vide {rdo_sparse_matrix.h} */
#define rdo_sparse_matrix_C_COPYRIGHT "Copyright © 2008 Danillo Pereira and J. Stolfi, UNICAMP"
/* Last edited on 2008-07-28 01:48:37 by stolfi */ 

#define _GNU_SOURCE
#include <stdio.h>
#include <stdlib.h>
#include <assert.h>

#include <vec.h>
#include <affirm.h>
#include <filefmt.h>
#include <nget.h>
#include <fget.h>

#include <sparse_matrix.h>
#include <sparse_matrix_io.h>
#include <sparse_matrix_linalg.h>

#include "rdo_geometry.h"
#include "rdo_sparse_matrix.h"

#define dmat_trivial_elem (0.0)
#define dmat_elem_is_trivial(X) ((X) == 0.0)

sparse_matrix_impl(dmat_t, dmat, double);

void dmat_elem_write(FILE *wr, double *valP) { fprintf(wr, "%24.16e", *valP); }
void dmat_elem_read(FILE *rd, double *valP) { (*valP) = fget_double(rd); }

sparse_matrix_io_impl(dmat_t, dmat, double);

#define dmat_elem_one (1.0)
#define dmat_elem_add(X,Y) ((X)+(Y))
#define dmat_elem_mul(X,Y) ((X)*(Y))

sparse_matrix_linalg_impl(dmat_t, dmat, double);