For PyLith versions 2.2.1 and 2.2.2, could you please comment on the following points:
I gather that PyLith have the option of reporting stresses at the quadrature points, which can be done by commenting out the
CellFilterAvgoption. However, I know that for problems with faults PyLith reports tractions at the fault nodes.
Thus, does PyLith perform any sort of mapping for the stresses when dealing with faults to solve the contact problem ? If so, how is this done ?
How does PyLith computes tractions at the fault nodes for problems involving fault slip ?
Thanks in advance,