Thyagarajulu Gollapalli

Unit vectors are not defined as normals

4

Hi @lmoresi, In the following snippet, normals are not defined: https://github.com/underworldcode/underworld3/blob/4fe6bf4e6916a6d7b392840b1a8c81780f2fe4b4/src/underworld3/meshing.py#L719 I have two questions regarding this Q1) Is this necessary to define normals? Normals are defined for all other...

Spherical refinement_callback is experiencing some issues

1

It appears that the refinement_callback is the culprit when it comes to issues in spherical models, especially when both refinement and refinement_callback are enabled. However, when I ran the model...

How to remove rotational modes (i.e., Null space) in sphere?

5

Here are the details of null space removal from the benchmark paper. https://gmd.copernicus.org/articles/14/1899/2021/#:~:text=Solving%20the%20linear%20system%20and%20dealing%20with%20null%20spaces Do we also need to remove zero modes from the approximate solution at every iteration?

Error encountered when trying to interpolate mesh variable

10

Hi, In serial: I encountered this error for StructuredQuadBox mesh, however, it works fine for UnstructuredSimplexBox. In parallel: both meshes encounter error. ``` [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC...

Hi, `mesh.CoordinateSystem.rRotN` should gives a matrix to convert coordinates in reference frame (i.e., spherical/Annulus) to native frame (i.e., xyz). However, currently it output a matrix that map N --> r...

RuntimeError while evaluating integration on mesh

3

The mesh requires at least a single variable for integration to function correctly. This is a PETSc limitation. Here is the code to reproduce the error. ```python3 import underworld3 as...

How can I create a 3D surface from point cloud data?

1

Hi, I'm working through the LaGriT tutorials and examples but haven't found a way or a LaGriT command to generate a surface in 3D from point cloud data. I have...

Unexpected kernel crash with essential boundary conditions

5

Hi @lmoresi @julesghub @knepley , Issue: It is currently not possible to implement a essential boundary condition using a mesh variable that is not available in the Stokes object. In...

Spherical benchmarks: mismatch between analytical and computed pressure values

8

Hi @knepley, I am running spherical benchmarks from the following papers: 1. [Paper1](https://gmd.copernicus.org/articles/14/1899/2021/) 2. [Paper2](https://se.copernicus.org/articles/8/1181/2017/) I successfully reproduced the benchmark results from the first paper, but I'm encountering difficulties with...