Develop and Download Open Source Software

Powder auto-indexing software Conograph

Browse Subversion Repository

/[conograph]/Conograph/trunk/src/lattice_symmetry/LatticeFigureOfMeritToCheckSymmetry.cc

Diff of /Conograph/trunk/src/lattice_symmetry/LatticeFigureOfMeritToCheckSymmetry.cc

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 32 by rtomiyasu, Mon Jul 7 02:35:51 2014 UTC revision 33 by rtomiyasu, Wed Sep 7 04:38:51 2016 UTC
# Line 168  bool LatticeFigureOfMeritToCheckSymmetry Line 168  bool LatticeFigureOfMeritToCheckSymmetry
168                  {                  {
169                          S_red = ans0;                          S_red = ans0;
170                          trans_mat2 = m_S_red.second;                          trans_mat2 = m_S_red.second;
171                          putBuergerReducedMonoclinicP(1, 2, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(1, 2, S_red, trans_mat2);
172                  }                  }
173                  else if( epg_new == C2h_Y )                  else if( epg_new == C2h_Y )
174                  {                  {
175                          S_red = transform_sym_matrix(put_matrix_YXZ(), ans0);                          S_red = transform_sym_matrix(put_matrix_YXZ(), ans0);
176                          trans_mat2 = mprod(m_S_red.second, put_matrix_YXZ());                          trans_mat2 = mprod(m_S_red.second, put_matrix_YXZ());
177                          putBuergerReducedMonoclinicP(0, 2, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(0, 2, S_red, trans_mat2);
178                  }                  }
179                  else // if( epg_new == C2h_Z )                  else // if( epg_new == C2h_Z )
180                  {                  {
181                            // (0 0 1)(0 1 0)
182                            // (1 0 0)(0 0 1)
183                            // (0 1 0)(1 0 0)
184                          S_red = transform_sym_matrix(put_matrix_YZX(), ans0);                          S_red = transform_sym_matrix(put_matrix_YZX(), ans0);
185                          trans_mat2 = mprod(m_S_red.second, put_matrix_ZXY());                          trans_mat2 = mprod(m_S_red.second, put_matrix_ZXY());
186                          putBuergerReducedMonoclinicP(0, 1, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(0, 1, S_red, trans_mat2);
187                  }                  }
188                  ans.insert( SymMat43_Double( S_red, trans_mat2) );                  ans.insert( SymMat43_Double( S_red, trans_mat2) );
189          }          }
# Line 193  bool LatticeFigureOfMeritToCheckSymmetry Line 196  bool LatticeFigureOfMeritToCheckSymmetry
196                  {                  {
197                          S_red = transform_sym_matrix(put_matrix_YXZ(), ans0);                          S_red = transform_sym_matrix(put_matrix_YXZ(), ans0);
198                          trans_mat2 = mprod(m_S_red.second, put_matrix_YXZ());                          trans_mat2 = mprod(m_S_red.second, put_matrix_YXZ());
199                          putBuergerReducedMonoclinicP(1, 2, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(1, 2, S_red, trans_mat2);
200                  }                  }
201                  else if( epg_new == C2h_Y )                  else if( epg_new == C2h_Y )
202                  {                  {
203                          S_red = ans0;                          S_red = ans0;
204                          trans_mat2 = m_S_red.second;                          trans_mat2 = m_S_red.second;
205                          putBuergerReducedMonoclinicP(0, 2, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(0, 2, S_red, trans_mat2);
206                  }                  }
207                  else // if( epg_new == C2h_Z )                  else // if( epg_new == C2h_Z )
208                  {                  {
209                          S_red = transform_sym_matrix(put_matrix_XZY(), ans0);                          S_red = transform_sym_matrix(put_matrix_XZY(), ans0);
210                          trans_mat2 = mprod(m_S_red.second, put_matrix_XZY());                          trans_mat2 = mprod(m_S_red.second, put_matrix_XZY());
211                          putBuergerReducedMonoclinicP(0, 1, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(0, 1, S_red, trans_mat2);
212                  }                  }
213                  ans.insert( SymMat43_Double( S_red, trans_mat2) );                  ans.insert( SymMat43_Double( S_red, trans_mat2) );
214          }          }
# Line 216  bool LatticeFigureOfMeritToCheckSymmetry Line 219  bool LatticeFigureOfMeritToCheckSymmetry
219          {          {
220                  if( epg_new == C2h_X )                  if( epg_new == C2h_X )
221                  {                  {
222                            // (0 1 0)(0 0 1)
223                            // (0 0 1)(1 0 0)
224                            // (1 0 0)(0 1 0)
225                          S_red = transform_sym_matrix(put_matrix_ZXY(), ans0);                          S_red = transform_sym_matrix(put_matrix_ZXY(), ans0);
226                          trans_mat2 = mprod(m_S_red.second, put_matrix_YZX());                          trans_mat2 = mprod(m_S_red.second, put_matrix_YZX());
227                          putBuergerReducedMonoclinicP(1, 2, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(1, 2, S_red, trans_mat2);
228                  }                  }
229                  else if( epg_new == C2h_Y )                  else if( epg_new == C2h_Y )
230                  {                  {
231                          S_red = transform_sym_matrix(put_matrix_XZY(), ans0);                          S_red = transform_sym_matrix(put_matrix_XZY(), ans0);
232                          trans_mat2 = mprod(m_S_red.second, put_matrix_XZY());                          trans_mat2 = mprod(m_S_red.second, put_matrix_XZY());
233                          putBuergerReducedMonoclinicP(0, 2, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(0, 2, S_red, trans_mat2);
234                  }                  }
235                  else // if( epg_new == C2h_Z )                  else // if( epg_new == C2h_Z )
236                  {                  {
237                          S_red = ans0;                          S_red = ans0;
238                          trans_mat2 = m_S_red.second;                          trans_mat2 = m_S_red.second;
239                          putBuergerReducedMonoclinicP(0, 1, S_red, trans_mat2);                          putBuergerReducedMonoclinicP<Double>(0, 1, S_red, trans_mat2);
240                  }                  }
241                  ans.insert( SymMat43_Double( S_red, trans_mat2) );                  ans.insert( SymMat43_Double( S_red, trans_mat2) );
242          }          }
Legend:
Removed from v.32  
changed lines
  Added in v.33
Back to OSDN">Back to OSDN
 ViewVC Help
Powered by ViewVC 1.1.26  
About OSDN
Find Software
Develop Software
Help