diff --git a/cbind/base/psb_c_cbase.h b/cbind/base/psb_c_cbase.h
index c6bea6b6..8cb29db2 100644
--- a/cbind/base/psb_c_cbase.h
+++ b/cbind/base/psb_c_cbase.h
@@ -8,11 +8,11 @@ extern "C" {
 
 typedef struct PSB_C_CVECTOR {
   void *cvector;
-} psb_c_cvector; 
+} psb_c_cvector;
 
 typedef struct PSB_C_CSPMAT {
   void *cspmat;
-} psb_c_cspmat; 
+} psb_c_cspmat;
 
 
 /* dense vectors */
@@ -43,7 +43,7 @@ psb_i_t    psb_c_cmat_get_ncols(psb_c_cspmat *mh);
 /* psb_i_t    psb_c_cspasb_opt(psb_c_cspmat *mh, psb_c_descriptor *cdh,  */
 /* 			const char *afmt, psb_i_t upd, psb_i_t dupl); */
 psb_i_t    psb_c_csprn(psb_c_cspmat *mh, psb_c_descriptor *cdh, _Bool clear);
-psb_i_t    psb_c_cmat_name_print(psb_c_cspmat *mh, char *name); 
+psb_i_t    psb_c_cmat_name_print(psb_c_cspmat *mh, char *name);
 
 /* psblas computational routines */
 psb_c_t psb_c_cgedot(psb_c_cvector *xh, psb_c_cvector *yh, psb_c_descriptor *cdh);
@@ -51,15 +51,25 @@ psb_s_t psb_c_cgenrm2(psb_c_cvector *xh, psb_c_descriptor *cdh);
 psb_s_t psb_c_cgeamax(psb_c_cvector *xh, psb_c_descriptor *cdh);
 psb_s_t psb_c_cgeasum(psb_c_cvector *xh, psb_c_descriptor *cdh);
 psb_s_t psb_c_cspnrmi(psb_c_cspmat *ah, psb_c_descriptor *cdh);
-psb_i_t psb_c_cgeaxpby(psb_c_t alpha, psb_c_cvector *xh, 
+psb_i_t psb_c_cgeaxpby(psb_c_t alpha, psb_c_cvector *xh,
 		       psb_c_t beta, psb_c_cvector *yh, psb_c_descriptor *cdh);
-psb_i_t psb_c_cspmm(psb_c_t alpha, psb_c_cspmat *ah, psb_c_cvector *xh, 
+psb_i_t psb_c_cspmm(psb_c_t alpha, psb_c_cspmat *ah, psb_c_cvector *xh,
 		    psb_c_t beta, psb_c_cvector *yh, psb_c_descriptor *cdh);
-psb_i_t psb_c_cspmm_opt(psb_c_t alpha, psb_c_cspmat *ah, psb_c_cvector *xh, 
+psb_i_t psb_c_cspmm_opt(psb_c_t alpha, psb_c_cspmat *ah, psb_c_cvector *xh,
 			psb_c_t beta, psb_c_cvector *yh, psb_c_descriptor *cdh,
 			char *trans, bool doswap);
-psb_i_t psb_c_cspsm(psb_c_t alpha, psb_c_cspmat *th, psb_c_cvector *xh, 
+psb_i_t psb_c_cspsm(psb_c_t alpha, psb_c_cspmat *th, psb_c_cvector *xh,
 		      psb_c_t beta, psb_c_cvector *yh, psb_c_descriptor *cdh);
+/* Additional computational routines */
+psb_i_t psb_c_cgemlt(psb_c_cvector *xh,psb_c_cvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_cgediv(psb_c_cvector *xh,psb_c_cvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_cgediv_check(psb_c_cvector *xh,psb_c_cvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_cgeinv(psb_c_cvector *xh,psb_c_cvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_cgeinv_check(psb_c_cvector *xh,psb_c_cvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_cgeabs(psb_c_cvector *xh,psb_c_cvector *yh,psb_c_cvector *cdh);
+psb_i_t psb_c_cgecmp(psb_c_cvector *xh,psb_c_cvector *ch,psb_c_cvector *zh,psb_c_descriptor *cdh);
+psb_s_t psb_c_cgenrm2_weight(psb_c_cvector *xh,psb_c_cvector *wh,psb_c_descriptor *cdh);
+psb_s_t psb_c_cgenrm2_weightmask(psb_c_cvector *xh,psb_c_cvector *wh,psb_c_cvector *idvh,psb_c_descriptor *cdh);
 #ifdef __cplusplus
 }
 #endif  /* __cplusplus */
diff --git a/cbind/base/psb_c_dbase.h b/cbind/base/psb_c_dbase.h
index 95baca5d..4dde00a8 100644
--- a/cbind/base/psb_c_dbase.h
+++ b/cbind/base/psb_c_dbase.h
@@ -8,11 +8,11 @@ extern "C" {
 
 typedef struct PSB_C_DVECTOR {
   void *dvector;
-} psb_c_dvector; 
+} psb_c_dvector;
 
 typedef struct PSB_C_DSPMAT {
   void *dspmat;
-} psb_c_dspmat; 
+} psb_c_dspmat;
 
 
 /* dense vectors */
@@ -43,7 +43,7 @@ psb_i_t    psb_c_dmat_get_ncols(psb_c_dspmat *mh);
 /* psb_i_t    psb_c_dspasb_opt(psb_c_dspmat *mh, psb_c_descriptor *cdh,  */
 /* 			const char *afmt, psb_i_t upd, psb_i_t dupl); */
 psb_i_t    psb_c_dsprn(psb_c_dspmat *mh, psb_c_descriptor *cdh, _Bool clear);
-psb_i_t    psb_c_dmat_name_print(psb_c_dspmat *mh, char *name); 
+psb_i_t    psb_c_dmat_name_print(psb_c_dspmat *mh, char *name);
 
 /* psblas computational routines */
 psb_d_t psb_c_dgedot(psb_c_dvector *xh, psb_c_dvector *yh, psb_c_descriptor *cdh);
@@ -51,15 +51,27 @@ psb_d_t psb_c_dgenrm2(psb_c_dvector *xh, psb_c_descriptor *cdh);
 psb_d_t psb_c_dgeamax(psb_c_dvector *xh, psb_c_descriptor *cdh);
 psb_d_t psb_c_dgeasum(psb_c_dvector *xh, psb_c_descriptor *cdh);
 psb_d_t psb_c_dspnrmi(psb_c_dvector *xh, psb_c_descriptor *cdh);
-psb_i_t psb_c_dgeaxpby(psb_d_t alpha, psb_c_dvector *xh, 
+psb_i_t psb_c_dgeaxpby(psb_d_t alpha, psb_c_dvector *xh,
 		       psb_d_t beta, psb_c_dvector *yh, psb_c_descriptor *cdh);
-psb_i_t psb_c_dspmm(psb_d_t alpha, psb_c_dspmat *ah, psb_c_dvector *xh, 
+psb_i_t psb_c_dspmm(psb_d_t alpha, psb_c_dspmat *ah, psb_c_dvector *xh,
 		    psb_d_t beta, psb_c_dvector *yh, psb_c_descriptor *cdh);
-psb_i_t psb_c_dspmm_opt(psb_d_t alpha, psb_c_dspmat *ah, psb_c_dvector *xh, 
+psb_i_t psb_c_dspmm_opt(psb_d_t alpha, psb_c_dspmat *ah, psb_c_dvector *xh,
 			psb_d_t beta, psb_c_dvector *yh, psb_c_descriptor *cdh,
 			char *trans, bool doswap);
-psb_i_t psb_c_dspsm(psb_d_t alpha, psb_c_dspmat *th, psb_c_dvector *xh, 
+psb_i_t psb_c_dspsm(psb_d_t alpha, psb_c_dspmat *th, psb_c_dvector *xh,
 		      psb_d_t beta, psb_c_dvector *yh, psb_c_descriptor *cdh);
+/* Additional computational routines */
+psb_i_t psb_c_dgemlt(psb_c_dvector *xh,psb_c_dvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_dgediv(psb_c_dvector *xh,psb_c_dvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_dgediv_check(psb_c_dvector *xh,psb_c_dvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_dgeinv(psb_c_dvector *xh,psb_c_dvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_dgeinv_check(psb_c_dvector *xh,psb_c_dvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_dgeabs(psb_c_dvector *xh,psb_c_dvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_dgecmp(psb_c_dvector *xh,psb_c_dvector *ch,psb_c_dvector *zh,psb_c_descriptor *cdh);
+psb_d_t psb_c_dgenrm2_weight(psb_c_dvector *xh,psb_c_dvector *wh,psb_c_descriptor *cdh);
+psb_d_t psb_c_dgenrm2_weightmask(psb_c_dvector *xh,psb_c_dvector *wh,psb_c_dvector *idvh,psb_c_descriptor *cdh);
+psb_i_t psb_c_dmask(psb_c_dvector *ch,psb_c_dvector *xh,psb_c_dvector *mh, bool t, psb_c_descriptor *cdh);
+
 #ifdef __cplusplus
 }
 #endif  /* __cplusplus */
diff --git a/cbind/base/psb_c_sbase.h b/cbind/base/psb_c_sbase.h
index 5f5c5234..fcdb0cbd 100644
--- a/cbind/base/psb_c_sbase.h
+++ b/cbind/base/psb_c_sbase.h
@@ -8,11 +8,11 @@ extern "C" {
 
 typedef struct PSB_C_SVECTOR {
   void *svector;
-} psb_c_svector; 
+} psb_c_svector;
 
 typedef struct PSB_C_SSPMAT {
   void *sspmat;
-} psb_c_sspmat; 
+} psb_c_sspmat;
 
 
 /* dense vectors */
@@ -43,7 +43,7 @@ psb_i_t    psb_c_smat_get_ncols(psb_c_sspmat *mh);
 /* psb_i_t    psb_c_sspasb_opt(psb_c_sspmat *mh, psb_c_descriptor *cdh,  */
 /* 			const char *afmt, psb_i_t upd, psb_i_t dupl); */
 psb_i_t    psb_c_ssprn(psb_c_sspmat *mh, psb_c_descriptor *cdh, _Bool clear);
-psb_i_t    psb_c_smat_name_print(psb_c_sspmat *mh, char *name); 
+psb_i_t    psb_c_smat_name_print(psb_c_sspmat *mh, char *name);
 
 /* psblas computational routines */
 psb_s_t psb_c_sgedot(psb_c_svector *xh, psb_c_svector *yh, psb_c_descriptor *cdh);
@@ -51,15 +51,26 @@ psb_s_t psb_c_sgenrm2(psb_c_svector *xh, psb_c_descriptor *cdh);
 psb_s_t psb_c_sgeamax(psb_c_svector *xh, psb_c_descriptor *cdh);
 psb_s_t psb_c_sgeasum(psb_c_svector *xh, psb_c_descriptor *cdh);
 psb_s_t psb_c_sspnrmi(psb_c_sspmat *ah, psb_c_descriptor *cdh);
-psb_i_t psb_c_sgeaxpby(psb_s_t alpha, psb_c_svector *xh, 
+psb_i_t psb_c_sgeaxpby(psb_s_t alpha, psb_c_svector *xh,
 		       psb_s_t beta, psb_c_svector *yh, psb_c_descriptor *cdh);
-psb_i_t psb_c_sspmm(psb_s_t alpha, psb_c_sspmat *ah, psb_c_svector *xh, 
+psb_i_t psb_c_sspmm(psb_s_t alpha, psb_c_sspmat *ah, psb_c_svector *xh,
 		    psb_s_t beta, psb_c_svector *yh, psb_c_descriptor *cdh);
-psb_i_t psb_c_sspmm_opt(psb_s_t alpha, psb_c_sspmat *ah, psb_c_svector *xh, 
+psb_i_t psb_c_sspmm_opt(psb_s_t alpha, psb_c_sspmat *ah, psb_c_svector *xh,
 			psb_s_t beta, psb_c_svector *yh, psb_c_descriptor *cdh,
 			char *trans, bool doswap);
-psb_i_t psb_c_sspsm(psb_s_t alpha, psb_c_sspmat *th, psb_c_svector *xh, 
+psb_i_t psb_c_sspsm(psb_s_t alpha, psb_c_sspmat *th, psb_c_svector *xh,
 		      psb_s_t beta, psb_c_svector *yh, psb_c_descriptor *cdh);
+/* Additional computational routines */
+psb_i_t psb_c_sgemlt(psb_c_svector *xh,psb_c_svector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_sgediv(psb_c_svector *xh,psb_c_svector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_sgediv_check(psb_c_svector *xh,psb_c_svector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_sgeinv(psb_c_svector *xh,psb_c_svector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_sgeinv_check(psb_c_svector *xh,psb_c_svector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_sgeabs(psb_c_svector *xh,psb_c_svector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_sgecmp(psb_c_svector *xh,psb_c_svector *ch,psb_c_svector *zh,psb_c_descriptor *cdh);
+psb_s_t psb_c_sgenrm2_weight(psb_c_svector *xh,psb_c_svector *wh,psb_c_descriptor *cdh);
+psb_s_t psb_c_sgenrm2_weightmask(psb_c_svector *xh,psb_c_svector *wh,psb_c_svector *idvh,psb_c_descriptor *cdh);
+psb_i_t psb_c_smask(psb_c_svector *ch,psb_c_svector *xh,psb_c_svector *mh, bool t, psb_c_descriptor *cdh);
 #ifdef __cplusplus
 }
 #endif  /* __cplusplus */
diff --git a/cbind/base/psb_c_zbase.h b/cbind/base/psb_c_zbase.h
index f61f64cf..e3e5e796 100644
--- a/cbind/base/psb_c_zbase.h
+++ b/cbind/base/psb_c_zbase.h
@@ -8,11 +8,11 @@ extern "C" {
 
 typedef struct PSB_C_ZVECTOR {
   void *zvector;
-} psb_c_zvector; 
+} psb_c_zvector;
 
 typedef struct PSB_C_ZSPMAT {
   void *zspmat;
-} psb_c_zspmat; 
+} psb_c_zspmat;
 
 
 /* dense vectors */
@@ -43,7 +43,7 @@ psb_i_t    psb_c_zmat_get_ncols(psb_c_zspmat *mh);
 /* psb_i_t    psb_c_zspasb_opt(psb_c_zspmat *mh, psb_c_descriptor *cdh,  */
 /* 			const char *afmt, psb_i_t upd, psb_i_t dupl); */
 psb_i_t    psb_c_zsprn(psb_c_zspmat *mh, psb_c_descriptor *cdh, _Bool clear);
-psb_i_t    psb_c_zmat_name_print(psb_c_zspmat *mh, char *name); 
+psb_i_t    psb_c_zmat_name_print(psb_c_zspmat *mh, char *name);
 
 /* psblas computational routines */
 psb_z_t psb_c_zgedot(psb_c_zvector *xh, psb_c_zvector *yh, psb_c_descriptor *cdh);
@@ -51,15 +51,25 @@ psb_d_t psb_c_zgenrm2(psb_c_zvector *xh, psb_c_descriptor *cdh);
 psb_d_t psb_c_zgeamax(psb_c_zvector *xh, psb_c_descriptor *cdh);
 psb_d_t psb_c_zgeasum(psb_c_zvector *xh, psb_c_descriptor *cdh);
 psb_d_t psb_c_zspnrmi(psb_c_zspmat *ah, psb_c_descriptor *cdh);
-psb_i_t psb_c_zgeaxpby(psb_z_t alpha, psb_c_zvector *xh, 
+psb_i_t psb_c_zgeaxpby(psb_z_t alpha, psb_c_zvector *xh,
 		       psb_z_t beta, psb_c_zvector *yh, psb_c_descriptor *cdh);
-psb_i_t psb_c_zspmm(psb_z_t alpha, psb_c_zspmat *ah, psb_c_zvector *xh, 
+psb_i_t psb_c_zspmm(psb_z_t alpha, psb_c_zspmat *ah, psb_c_zvector *xh,
 		    psb_z_t beta, psb_c_zvector *yh, psb_c_descriptor *cdh);
-psb_i_t psb_c_zspmm_opt(psb_z_t alpha, psb_c_zspmat *ah, psb_c_zvector *xh, 
+psb_i_t psb_c_zspmm_opt(psb_z_t alpha, psb_c_zspmat *ah, psb_c_zvector *xh,
 			psb_z_t beta, psb_c_zvector *yh, psb_c_descriptor *cdh,
 			char *trans, bool doswap);
-psb_i_t psb_c_zspsm(psb_z_t alpha, psb_c_zspmat *th, psb_c_zvector *xh, 
+psb_i_t psb_c_zspsm(psb_z_t alpha, psb_c_zspmat *th, psb_c_zvector *xh,
 		      psb_z_t beta, psb_c_zvector *yh, psb_c_descriptor *cdh);
+/* Additional computational routines */
+psb_i_t psb_c_zgemlt(psb_c_zvector *xh,psb_c_zvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_zgediv(psb_c_zvector *xh,psb_c_zvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_zgediv_check(psb_c_zvector *xh,psb_c_zvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_zgeinv(psb_c_zvector *xh,psb_c_zvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_zgeinv_check(psb_c_zvector *xh,psb_c_zvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_zgeabs(psb_c_zvector *xh,psb_c_zvector *yh,psb_c_descriptor *cdh);
+psb_i_t psb_c_zgecmp(psb_c_zvector *xh,psb_c_zvector *ch,psb_c_zvector *zh,psb_c_descriptor *cdh);
+psb_d_t psb_c_zgenrm2_weight(psb_c_zvector *xh,psb_c_zvector *wh,psb_c_descriptor *cdh);
+psb_d_t psb_c_zgenrm2_weightmask(psb_c_zvector *xh,psb_c_zvector *wh,psb_c_zvector *idvh,psb_c_descriptor *cdh);
 #ifdef __cplusplus
 }
 #endif  /* __cplusplus */