Changeset 2839


Ignore:
Timestamp:
May 24, 2017 3:10:28 PM (5 years ago)
Author:
vondreele
Message:

fix numpy divide issues inquaternions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIImath.py

    r2838 r2839  
    44744474    d = nl.norm(np.array(V))
    44754475    if d:
    4476         V /= d
     4476        V = V/d
    44774477        if not A:       #==0.
    44784478            A = 2.*np.pi
     
    44934493        A = 360.
    44944494    if d:
    4495         V /= d
     4495        V = V/d
    44964496        p = A/2.
    44974497        Q[0] = cosd(p)
     
    45074507    A = 2.*acosd(Q[0])
    45084508    V = np.array(Q[1:])
    4509     V /= sind(A/2.)
     4509    V = V/sind(A/2.)
    45104510    return A,V
    45114511   
     
    45164516    A = 2.*np.arccos(Q[0])
    45174517    V = np.array(Q[1:])
    4518     V /= np.sin(A/2.)
     4518    V = V/np.sin(A/2.)
    45194519    return A,V
    45204520   
     
    45484548    V2 = np.cross(B,C)
    45494549    if nl.norm(V1)*nl.norm(V2):
    4550         V1 /= nl.norm(V1)
    4551         V2 /= nl.norm(V2)
     4550        V1 = V1/nl.norm(V1)
     4551        V2 = V2/nl.norm(V2)
    45524552        V3 = np.cross(V1,V2)
    45534553    else:
     
    45564556    D = 0.
    45574557    if nl.norm(V3):
    4558         V3 /= nl.norm(V3)
     4558        V3 = V3/nl.norm(V3)
    45594559        D1 = min(1.0,max(-1.0,np.vdot(V1,V2)))
    45604560        D = np.arccos(D1)/2.0
Note: See TracChangeset for help on using the changeset viewer.