palinear can't process more than 2 millions SNPs in a dosage file

[lckarssen]

This is a 'liftover' of bug #2437 from the R-forge bug tracker.

Error message:

mematrix(nr,nc): cannot allocate memory (1148,2582783)

Input: 1148 individuals, dosage file with 2582783 SNPs.

Cause:
In mematr1.h file lines:

20: data = new (nothrow) DT [ncol*nrow];
165: data = new (nothrow) DT [ncol*nrow];

The ncol and nrow are of type int. The result ncol*nrow is stored in int, which is not enough to represent big number and becomes negative. The allocation fails.

When changed to data = new (nothrow) DT [(unsigned int)ncol*nrow], the memory is allocated successfully.

Another crucial thing is to change nelements variable from int to unsigned int.