/* See double_spmat.h */
/* Last edited on 2009-01-17 19:04:03 by stolfi */

#define _GNU_SOURCE
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <string.h>
#include <limits.h>
#include <assert.h>

#include <affirm.h>
#include <fget.h>
#include <nget.h>
#include <jsfile.h>

#include <spmat.h>
#include <spmat_io.h>
#include <spmat_linalg.h>

#include <double_spmat.h>

/* IMPLEMENTATIONS */

spmat_impl(double_spmat_t, double_spmat, double);

spmat_io_impl(double_spmat_t, double_spmat, double);
  
void double_spmat_elem_write(FILE *wr, double *valP)
  { fprintf(wr, "%24.16e", *valP); }

void double_spmat_elem_read(FILE *rd, double *valP)
  { (*valP) = fget_double(rd); }

spmat_linalg_impl(double_spmat_t, double_spmat, double);