As far as I am aware ufl.cross is only defined in 3D , ref Why is 'fem.assemble_scalar(fem.form(ufl.cross(ufl_vec, sigma(u)*n)[1]*ds(ds_section)))' throwing "ValueError: Index out of bounds." - #4 by dokken
1 Like