stimulate_shear.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python
	import numpy as np


	def compute(input=None, constants=None, **kwargs):
	pos = input
	disp = np.zeros(pos.shape)
	disp[:, 0] = pos[:, 1]*constants["shearstrain"]
	print("\nadditional displacement for {0}: sheared by tau = {1}".format(
	"disp", constants["shearstrain"]))
	print(" {0:e} <= disp <= {1:e} and {2:e} <= position <= {3:e} ".format(
	disp[:, 0].min(), disp[:, 0].max(),
	pos[:, 1].min(), pos[:, 1].max()))
	return disp


	if __name__ == "__main__":
	pass