Unable to compute the curl of a vector field

[dokken]

@raboynics

I’m sorry dokken, even with your code, I am having a hard time. I do not want to do exactly what you’ve done, I believe. I want to set “temp” and “volt” to constant values over the whole mesh, not on particular submeshes.

But why dont you the just do interpolation over the whole mesh.

I suggested compute_incident_entities because you had

@raboynics

def initial_guess_temp():
    return np.asarray([100.0])
    
def initial_guess_volt():
    return np.asarray([0.0])
TempVolt.sub(0).interpolate(initial_guess_temp, facet_indices)
TempVolt.sub(1).interpolate(initial_guess_volt, facet_indices)
# Solve the PDE.
problem = NonlinearProblem(weak_form, TempVolt, J = Jac, bcs=bcs)

in your code, which is simply not correct.

What is wrong with doing
TempVolt.sub(0).interpolate(initial_guess_temp)?

@raboynics

Also, you seem to have altered the order of tdim and fdim in the arguments of compute_incident_entities. Which order would be correct?

The order I have in my code is correct. First dimension is the dimension of the input entities, second dimension is output dimension.
See for instance: What is the difference between compute_incident_entities and entities_to_geometry?