Have you considered the following tips and ticks?
You have the line:
w = (E[1, 1] - E[0, 0] + delta) / (E[0, 1] )
If E is some kind of rate of strain tensor then your initial guess needs to be \in C^1 in order to not yield a singularity.
Have you considered the following tips and ticks?
You have the line:
w = (E[1, 1] - E[0, 0] + delta) / (E[0, 1] )
If E is some kind of rate of strain tensor then your initial guess needs to be \in C^1 in order to not yield a singularity.