Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
wppm_macros [2020/03/04 07:21]
rowlesmr3 doubled up the comment command to help with the automatic syntax colouring
wppm_macros [2020/07/16 11:29] (current)
Line 11: Line 11:
 <code topas> <code topas>
  
-'' WPPM macros by P. Scardi, C. Perez-Demydenko & C.L. Azanza Ricardo +' WPPM macros by P. Scardi, C. Perez-Demydenko & C.L. Azanza Ricardo 
-'' See also: P. Scardi, C. Perez-Demydenko,​ C.L. Azanza Ricardo & A. Coelho, Journal of Applied Crystallography,​ 2018 +' See also: P. Scardi, C. Perez-Demydenko,​ C.L. Azanza Ricardo & A. Coelho, Journal of Applied Crystallography,​ 2018 
-'' The following macros can be used with any space group; additional macros specific to cubic phases are available +' The following macros can be used with any space group; additional macros specific to cubic phases are available 
-'' email: Paolo.Scardi@unitn.it+' email: Paolo.Scardi@unitn.it
  
 macro WPPM_Sphere(RRc,​ RRv) macro WPPM_Sphere(RRc,​ RRv)
Line 21: Line 21:
       WPPM_ft_conv = 1 - 0.75*WPPM_L/​RRv + 0.0625*(WPPM_L/​RRv)^3 ;       WPPM_ft_conv = 1 - 0.75*WPPM_L/​RRv + 0.0625*(WPPM_L/​RRv)^3 ;
  
-      ​''​WPPM_break_on_small = 1e-7;+      '​WPPM_break_on_small = 1e-7;
       WPPM_L_max = 2*CeV(RRc, RRv);       WPPM_L_max = 2*CeV(RRc, RRv);
       WPPM_th2_range = 55;       WPPM_th2_range = 55;
Line 51: Line 51:
 macro getInvariant(sg,​ E1, E2, E3, E4, E5, E6, E7, E8, E9, E10, E11, E12, E13, E14, E15) macro getInvariant(sg,​ E1, E2, E3, E4, E5, E6, E7, E8, E9, E10, E11, E12, E13, E14, E15)
    {    {
-      #if Or(sg == 1, sg == 2); ''​Triclinic,​ -1+      #if Or(sg == 1, sg == 2); '​Triclinic,​ -1
          ​E1*H^4 + E2*K^4 + E3*L^4 + 2*(E4*H*H*K*K + E5*K*K*L*L + E6*H*H*L*L) + 4*(E7*(H^3)*K + E8*(H^3)*L + E9*H*(K^3) + E10*(K^3)*L + E11*H*(L^3) + E12*K*(L^3) + E13*H*H*K*L + E14*H*K*K*L + E15*H*K*L*L)          ​E1*H^4 + E2*K^4 + E3*L^4 + 2*(E4*H*H*K*K + E5*K*K*L*L + E6*H*H*L*L) + 4*(E7*(H^3)*K + E8*(H^3)*L + E9*H*(K^3) + E10*(K^3)*L + E11*H*(L^3) + E12*K*(L^3) + E13*H*H*K*L + E14*H*K*K*L + E15*H*K*L*L)
       ​       ​
-      #elseif And(sg >= 3, sg <= 15); ''​Monoclinic,​ 2/m +      #elseif And(sg >= 3, sg <= 15); '​Monoclinic,​ 2/m 
-         If (And(Constant(Get(al)) == 90, Constant(Get(ga)) == 90), '' unique axis b +         If (And(Constant(Get(al)) == 90, Constant(Get(ga)) == 90), ' unique axis b 
             E1*H^4 + E2*K^4 + E3*L^4 + 2*(E4*H*H*L*L + E5*K*K*L*L + E6*H*H*K*K) + 4*(E7*(H^3)*L + E8*H*(L^3) + E9*H*K*K*L),​             E1*H^4 + E2*K^4 + E3*L^4 + 2*(E4*H*H*L*L + E5*K*K*L*L + E6*H*H*K*K) + 4*(E7*(H^3)*L + E8*H*(L^3) + E9*H*K*K*L),​
-            If (And(Constant(Get(al)) == 90, Constant(Get(be)) == 90), '' unique axis c +            If (And(Constant(Get(al)) == 90, Constant(Get(be)) == 90), ' unique axis c 
                ​E1*H^4 + E2*K^4 + E3*L^4 + 2*(E4*H*H*K*K + E5*K*K*L*L + E6*H*H*L*L) + 4*(E7*(H^3)*K + E8*H*(K^3) + E9*H*K*L*L),​                ​E1*H^4 + E2*K^4 + E3*L^4 + 2*(E4*H*H*K*K + E5*K*K*L*L + E6*H*H*L*L) + 4*(E7*(H^3)*K + E8*H*(K^3) + E9*H*K*L*L),​
                "​Invalid Monoclinic lattice parameters"​                "​Invalid Monoclinic lattice parameters"​
             )             )
          )          )
-      #elseif And(sg >= 16, sg <= 74); ''​Orthorhombic,​ mmm+      #elseif And(sg >= 16, sg <= 74); '​Orthorhombic,​ mmm
          ​E1*H^4 + E2*K^4 + E3*L^4 + 2*(E4*H*H*K*K + E5*K*K*L*L + E6*H*H*L*L)          ​E1*H^4 + E2*K^4 + E3*L^4 + 2*(E4*H*H*K*K + E5*K*K*L*L + E6*H*H*L*L)
  
-      #elseif And(sg >= 75, sg <= 88); ''​Tetragonal,​ 4/m+      #elseif And(sg >= 75, sg <= 88); '​Tetragonal,​ 4/m
          ​E1*(H^4 + K^4) + E2*L^4 + 2*E3*H*H*K*K + 2*E4*(H*H + K*K)*L*L + 4*E5*H*K*(H*H - K*K)          ​E1*(H^4 + K^4) + E2*L^4 + 2*E3*H*H*K*K + 2*E4*(H*H + K*K)*L*L + 4*E5*H*K*(H*H - K*K)
                
-      #elseif And(sg >= 89, sg <= 142); ''​Tetragonal,​ 4/mmm+      #elseif And(sg >= 89, sg <= 142); '​Tetragonal,​ 4/mmm
          ​E1*(H^4 + K^4) + E2*L^4 + 2*E3*H*H*K*K + 2*E4*(H*H + K*K)*L*L          ​E1*(H^4 + K^4) + E2*L^4 + 2*E3*H*H*K*K + 2*E4*(H*H + K*K)*L*L
  
-      #elseif And(sg >= 143, sg <= 148); ''​Trigonal,​ -3 +      #elseif And(sg >= 143, sg <= 148); '​Trigonal,​ -3 
-         #if Or(sg == 146, sg == 148); ''​Trigonal space groups supporting also rhombohedral axes      ​+         #if Or(sg == 146, sg == 148); '​Trigonal space groups supporting also rhombohedral axes      ​
             If (Constant(Get(ga)) == 120,             If (Constant(Get(ga)) == 120,
                ​E1*(H*H - H*K + K*K)^2 + 2*E2*(H*H - H*K + K*K)*L*L + E3*L^4 + 4*E5*(H^3 - 3*H*K*K + (K^3))*L + 4*E4*H*(H - K)*K*L,                ​E1*(H*H - H*K + K*K)^2 + 2*E2*(H*H - H*K + K*K)*L*L + E3*L^4 + 4*E5*(H^3 - 3*H*K*K + (K^3))*L + 4*E4*H*(H - K)*K*L,
                If (Constant(Get(ga)) == 60,                If (Constant(Get(ga)) == 60,
-                  E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + (4/​3)*E4*L*(H^3 + 3*H*H*K - K^3) + (4/​3)*E5*L*(-H^3 + 3*H*K*K + K^3), '' 60deg+                  E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + (4/​3)*E4*L*(H^3 + 3*H*H*K - K^3) + (4/​3)*E5*L*(-H^3 + 3*H*K*K + K^3), ' 60deg
                   E1*(H^4 + K^4 + L^4) + 2*E2*(H*H*K*K + H*H*L*L + K*K*L*L) + 4*E3*((H^3)*K + (K^3)*L + H*(L^3)) + 4*E4*(H*(K^3) + (H^3)*L + K*(L^3))+ 4*E5*H*K*L*(H + K + L)                   E1*(H^4 + K^4 + L^4) + 2*E2*(H*H*K*K + H*H*L*L + K*K*L*L) + 4*E3*((H^3)*K + (K^3)*L + H*(L^3)) + 4*E4*(H*(K^3) + (H^3)*L + K*(L^3))+ 4*E5*H*K*L*(H + K + L)
                )                )
Line 83: Line 83:
             If (Constant(Get(ga)) == 120,             If (Constant(Get(ga)) == 120,
                ​E1*(H*H - H*K + K*K)^2 + 2*E2*(H*H - H*K + K*K)*L*L + E3*L^4 + 4*E5*(H^3 - 3*H*K*K + (K^3))*L + 4*E4*H*(H - K)*K*L,                ​E1*(H*H - H*K + K*K)^2 + 2*E2*(H*H - H*K + K*K)*L*L + E3*L^4 + 4*E5*(H^3 - 3*H*K*K + (K^3))*L + 4*E4*H*(H - K)*K*L,
-               ​E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + (4/​3)*E4*L*(H^3 + 3*H*H*K - K^3) + (4/​3)*E5*L*(-H^3 + 3*H*K*K + K^3) '' 60+               ​E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + (4/​3)*E4*L*(H^3 + 3*H*H*K - K^3) + (4/​3)*E5*L*(-H^3 + 3*H*K*K + K^3) ' 60
             )             )
          #​endif          #​endif
       ​       ​
-      #elseif And(sg >= 149, sg <= 167); ''​Trigonal,​ -3m +      #elseif And(sg >= 149, sg <= 167); '​Trigonal,​ -3m 
-         #if Or(sg ==150, sg ==152, sg ==154, sg ==155, sg ==156, sg ==158, sg ==160, sg ==161, sg ==164, sg ==165 , sg ==166, sg ==167); ​''​Trigonal,​ -3m1 +         #if Or(sg ==150, sg ==152, sg ==154, sg ==155, sg ==156, sg ==158, sg ==160, sg ==161, sg ==164, sg ==165 , sg ==166, sg ==167); '​Trigonal,​ -3m1 
-            #if Or(sg == 155, sg == 160, sg == 161, sg == 166, sg == 167); ''​Trigonal space groups supporting also rhombohedral axes      ​+            #if Or(sg == 155, sg == 160, sg == 161, sg == 166, sg == 167); '​Trigonal space groups supporting also rhombohedral axes      ​
                If (Constant(Get(ga)) == 120,                If (Constant(Get(ga)) == 120,
                   E1*(H*H - H*K + K*K)^2 + 2*E2*L*L*(H*H - H*K + K*K) + E3*L^4 + 4*E4*H*K*L*(H - K),                   E1*(H*H - H*K + K*K)^2 + 2*E2*L*L*(H*H - H*K + K*K) + E3*L^4 + 4*E4*H*K*L*(H - K),
                   If (Constant(Get(ga)) == 60,                   If (Constant(Get(ga)) == 60,
-                     ​E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + 4*E4*H*K*L*(H + K), '' 60deg+                     ​E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + 4*E4*H*K*L*(H + K), ' 60deg
                      ​E1*(H^4 + K^4 + L^4) + 2*E2*(H*H*K*K + H*H*L*L + K*K*L*L) + 4*E3*(H*K*(H*H + K*K) + H*L*(H*H + L*L) + K*L*(K*K + L*L)) + 4*E4*H*K*L*(H + K + L)                      ​E1*(H^4 + K^4 + L^4) + 2*E2*(H*H*K*K + H*H*L*L + K*K*L*L) + 4*E3*(H*K*(H*H + K*K) + H*L*(H*H + L*L) + K*L*(K*K + L*L)) + 4*E4*H*K*L*(H + K + L)
                   )                   )
Line 100: Line 100:
                If (Constant(Get(ga)) == 120,                If (Constant(Get(ga)) == 120,
                   E1*(H*H - H*K + K*K)^2 + 2*E2*L*L*(H*H - H*K + K*K) + E3*L^4 + 4*E4*H*K*L*(H - K),                   E1*(H*H - H*K + K*K)^2 + 2*E2*L*L*(H*H - H*K + K*K) + E3*L^4 + 4*E4*H*K*L*(H - K),
-                  E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + 4*E4*H*K*L*(H + K) '' 60deg+                  E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + 4*E4*H*K*L*(H + K) ' 60deg
                )                )
             #​endif ​           ​             #​endif ​           ​
-         #​else ​''​Trigonal,​ -31m+         #else '​Trigonal,​ -31m
             If (Constant(Get(ga)) == 120,             If (Constant(Get(ga)) == 120,
                ​E1*(H*H - H*K + K*K)^2 + 2*E2*L*L*(H*H - H*K + K*K) + E3*L^4 + E4*L*(4*(H^3) - 6*H*H*K - 6*H*K*K + 4*(K^3)),                ​E1*(H*H - H*K + K*K)^2 + 2*E2*L*L*(H*H - H*K + K*K) + E3*L^4 + E4*L*(4*(H^3) - 6*H*H*K - 6*H*K*K + 4*(K^3)),
-               ​E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + (4/​3)*E4*L*(2*(H^3) + 3*H*H*K - 3*H*K*K - 2*(K^3)) ​'' 60deg+               ​E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 + (4/​3)*E4*L*(2*(H^3) + 3*H*H*K - 3*H*K*K - 2*(K^3)) ' 60deg
             )             )
          #​endif          #​endif
       ​       ​
-      #elseif And(sg >= 168, sg <= 194); ''​Hexagonal,​ 6/m and 6/mmm+      #elseif And(sg >= 168, sg <= 194); '​Hexagonal,​ 6/m and 6/mmm
          If (Constant(Get(ga)) == 120,          If (Constant(Get(ga)) == 120,
             E1*(H*H - H*K + K*K)^2 + 2*E2*L*L*(H*H - H*K + K*K) + E3*L^4,             E1*(H*H - H*K + K*K)^2 + 2*E2*L*L*(H*H - H*K + K*K) + E3*L^4,
-            E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 ​'' 60deg+            E1*(H*H + H*K + K*K)^2 + 2*E2*L*L*(H*H + H*K + K*K) + E3*L^4 ' 60deg
          )          )
-      #elseif And(sg >= 195, sg <= 230); ''​Cubic,​ m-3 and m-3m+      #elseif And(sg >= 195, sg <= 230); '​Cubic,​ m-3 and m-3m
          ​E1*(H^4 + K^4 + L^4) + 2*E2*(H*H*K*K + H*H*L*L + K*K*L*L)          ​E1*(H^4 + K^4 + L^4) + 2*E2*(H*H*K*K + H*H*L*L + K*K*L*L)
       #else       #else
Line 122: Line 122:
 macro Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5, E6, E7, E8, E9)    ​ macro Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5, E6, E7, E8, E9)    ​
    {    {
-      ​'' E1 to E9 for space groups 3 to 15 (2/m)+      ' E1 to E9 for space groups 3 to 15 (2/m)
       Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5, E6, E7, E8, E9,,,,,,​) ​       Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5, E6, E7, E8, E9,,,,,,​) ​
    }    }
 macro Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5, E6)  macro Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5, E6) 
    {    {
-      ​'' E1 to E6 for space groups 16 to 64 (mmm)+      ' E1 to E6 for space groups 16 to 64 (mmm)
       Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5, E6,,,,,,,,,​) ​       Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5, E6,,,,,,,,,​) ​
    }    }
Line 133: Line 133:
 macro Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5)  macro Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5) 
    {    {
-      ​'' E1 to E5 for space groups 75 to 88 (4/m) and 143 to 148 (-3)+      ' E1 to E5 for space groups 75 to 88 (4/m) and 143 to 148 (-3)
       Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5,,,,,,,,,,​) ​       Out_Invariant(file,​sg,​ E1, E2, E3, E4, E5,,,,,,,,,,​) ​
    }    }
 macro Out_Invariant(file,​sg,​ E1, E2, E3, E4)  macro Out_Invariant(file,​sg,​ E1, E2, E3, E4) 
    {    {
-      ​'' E1 to E4 for space groups 89 to 142 (4/mmm) and 149 to 167 (-3m)+      ' E1 to E4 for space groups 89 to 142 (4/mmm) and 149 to 167 (-3m)
       Out_Invariant(file,​sg,​ E1, E2, E3, E4,,,,,,,,,,,​) ​       Out_Invariant(file,​sg,​ E1, E2, E3, E4,,,,,,,,,,,​) ​
    }    }
 macro Out_Invariant(file,​sg,​ E1, E2, E3)  macro Out_Invariant(file,​sg,​ E1, E2, E3) 
    {    {
-      ​'' E1 to E3 for space groups 168 to 194 (6/m and 6/mmm)+      ' E1 to E3 for space groups 168 to 194 (6/m and 6/mmm)
       Out_Invariant(file,​sg,​ E1, E2, E3,,,,,,,,,,,,​) ​       Out_Invariant(file,​sg,​ E1, E2, E3,,,,,,,,,,,,​) ​
    }    }
 macro Out_Invariant(file,​sg,​ E1, E2)  macro Out_Invariant(file,​sg,​ E1, E2) 
    {    {
-      ​'' E1 to E2 for space groups 195 to 230 (m-3 and m-3m)+      ' E1 to E2 for space groups 195 to 230 (m-3 and m-3m)
       Out_Invariant(file,​sg,​ E1, E2,,,,,,,,,,,,,​) ​       Out_Invariant(file,​sg,​ E1, E2,,,,,,,,,,,,,​) ​
    }    }
Line 179: Line 179:
 macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5, E6, E7, E8, E9)  macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5, E6, E7, E8, E9) 
    {    {
-      ​'' E1 to E9 for space groups 3 to 15 (2/m)+      ' E1 to E9 for space groups 3 to 15 (2/m)
       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5, E6, E7, E8, E9,,,,,,​) ​       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5, E6, E7, E8, E9,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5, E6)  macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5, E6) 
    {    {
-      ​'' E1 to E6 for space groups 16 to 64 (mmm)+      ' E1 to E6 for space groups 16 to 64 (mmm)
       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5, E6,,,,,,,,,​) ​       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5, E6,,,,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5)  macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5) 
    {    {
-      ​'' E1 to E5 for space groups 75 to 88 (4/m) and 143 to 148 (-3)+      ' E1 to E5 for space groups 75 to 88 (4/m) and 143 to 148 (-3)
       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5,,,,,,,,,,​) ​       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4, E5,,,,,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4)  macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4) 
    {    {
-      ​'' E1 to E4 for space groups 89 to 142 (4/mmm) and 149 to 167 (-3m)+      ' E1 to E4 for space groups 89 to 142 (4/mmm) and 149 to 167 (-3m)
       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4,,,,,,,,,,,​) ​       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3, E4,,,,,,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3)  macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3) 
    {    {
-      ​'' E1 to E3 for space groups 168 to 194 (6/m and 6/mmm)+      ' E1 to E3 for space groups 168 to 194 (6/m and 6/mmm)
       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3,,,,,,,,,,,,​) ​       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2, E3,,,,,,,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2)  macro WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2) 
    {    {
-      ​'' E1 to E2 for space groups 195 to 230 (m-3 and m-3m)+      ' E1 to E2 for space groups 195 to 230 (m-3 and m-3m)
       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2,,,,,,,,,,,,,​) ​       WPPM_Strain_InvariantWilkens(rhoc,​ rhov, Rec, Rev, burger, wsg, E1, E2,,,,,,,,,,,,,​) ​
    }    }
Line 234: Line 234:
 macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5, E6, E7, E8, E9)  macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5, E6, E7, E8, E9) 
    {    {
-      ​'' E1 to E9 for space groups 3 to 15 (2/m)+      ' E1 to E9 for space groups 3 to 15 (2/m)
       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5, E6, E7, E8, E9,,,,,,​) ​       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5, E6, E7, E8, E9,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5, E6)  macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5, E6) 
    {    {
-      ​'' E1 to E6 for space groups 16 to 64 (mmm)+      ' E1 to E6 for space groups 16 to 64 (mmm)
       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5, E6,,,,,,,,,​) ​       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5, E6,,,,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5)  macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5) 
    {    {
-      ​'' E1 to E5 for space groups 75 to 88 (4/m) and 143 to 148 (-3)+      ' E1 to E5 for space groups 75 to 88 (4/m) and 143 to 148 (-3)
       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5,,,,,,,,,,​) ​       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4, E5,,,,,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4)  macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4) 
    {    {
-      ​'' E1 to E4 for space groups 89 to 142 (4/mmm) and 149 to 167 (-3m)+      ' E1 to E4 for space groups 89 to 142 (4/mmm) and 149 to 167 (-3m)
       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4,,,,,,,,,,,​) ​       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3, E4,,,,,,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3)  macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3) 
    {    {
-      ​'' E1 to E3 for space groups 168 to 194 (6/m and 6/mmm)+      ' E1 to E3 for space groups 168 to 194 (6/m and 6/mmm)
       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3,,,,,,,,,,,,​) ​       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2, E3,,,,,,,,,,,,​) ​
    }    }
 macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2)  macro WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2) 
    {    {
-      ​'' E1 to E2 for space groups 195 to 230 (m-3 and m-3m)+      ' E1 to E2 for space groups 195 to 230 (m-3 and m-3m)
       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2,,,,,,,,,,,,,​) ​       WPPM_Strain_InvariantPAH(ac,​ av, bc, bv, wsg, E1, E2,,,,,,,,,,,,,​) ​
    }    }

Personal Tools