void mexFunction ( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[] );
void assemble ( double *p, int *t, double *u, int nt, int np,
  double *Pr, int *Ir, int *Jc, double *L, double mu );
void assemble_constr ( double *e, double *u, int np, int ne, int n0,
  double *Pr, int *Ir, int *Jc, double *L );
void init_shape ( void );
void i4vec_heap_d ( int n, int a[] );
void i4vec_sort_heap_a ( int n, int a[] );
void localKL ( double *p, int *tt, double u0[], int np, double lK[15][15],
  double lL[15], double mu );
mxArray* sparse_create ( int *t, double *e, int nt, int np, int np0, int ne,
  int ndof );
void sparse_set ( double *Pr, int *Ir, int *Jc, int ri, int rj, double val );