mat_wpow Subroutine

subroutine mat_wpow(in_real,in_imag,out_real,out_imag,power_real,power_imag)

! ident_21="@(#)M_LA::mat_wpow(3fp): y = x**p

doubleprecision,intent(in)  :: in_real
doubleprecision,intent(in)  :: in_imag
doubleprecision,intent(in)  :: power_real
doubleprecision,intent(in)  :: power_imag
doubleprecision,intent(out) :: out_real
doubleprecision,intent(out) :: out_imag
complex(kind=real64)        :: t
   ! placeholder method, just using Fortran
   t=cmplx(in_real,in_imag,kind=real64)**cmplx(power_real,power_imag,kind=real64)
   out_real=real(t,kind=real64)
   out_imag=aimag(t)
end subroutine mat_wpow