indexing bugfix in weave

This commit is contained in:
James Hensman 2014-10-22 16:43:30 +01:00
parent a87f76df4f
commit f6f36234c6

View file

@ -200,7 +200,7 @@ class Stationary(Kern):
gradq += tmp(n,m)*(X(n,q)-X2(m,q))*(X(n,q)-X2(m,q)); gradq += tmp(n,m)*(X(n,q)-X2(m,q))*(X(n,q)-X2(m,q));
} }
} }
grads[q] = gradq; grads(q) = gradq;
} }
""" """
weave.inline(code, ['tmp', 'X', 'X2', 'grads', 'N', 'M', 'Q'], type_converters=weave.converters.blitz, support_code="#include <math.h>") weave.inline(code, ['tmp', 'X', 'X2', 'grads', 'N', 'M', 'Q'], type_converters=weave.converters.blitz, support_code="#include <math.h>")