def function(out, val_qp, mat, vg, fmode): if fmode == 2: out[:] = dot_sequences(mat, val_qp) status = 0 else: status = terms.de_cauchy_stress(out, val_qp, mat, vg, fmode) out *= -1.0 return status
def function(out, coef, strain, mat, vg, fmode): if fmode == 2: out[:] = dot_sequences(mat, strain) status = 0 else: status = terms.de_cauchy_stress(out, strain, mat, vg, fmode) if coef is not None: out *= coef return status