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	564	January 29, 2021
Matrix multiplication between ufl shape (2,2) and (2,) Linear Algebra	1	151	December 29, 2023
Map ufl expression? variational formulation	13	762	August 10, 2020
Normalized vector of trial function variational formulation	8	1238	June 16, 2019
Index notation on trial function?	6	355	February 1, 2021

Checking 2 forms are equal

Related Topics