70
70
if(!is.complex(y)) y[] <- as.complex(y)
71
71
coef <- matrix(NA_complex_, nrow = p, ncol = ny)
72
coef[qr$pivot,] <- .Call("qr_coef_cmplx", qr, y, PACKAGE = "base")
72
coef[qr$pivot,] <- .Call("qr_coef_cmplx", qr, y, PACKAGE = "base")[1:p,]
73
73
return(if(im) coef else c(coef))