Hi, degree function is erroneous

[lynochka]

Hi,

I am sorry I believe your bearing function is erroneous

p1 = [1.0, 0.0]
p0 = [2.0, -0.1]
from geostatsmodels import utilities
utilities.bearing( p0, p1 )

should be close to 90 degrees but results in 275.71.

I also know that angle tolarance function (lagindices) is wrong but I could not tell exactly where. The reason is, if you do not exclude indices as you do by "indices = np.array([ i for i in indices if i[1] > i[0] ])", the resulting pairs of indicies might be assymetric, check with np.allclose(np.sort(indices[:,0]), np.sort(indices[:,1])). Tell me if you need more information. Do we use the same definition for tolerance?

[wenouyang]

Hi lyonchka,

Have you tested the test_simple_kriging.py? Running kv = kriging.krige(data,kriging.spherical,hs,bw,u,N) gives the following error message. Have you meet the following problem? Thanks.