subroutine get_pfcoef_inv_sqrt(this,m,mpow)
implicit none
type(partial_fraction_coef), intent(inout) :: this
integer, intent(in) :: m
integer, intent(in) :: mpow
integer :: ikey
call new_pfcoef_inv_sqrt
if (allocated(this%alp)) deallocate(this%alp)
if (allocated(this%bet)) deallocate(this%bet)
allocate(this%alp(0:m))
allocate(this%bet(1:m))
ikey = mpow + (m-npmin)*npowmax
this%npow = mpow
this%n = m
this%minimaxerror = minimaxerror(ikey)
this%alp(0:m) = alp(ikey,0:m)
this%bet(1:m) = bet(ikey,1:m)
return
end subroutine