I'm trying to get DOLFIN to build against the new UFC branch (the
foo-ufc-geometry branches) and it works mostly ok now, but for a few
demos, something strange happens in the optimized quadrature code.
One example that fails is AdaptivePoisson.ufl/h in
demo/undocumented/auto-adaptive-poisson/. In the new code, the
following fails to build
double G[6];
G[0] = - F1*det*(K[0]*K[2] + K[1]*K[3]);
G[1] = - F1*det*(K[2]*K[2] + K[3]*K[3]);
G[2] = - F1*det*(K[0]*K[0] + K[1]*K[1]);
G[3] = F1*det;
G[4] = - det*(K[0]*(F4*K[0] + F5*K[2]) + K[1]*(F4*K[1] + F5*K[3]));
G[5] = - det*(K[2]*(F4*K[0] + F5*K[2]) + K[3]*(F4*K[1] + F5*K[3]));
because F1 has not been defined.
In the old code, that same piece of code is instead
double G[3];
G[0] = - det*(K_00*K_10 + K_01*K_11);
G[1] = - det*(K_10*K_10 + K_11*K_11);
G[2] = - det*(K_00*K_00 + K_01*K_01);
Any ideas what might have gone wrong?
--
Anders
_______________________________________________
Mailing list: https://launchpad.net/~ffc
Post to : ffc@lists.launchpad.net
Unsubscribe : https://launchpad.net/~ffc
More help : https://help.launchpad.net/ListHelp
