--- a/oct/comp_pgauss.cc
+++ b/oct/comp_pgauss.cc
@@ -22,7 +22,7 @@
   else
   {
      ComplexMatrix g(L,1);
-     pgauss_cmplx_d(L, w, c_t, c_f, (double (*)[2])g.fortran_vec());
+     pgauss_cmplx_d(L, w, c_t, c_f, (double _Complex*)g.fortran_vec());
   
      return octave_value (g);
   }