R/func.R

In VivianePhilipps/marqLevAlgParallel: A Parallelized General-Purpose Optimization Based on Marquardt-Levenberg Algorithm

#' @noRd
func <- function(v,fu,m,nfmax,da,ga,tr,nql)
{
	for(i in 1:m){
		ii <- i*(i+1)/2
		if(v[ii] != 0){
			fu[ii] <- v[ii]+da*((1-ga)*abs(v[ii])+ga*tr)
		}else{
			fu[ii] <- da*ga*tr
		}
	}

	dchole <- .Fortran(C_dchole,fu=as.double(fu),as.integer(m),as.integer(nql),idpos=as.integer(0))
	idpos <- dchole$idpos
	fu <-dchole$fu
	result <- list(idpos=idpos,fu=fu)
	return(result)
}