sparse L2,PL2; long *triplet_I,*triplet_J; double *triplet_VAL; long triplet_nz=0,triplet_allocated=0; vector bc_u,bc_v,bc_p,bc;