search for: cs_lu

...mfit package performs well when using no dummy variables (e.g. for year, segment, firm) - When adding dummy variables (year / segment) the systemfit package still provides results; however, if we remove the segment dummy variable and add a firm dummy variable we experience the error message: "cs_lu(A) failed: near-singular A (or out of memory)" - We also tested the systemfit package with firm and year dummy variable using a sub-sample of our data (ca. 10 companies) and the package worked without any error messages - Afterwards we tested the systemfit package using a server (64 gb ram,...

...ations is the same value, and the same for all the other coefficients). When I try to run the model without the restricting equations (C2, C3) it runs just fine, but when I add these two equations I get the error: "Error in solve(xtOmegaInv %*% xMat2, xtOmegaInv %*% yVec, tol = solvetol) : cs_lu(A) failed: near-singular A (or out of memory)" Any ideas on what the problem might be? All the best, Rui Neiva P.S.: I've also posted this question on the Matrix help forum, but since I do not know how active that forum is I've decided to see if anyone in the mailing list would be a...

...there a solver in R that I can use in order to get x without any errors given the structure of A? Thanks for your time. #CODE STARTS HEREA = as(matrix(c(1.5,-1.5,0,-1.5,2.5,-1,0,-1,1),nrow=3,ncol=3),"sparseMatrix")b = matrix(c(-30,40,-10),nrow=3,ncol=1) #solve for x, Error in LU.dgC(a) : cs_lu(A) failed: near-singular A (or out of memory)solve(A,b,sparse=TRUE,tol=.Machine$double.eps) #one x that happens to solve Ax = bx = matrix(c(-10,10,0),nrow=3,ncol=1)A %*% x #Error in lsfit(A, b) : only 3 cases, but 4 variableslsfit(A,b)#solves the system, but fails belowsolve(qr(A, LAPACK=TRUE),b)#E...

...olver in R that I can use in order to get x without any errors given the structure of A? Thanks for your time. > #CODE STARTS HEREA = as(matrix(c(1.5,-1.5,0,-1.5,2.5,-1,0,-1,1),nrow=3,ncol=3),"sparseMatrix")b = matrix(c(-30,40,-10),nrow=3,ncol=1) > #solve for x, Error in LU.dgC(a) : cs_lu(A) failed: near-singular A (or out of memory)solve(A,b,sparse=TRUE,tol=.Machine$double.eps) > #one x that happens to solve Ax = bx = matrix(c(-10,10,0),nrow=3,ncol=1)A %*% x > #Error in lsfit(A, b) : only 3 cases, but 4 variableslsfit(A,b)#solves the system, but fails belowsolve(qr(A, LAPACK=...

...I can use in order to get x > without any errors given the structure of A? Thanks for your time. > #CODE STARTS HEREA = > as(matrix(c(1.5,-1.5,0,-1.5,2.5,-1,0,-1,1),nrow=3,ncol=3),"sparseMatrix")b > = matrix(c(-30,40,-10),nrow=3,ncol=1) > #solve for x, Error in LU.dgC(a) : cs_lu(A) failed: near-singular A (or > out of memory)solve(A,b,sparse=TRUE,tol=.Machine$double.eps) > #one x that happens to solve Ax = bx = matrix(c(-10,10,0),nrow=3,ncol=1)A > %*% x > #Error in lsfit(A, b) : only 3 cases, but 4 variableslsfit(A,b)#solves the > system, but fails belowsolv...

...t I can >use in order to get x without any errors given the structure of A? >Thanks for your time. >#CODE STARTS HEREA = >as(matrix(c(1.5,-1.5,0,-1.5,2.5,-1,0,-1,1),nrow=3,ncol=3),"sparseMatrix")b >= matrix(c(-30,40,-10),nrow=3,ncol=1) >#solve for x, Error in LU.dgC(a) : cs_lu(A) failed: near-singular A (or >out of memory)solve(A,b,sparse=TRUE,tol=.Machine$double.eps) >#one x that happens to solve Ax = bx = >matrix(c(-10,10,0),nrow=3,ncol=1)A %*% x >#Error in lsfit(A, b) : only 3 cases, but 4 variableslsfit(A,b)#solves >the system, but fails belowsolve(qr(...

...olver in R that I can use in order to get x without any errors given the structure of A? Thanks for your time. > #CODE STARTS HEREA = as(matrix(c(1.5,-1.5,0,-1.5,2.5,-1,0,-1,1),nrow=3,ncol=3),"sparseMatrix")b = matrix(c(-30,40,-10),nrow=3,ncol=1) > #solve for x, Error in LU.dgC(a) : cs_lu(A) failed: near-singular A (or out of memory)solve(A,b,sparse=TRUE,tol=.Machine$double.eps) > #one x that happens to solve Ax = bx = matrix(c(-10,10,0),nrow=3,ncol=1)A %*% x > #Error in lsfit(A, b) : only 3 cases, but 4 variableslsfit(A,b)#solves the system, but fails belowsolve(qr(A, LAPACK=...

...nts.o conversion.o cores.o cs_add.o cs_amd.o cs_chol.o cs_cholsol.o cs_compress.o cs_counts.o cs_cumsum.o cs_dfs.o cs_dmperm.o cs_droptol.o cs_dropzeros.o cs_dupl.o cs_entry.o cs_ereach.o cs_etree.o cs_fkeep.o cs_gaxpy.o cs_happly.o cs_house.o cs_ipvec.o cs_leaf.o cs_load.o cs_lsolve.o cs_ltsolve.o cs_lu.o cs_lusol.o cs_malloc.o cs_maxtrans.o cs_multiply.o cs_norm.o cs_permute.o cs_pinv.o cs_post.o cs_print.o cs_pvec.o cs_qr.o cs_qrsol.o cs_randperm.o cs_reach.o cs_scatter.o cs_scc.o cs_schol.o cs_spsolve.o cs_sqr.o cs_symperm.o cs_tdfs.o cs_transpose.o cs_updown.o cs_usolve.o cs_util.o cs_utsolve....

Hi All, Users have contacted me because they can not build MCMCglmm from source. All are using R 3.3.0 on various machines with different compilers gcc (Ubuntu 5.4.0-6ubuntu1~16.04.2) 5.4.0 g++ (Ubuntu 4.8.4-2ubuntu1~14.04.3) 4.8.4 Mac OS X El Capitan (version/compiler unspecified) The issue seems to be with mixing C/C++ with the repeated error: /usr/include/c++/5/bits/cpp_type_traits.h:118:3:

...s -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_ltsolve.c -o cs_ltsolve.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2 -fexceptions -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_lu.c -o cs_lu.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2 -fexceptions -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_lusol.c -o cs_lusol.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -...

...s -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_ltsolve.c -o cs_ltsolve.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2 -fexceptions -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_lu.c -o cs_lu.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2 -fexceptions -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_lusol.c -o cs_lusol.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -...

...solve.o > gcc -I/usr/share/R/include -DNDEBUG??????-fpic??-g -O2 -fstack- > protector-strong -Wformat -Werror=format-security -Wdate-time > -D_FORTIFY_SOURCE=2 -g??-O3 -Wall -pipe -pedantic -std=gnu99??-O3 > -pipe -std=gnu99 -Wno-maybe-uninitialized -Wno-unused-but-set- > variable -c cs_lu.c -o cs_lu.o > gcc -I/usr/share/R/include -DNDEBUG??????-fpic??-g -O2 -fstack- > protector-strong -Wformat -Werror=format-security -Wdate-time > -D_FORTIFY_SOURCE=2 -g??-O3 -Wall -pipe -pedantic -std=gnu99??-O3 > -pipe -std=gnu99 -Wno-maybe-uninitialized -Wno-unused-but-set- > variab...

...-fexceptions -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_ltsolve.c -o cs_ltsolve.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2 -fexceptions -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_lu.c -o cs_lu.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2 -fexceptions -fstack-protector --param=ssp-buffer-size=4 -m64 -mtune=generic -c cs_lusol.c -o cs_lusol.o gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/local/include -fpic -O...

...le -c cs_ltsolve.c -o cs_ltsolve.o gcc -I/usr/share/R/include -DNDEBUG -fpic -g -O2 -fstack-protector-strong -Wformat -Werror=format-security -Wdate-time -D_FORTIFY_SOURCE=2 -g -O3 -Wall -pipe -pedantic -std=gnu99 -O3 -pipe -std=gnu99 -Wno-maybe-uninitialized -Wno-unused-but-set-variable -c cs_lu.c -o cs_lu.o gcc -I/usr/share/R/include -DNDEBUG -fpic -g -O2 -fstack-protector-strong -Wformat -Werror=format-security -Wdate-time -D_FORTIFY_SOURCE=2 -g -O3 -Wall -pipe -pedantic -std=gnu99 -O3 -pipe -std=gnu99 -Wno-maybe-uninitialized -Wno-unused-but-set-variable -c cs_lusol.c -o cs_luso...