Thyagarajulu Gollapalli

How about also adding `variable_type` in the `add_bc` function? This way, users can easily identify whether the boundary condition is applied to a velocity or pressure variable. In a way,...

@julesghub @lmoresi Can we also add surface integration feature? Currently we have volume integration over the entire mesh. This can help to compute average quantities on inner/outer surfaces in Annulus/Spherical...

I am using following Stokes settings ```python3 # Stokes settings stokes.tolerance = stokes_tol stokes.petsc_options["ksp_monitor"] = None stokes.petsc_options["snes_type"] = "newtonls" stokes.petsc_options["ksp_type"] = "fgmres" # stokes.petsc_options.setValue("fieldsplit_velocity_pc_type", "mg") stokes.petsc_options.setValue("fieldsplit_velocity_pc_mg_type", "kaskade") stokes.petsc_options.setValue("fieldsplit_velocity_pc_mg_cycle_type", "w") stokes.petsc_options["fieldsplit_velocity_mg_coarse_pc_type"]...

Agree

I was doing something similar to this in #142. I am closing #142. We can track this one instead.

@bknight1, I agree. I referred another issue here because these issues fall under a broader category known as "the scenario where a processor contains zero particles." The existing built-in particle...

Hi @bknight1 @lmoresi , I made a small change in the script, replaced `swarm.data` with `swarm.data[1:]` then it worked. Note: evaluate work for all other np. array datasets except `swarm.data`...

Thanks @bknight1, `uw.function.evalf` works fine both in serial and parallel. However, the numerical solution it provides diverges significantly from the analytical solution. This cause issue in converting some of the...

@knepley could you please point me to the tests you are referring? Thanks

@lmoresi Yes, the issue is still present. I'll investigate further and provide an update on where the problem originates.