#define precision_matrix_generation 
#ifdef precision_matrix_generation
      double complex
#else
      complex
#endif
     a a
      
      a = dcmplx(0.d0, 0.d0)
      print *,a
      
      stop
      end