matlab - Debugging 24 symmetric operators -


i have written following code calculate disorientation between 2 points in large dataset using 24 crystal symmetry operators. code seems work fine though end result not right. have huge dataset containing euler angles each points. find misorientation between 1st point , 2nd point , between 2 , 3rd , on.

i want end result contain corresponding misorientaion angles each 2 data points. following full code better understanding of requirement. 

la=[phi1 phi phi2];  function gi=get_gi(pvec) %pvec 1x3 vector of [phi1 phi phi2] g_11=((cosd(pvec(1)).*cosd(pvec(3)))-(sind(pvec(1)).*sind(pvec(3)).*cosd(pvec(2)))); g_12=((sind(pvec(1)).*cosd(pvec(3)))+(cosd(pvec(1)).*sind(pvec(3)).*cosd(pvec(2)))); g_13= (sind(pvec(3)).*sind(pvec(2))); g_21 =((-cosd(pvec(1)).*sind(pvec(3)))-(sind(pvec(1)).*cos(pvec(3)).*cos(pvec(2))));  g_22 = ((-sin(pvec(1)).*sind(pvec(3)))+(cosd(pvec(1)).*cosd(pvec(3)).*cosd(pvec(2)))); g_23 = (cosd(pvec(3)).*sind(pvec(2))); g_31 = (sind(pvec(1)).* sind(pvec(2))); g_32 = -cosd(pvec(1)).* sind(pvec(2)); g_33 = cosd(pvec(2));   gi =[g_11 g_12 g_13;g_21 g_22 g_23;g_31 g_32 g_33];   f = [1 1 1 -1 1 -1 -1 -1 1 1 -1 -1]; l= [1 1 1]; i=1:3:10     l1= [f(i) 0 0;0 f(i+1) 0;0 0 f(i+2)]; l2= [0 f(i) 0;0 0 f(i+1);f(i+2) 0 0]; l3= [0 0 f(i);f(i+1) 0 0;0 f(i+2) 0]; l4= -[0 0 f(i);0 f(i+1) 0;f(i+2) 0 0]; l5= -[0 f(i) 0;f(i+1) 0 0;0 0 f(i+2)]; l6= -[f(i) 0 0;0 0 f(i+1);0 f(i+2) 0]; l=[l;l1;l2;l3;l4;l5;l6]; end k=1; t=1;   m=1:(length(a)-1) i=2:3:71     j=2:3:71         g_r= (get_gi(la(m,:)*l(i:i+2,1:3))*(inv(get_gi(la(m+1,:)))*inv(l(j:j+2,1:3))));         tr(k)= g_r(1,1) + g_r(2,2) +g_r(3,3);         angle(k) = real(acosd((tr(k)-1)/2));         k=k+1;     end     angle(angle==0)=360;     del_theta=min(angle)     del(t)=del_theta;     t=t+1; end


Comments

Post a Comment

Popular posts from this blog

SPSS keyboard combination alters encoding -

my question not have programming practical issue on spss. accident keyword combination cannot remember because did not intend press changed encoding in spss , text in data in different language english shown small boxes. nevertheless, when kind of analysis performed different language text shown in output. does know reason? tried change encoding did not produce result. first, @ edit > options , see whether setting unicode or locale. might fix it. if have text in multiple languages, want unicode. otherwise, may have specific locale set wrong. can change set locale, need know locale want use.
Read more

Add new record to the table by click on the button in Microsoft Access -

in ms access form want implement separate button, adds new record table. in order added button , attached button event: private sub btnaddrec_click() refresh codecontextobject on error resume next docmd.gotorecord , , acnewrec if err.number <> 0 btnaddrec.enabled = false end if end end sub everything ok when open window , click btnaddrec button, problem when first of perform navigation through existed records , after click on button. got runtime error: 2105: «you can't go specified record. may @ end of recordset» . how solve issue, need have ability add new record on click on specific button, no matter, have walked or not walked through records before. thanks. i created simple form field call description (and autonumber) , created button code follows click event. filled number of records , navigated through them, clicked addnewrec button. form navigated new record without issues. able click addnewrec button directly after
Read more

javascript - jQuery .height() return 0 when visible but non-0 when hidden -

i have problem here right opposite others: in jquery cannot height of element when it's visible, when hidden. have html code: <div id="filtering-filter-holder"> <div id="filtering-filter-holder-inner"></div> </div> basic css it: #filtering-filter-holder { display: none; position: absolute; width: 420px; min-height: 100px; padding: 5px 5px 5px 5px; background-color: #ffffff; border: 1px solid #164372; } and finally, js in document.ready(): $('div.filtering-filter-headline').mouseenter(function() { var el = $(this).parent().find('div.filtering-filter-content'); if ($('div#filtering-filter-holder').is(':visible')) { $('div#filtering-filter-holder-inner').stop(true, true).fadeout(300, function() { $('div#filtering-filter-holder-inner').html($(el).html()); var height = $('div#filtering-filter-holder-inner&#
Read more