Hi all, I want to check if two forms are equal. Can someone explain why this doesn’t work ? from fenics import * mesh = UnitIntervalMesh(10) V = FunctionSpace(mesh, 'P', 1) u = Function(V) v = TestFunction(V) F1 = dot(2*grad(u), grad(v))*dx F2 = dot(2*grad(u), grad(v))*dx print(F1) print(F2) #…

The reason is that the object created by the scalar–vector product 2*grad(u) has a free index. Every time an object with a free index is created, the global index counter is incremented, to avoid naming conflicts. The reason the index names start at i_8 and not i_0 is that eight indices are alread…

Checking 2 forms are equal

variational formulation

dokken January 21, 2024, 2:53pm 6

Resolved in: Compare two forms - #3 by dokken

Topic		Replies	Views
UFL index notation	1	624	January 29, 2021
Elementwise multiplication in weak formulation variational formulation	4	951	March 1, 2021
Attempting to index an object that is already indexed variational formulation	2	1232	December 23, 2019
How To Matrix of TestFunctions variational formulation	2	360	November 23, 2021
Matrix multiplication between ufl shape (2,2) and (2,) Linear Algebra	1	257	December 29, 2023

Checking 2 forms are equal

Related topics