︠8257871f-246a-4a72-ac10-29ef2fb1f88fi︠
%html
<h2>See <a href="http://groups.google.com/group/sage-devel/browse_thread/thread/95bcfe9c55038d6c/">this</a> thread on <a href="http://groups.google.com/group/sage-devel/">[sage-devel]</a> for details.</h2>

︡69f7c9fc-a4c3-4a20-a6fd-6b8841a07b29︡{"html": "<h2>See <a href=\"http://groups.google.com/group/sage-devel/browse_thread/thread/95bcfe9c55038d6c/\">this</a> thread on <a href=\"http://groups.google.com/group/sage-devel/\">[sage-devel]</a> for details.</h2>"}︡
︠b00eba40-12c9-4d4f-8a4d-3c4ed4c4d998︠
import scipy
from scipy import io
︡98114ff7-d932-4b84-86e3-d2de024241a9︡︡
︠d7544290-bb19-4218-bacb-9eff3048aa4e︠
x=io.loadmat('%s/yodapose.mat'%DATA)#you can change it to yodapose to use the high res version
︡df36aeec-0206-443a-a00e-fc4439a342e0︡︡
︠9d68b25a-b1eb-4015-b8ba-4e2cc99b65ce︠
from sage.plot.plot3d.index_face_set import IndexFaceSet
︡bcacfd66-9c49-4f48-b1d1-4921708f415b︡︡
︠3ed6eace-ec94-44fb-907e-4a8b5bed411f︠
V=x['V']
F3=x['F3']-1
F4=x['F4']-1
︡5dbf97d6-73f9-40eb-8e9d-af1438d3b1a6︡︡
︠1ddd16bc-ae28-4b01-910d-40def9246cba︠
yoda=IndexFaceSet(F3,V,color='green')+IndexFaceSet(F4,V,color='green')
︡d7bd0d6e-e216-4dbe-89f5-ab7efa5af386︡︡
︠d80a63ed-b5bb-4e8d-89c9-0a01f851ee1a︠
yoda.show()
︡c322bd81-6812-44e5-9972-d99c3c18fad9︡︡
︠ee909b3b-1a5c-42a6-8d2e-93d7bd30abee︠
theta = -7*pi/24
R_y = scipy.array([[cos(theta), 0, -sin(theta)],[0, 1, 0],[ sin(theta), 0, cos(theta)]],dtype=float)
Vr = scipy.dot(V,R_y)
︡09ca5f18-0287-47f9-8abe-f5e33a72c737︡︡
︠a612e2a6-8b5e-4f64-b303-cde4ec9ded19︠
yoda_r=IndexFaceSet(F3,Vr,color='green')+IndexFaceSet(F4,Vr,color='green')
︡76be6293-f2f7-4211-a738-60e364bdeeb1︡︡
︠641e76e1-7805-4e50-ae23-36ec11bf6762︠
yoda_r.show()
︡0a73055b-f782-4ede-bffd-5d406817d028︡︡