underworld3
underworld3 copied to clipboard
I encountered an error indicating a mismatch in particle numbers when using the swarm.advection function. I am using the global_evaluation_rewrite branch with Petsc3.21.5. Here is the error that occurs when...
I’m testing the newest version of Underworld3 and noticed that the Stokes solver seems to enforce only the standard incompressible form of the continuity equation: ∇⋅u=0 In earlier UW3 examples...
I'm getting the following error when the code is doing an uw.function.evaluate: ``` (uw) benknight@MacBook-Pro UW3_timing % mpiexec -np 4 python3 eval_error_test.py [2]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [2]PETSC ERROR:...
This brings the README updates (these include the JOSS publication information). The nasty bug that I introduced into the kdtree algorithm - a fix for this, I mean ! That's...
I have been playing around with scaling in sympy and have got the functionality pretty close to our pint implementation. Before I make it backwards compatible, do we want to...
What should we define as the output of .points_in_domain() when the input is empty (i.e. [])? Currently, it returns False. Do we stick with this behaviour? This is relevant when...
Suppose we add three particles using ```add_particles_with_global_coordinates()```: ```python xy_st = np.array([[0.5, 0.33], [1.5, 0.33], [1.5, 0.50]]) basic_swarm.add_particles_with_coordinates(xy_st, migrate = True) ``` In serial, this is not a problem. In parallel,...
The `ddt.py` routines are Material Derivative classes. There are three of them which are designed to handle the various end member cases: 1. Lagrangian: the classic swarm-based history manager in...
What is the recommended PETSc way to project a field from one DMPlex to another (same bounding box)?
Hi @knepley, I'm trying to project a field from one DMPlex to another using PETSc within Underworld3. Both meshes cover the same bounding box, but have different resolutions. Here's a...