Errors can occur in the gstat code or during the matrix operations in the meschach matrix library. The cause of and possible solutions to the latter are explained in the second part of this section. Following is a list of the gstat error messages (with exit values).
Other possible errors are:
variogram(a): 1sph(2) + 1sph(2);
then the two models are linearly dependent--fitting their sill will lead to a singularity. Also, if, at a nonlinear fit, the range of a model tends to infinity, the true model may have to be a linear model (having one parameter), but two parameters are being fit for it--they will then be linearly dependent and lead to a singularity during fit. Solution: simplify the variogram model.
matrix not positive definite in function CHfactor()
or
singular matrix in function LDLfactor()
These two error messages may occur (a) during simulation, when an observation falls almost exactly at a simulation location (b) when two observations occur at identical location occur and noaverage was defined in its data definition. Solution to (a): increase the value of zero, to (b): remove the noaverage. Read also the next section.