[R] quadprog error?

[none]

I forgot to attach the problem data, 'quadprog.Rdata' file, in my prior 
email.


I want to report a following error with quadprog. The solve.QP function
finds a solution to the problem below that violates the last equality
constraint. I tried to solve the same problem using ipop from kernlab
package and get the solution in which all equality constraints are
enforced. I also tried an old version of quadprog, Version: 1.4-11,
Date: 2007-07-12 and my problem is solved correctly.

I have tried to contact Berwin A. Turlach <Berwin.Turlach at gmail.com>
(maintainer for quadprog package) a week ago, with no success.


##############################################################
load(file='quadprog.Rdata')		

	# solve QP using quadprog	
	require(quadprog)
	sol = solve.QP(Dmat, dvec, Amat, bvec, meq)
		x = sol$solution
		check = x %*% Amat - bvec	
		# for some reason last equality constraint is violated
		round(check[1:meq], 4)


	# solve QP using kernlab
	require(kernlab)	
	n = nrow(Amat)
	sv = ipop(c = matrix(dvec), H = Dmat, A = t(Amat[,1:meq]),
			b = bvec[1:meq], l = rep(-1000, n),
			u = rep(1000, n), r = rep(0,meq))
				
		x = primal(sv)
		check = x %*% Amat - bvec	
		# all constraints are ok
		round(check[1:meq], 4)