Lines Matching refs:sops
31 …ns(msym_equivalence_set_t *es, msym_thresholds_t *t, int *lsops, msym_symmetry_operation_t **sops);
42 msym_error_t reduceSymmetry(int sopsl, msym_symmetry_operation_t sops[sopsl], msym_thresholds_t *th…
43 msym_error_t filterSymmetryOperations(int sopsl, msym_symmetry_operation_t sops[sopsl], msym_thresh…
45 …sym_equivalence_set_t es[esl], msym_thresholds_t *t, int *lsops, msym_symmetry_operation_t **sops){ in findSymmetryOperations() argument
65 *sops = rsops; in findSymmetryOperations()
69 *sops = NULL; in findSymmetryOperations()
74 …ns(msym_equivalence_set_t *es, msym_thresholds_t *t, int *lsops, msym_symmetry_operation_t **sops){ in findEquivalenceSetSymmetryOperations() argument
117 if(*sops == NULL){ in findEquivalenceSetSymmetryOperations()
118 *sops = fsops; in findEquivalenceSetSymmetryOperations()
121 if(MSYM_SUCCESS != (ret = reduceSymmetry(lfsops, fsops, t, lsops, sops))) goto err; in findEquivalenceSetSymmetryOperations()
146 msym_symmetry_operation_t *sops = NULL; in findSymmetryLinear() local
161 sops = malloc(sopsl*sizeof(msym_symmetry_operation_t)); in findSymmetryLinear()
162 vcopy(v,sops[0].v); in findSymmetryLinear()
163 vcopy(v,sops[1].v); in findSymmetryLinear()
164 sops[0].type = PROPER_ROTATION; in findSymmetryLinear()
165 sops[0].order = 0; in findSymmetryLinear()
166 sops[0].power = 1; in findSymmetryLinear()
167 sops[1].type = REFLECTION; in findSymmetryLinear()
168 sops[1].order = 1; in findSymmetryLinear()
169 sops[1].power = 1; in findSymmetryLinear()
170 sops[2].type = INVERSION; in findSymmetryLinear()
171 sops[2].order = 1; in findSymmetryLinear()
172 sops[2].power = 1; in findSymmetryLinear()
173 sops[2].v[0] = sops[2].v[1] = sops[2].v[2] = 0; in findSymmetryLinear()
177 sops = malloc(sopsl*sizeof(msym_symmetry_operation_t)); in findSymmetryLinear()
178 vcopy(cm,sops[0].v); in findSymmetryLinear()
179 vnorm(sops[0].v); in findSymmetryLinear()
180 vcopy(ev[prim],sops[1].v); in findSymmetryLinear()
181 vnorm(sops[1].v); in findSymmetryLinear()
182 vcrossnorm(sops[0].v,sops[1].v,sops[2].v); in findSymmetryLinear()
183 sops[0].type = PROPER_ROTATION; in findSymmetryLinear()
184 sops[0].order = 2; in findSymmetryLinear()
185 sops[0].power = 1; in findSymmetryLinear()
186 sops[1].type = REFLECTION; in findSymmetryLinear()
187 sops[2].type = REFLECTION; in findSymmetryLinear()
190 *rsops = sops; in findSymmetryLinear()
204 msym_symmetry_operation_t *sops = NULL; in findSymmetryPlanarRegular() local
236 sops = malloc(sopsl*sizeof(msym_symmetry_operation_t)); in findSymmetryPlanarRegular()
239 sops[n].type = PROPER_ROTATION; in findSymmetryPlanarRegular()
240 sops[n].order = div[n]; in findSymmetryPlanarRegular()
241 sops[n].power = 1; in findSymmetryPlanarRegular()
242 vcopy(ev[prim],sops[n].v); in findSymmetryPlanarRegular()
245 sops[n].type = REFLECTION; in findSymmetryPlanarRegular()
246 vcopy(ev[prim],sops[n].v); in findSymmetryPlanarRegular()
250 sops[n].type = IMPROPER_ROTATION; in findSymmetryPlanarRegular()
251 sops[n].order = div[s]; in findSymmetryPlanarRegular()
252 sops[n].power = 1; in findSymmetryPlanarRegular()
253 vcopy(ev[prim],sops[n].v); in findSymmetryPlanarRegular()
260 sops[n].order = 1; in findSymmetryPlanarRegular()
261 sops[n].power = 1; in findSymmetryPlanarRegular()
262 sops[n].type = INVERSION; in findSymmetryPlanarRegular()
263 sops[n].v[0] = sops[n].v[1] = sops[n].v[2] = 0; in findSymmetryPlanarRegular()
272 vcrossnorm(r,ev[prim],sops[n].v); in findSymmetryPlanarRegular()
273 sops[n].type = REFLECTION; in findSymmetryPlanarRegular()
275 if(!findSymmetryOperation(&sops[n], sops, n, thresholds)){ in findSymmetryPlanarRegular()
278 vcopy(r,sops[n].v); in findSymmetryPlanarRegular()
279 sops[n].type = PROPER_ROTATION; in findSymmetryPlanarRegular()
280 sops[n].order = 2; in findSymmetryPlanarRegular()
281 sops[n].power = 1; in findSymmetryPlanarRegular()
294 *rsops = sops; in findSymmetryPlanarRegular()
300 free(sops); in findSymmetryPlanarRegular()
307 msym_symmetry_operation_t *sops = NULL; in findSymmetryPlanarIrregular() local
320 sops = malloc(sopsl*sizeof(msym_symmetry_operation_t)); in findSymmetryPlanarIrregular()
324 sops = malloc(sopsl*sizeof(msym_symmetry_operation_t)); in findSymmetryPlanarIrregular()
329 vcopy(ev[2],sops[0].v); in findSymmetryPlanarIrregular()
330 vnorm(sops[0].v); in findSymmetryPlanarIrregular()
331 sops[0].type = PROPER_ROTATION; in findSymmetryPlanarIrregular()
332 sops[0].order = 2; in findSymmetryPlanarIrregular()
333 sops[0].power = 1; in findSymmetryPlanarIrregular()
335 vcopy(ev[1],sops[1].v); in findSymmetryPlanarIrregular()
336 vnorm(sops[1].v); in findSymmetryPlanarIrregular()
337 sops[1].type = REFLECTION; in findSymmetryPlanarIrregular()
339 vcopy(ev[0],sops[2].v); in findSymmetryPlanarIrregular()
340 vnorm(sops[2].v); in findSymmetryPlanarIrregular()
341 sops[2].type = REFLECTION; in findSymmetryPlanarIrregular()
344 vcopy(sops[0].v, sops[3].v); in findSymmetryPlanarIrregular()
345 sops[3].type = REFLECTION; in findSymmetryPlanarIrregular()
347 vcopy(sops[1].v, sops[4].v); in findSymmetryPlanarIrregular()
348 sops[4].type = PROPER_ROTATION; in findSymmetryPlanarIrregular()
349 sops[4].order = 2; in findSymmetryPlanarIrregular()
350 sops[4].power = 1; in findSymmetryPlanarIrregular()
352 vcopy(sops[2].v, sops[5].v); in findSymmetryPlanarIrregular()
353 sops[5].type = PROPER_ROTATION; in findSymmetryPlanarIrregular()
354 sops[5].order = 2; in findSymmetryPlanarIrregular()
355 sops[5].power = 1; in findSymmetryPlanarIrregular()
357 sops[6].type = INVERSION; in findSymmetryPlanarIrregular()
358 sops[6].order = 1; in findSymmetryPlanarIrregular()
359 sops[6].power = 1; in findSymmetryPlanarIrregular()
360 sops[6].v[0] = sops[6].v[1] = sops[6].v[2] = 0; in findSymmetryPlanarIrregular()
365 *rsops = sops; in findSymmetryPlanarIrregular()
383 msym_symmetry_operation_t *sops = NULL; in findSymmetrySymmetricPolyhedron() local
470 sops = malloc(sopsl*sizeof(msym_symmetry_operation_t)); in findSymmetrySymmetricPolyhedron()
477 sops[n].type = PROPER_ROTATION; in findSymmetrySymmetricPolyhedron()
478 sops[n].order = div[n]; in findSymmetrySymmetricPolyhedron()
479 sops[n].power = 1; in findSymmetrySymmetricPolyhedron()
480 vcopy(ev[prim],sops[n].v); in findSymmetrySymmetricPolyhedron()
484 sops[n].type = REFLECTION; in findSymmetrySymmetricPolyhedron()
485 vcopy(ev[prim],sops[n].v); in findSymmetrySymmetricPolyhedron()
489 sops[n].type = IMPROPER_ROTATION; in findSymmetrySymmetricPolyhedron()
490 sops[n].order = div[s]; in findSymmetrySymmetricPolyhedron()
491 sops[n].power = 1; in findSymmetrySymmetricPolyhedron()
492 vcopy(ev[prim],sops[n].v); in findSymmetrySymmetricPolyhedron()
499 sops[n].type = INVERSION; in findSymmetrySymmetricPolyhedron()
500 sops[n].order = 1; in findSymmetrySymmetricPolyhedron()
501 sops[n].power = 1; in findSymmetrySymmetricPolyhedron()
502 sops[n].v[0] = sops[n].v[1] = sops[n].v[2] = 0; in findSymmetrySymmetricPolyhedron()
513 sops[n].type = IMPROPER_ROTATION; in findSymmetrySymmetricPolyhedron()
514 sops[n].order = 2*max; in findSymmetrySymmetricPolyhedron()
515 sops[n].power = 1; in findSymmetrySymmetricPolyhedron()
516 vcopy(ev[prim],sops[n].v); in findSymmetrySymmetricPolyhedron()
522 vrotate(theta_C2 + i*theta, v_init, ev[prim], sops[n].v); in findSymmetrySymmetricPolyhedron()
523 vnorm(sops[n].v); //Not really nessecary in findSymmetrySymmetricPolyhedron()
524 sops[n].type = PROPER_ROTATION; in findSymmetrySymmetricPolyhedron()
525 sops[n].order = 2; in findSymmetrySymmetricPolyhedron()
526 sops[n].power = 1; in findSymmetrySymmetricPolyhedron()
530 vrotate(i*theta, v_init,ev[prim], sops[n].v); in findSymmetrySymmetricPolyhedron()
531 vcrossnorm(sops[n].v,ev[prim],sops[n].v); in findSymmetrySymmetricPolyhedron()
532 sops[n].type = REFLECTION; in findSymmetrySymmetricPolyhedron()
546 *rsops = sops; in findSymmetrySymmetricPolyhedron()
551 free(sops); in findSymmetrySymmetricPolyhedron()
561 msym_symmetry_operation_t *sops = NULL; in findSymmetryAsymmetricPolyhedron() local
580 sops = malloc(sopsl*sizeof(msym_symmetry_operation_t)); in findSymmetryAsymmetricPolyhedron()
581 vcopy(ev[0], sops[0].v); in findSymmetryAsymmetricPolyhedron()
582 vcopy(ev[1], sops[1].v); in findSymmetryAsymmetricPolyhedron()
583 vcopy(ev[2], sops[2].v); in findSymmetryAsymmetricPolyhedron()
584 vnorm(sops[0].v); in findSymmetryAsymmetricPolyhedron()
585 vnorm(sops[1].v); in findSymmetryAsymmetricPolyhedron()
586 vnorm(sops[2].v); in findSymmetryAsymmetricPolyhedron()
587 sops[0].type = PROPER_ROTATION; in findSymmetryAsymmetricPolyhedron()
588 sops[0].order = 2; in findSymmetryAsymmetricPolyhedron()
589 sops[0].power = 1; in findSymmetryAsymmetricPolyhedron()
590 sops[1].type = PROPER_ROTATION; in findSymmetryAsymmetricPolyhedron()
591 sops[1].order = 2; in findSymmetryAsymmetricPolyhedron()
592 sops[1].power = 1; in findSymmetryAsymmetricPolyhedron()
593 sops[2].type = PROPER_ROTATION; in findSymmetryAsymmetricPolyhedron()
594 sops[2].order = 2; in findSymmetryAsymmetricPolyhedron()
595 sops[2].power = 1; in findSymmetryAsymmetricPolyhedron()
599 vcopy(sops[0].v,sops[3].v); in findSymmetryAsymmetricPolyhedron()
600 vcopy(sops[1].v,sops[4].v); in findSymmetryAsymmetricPolyhedron()
601 vcopy(sops[2].v,sops[5].v); in findSymmetryAsymmetricPolyhedron()
602 sops[3].type = REFLECTION; in findSymmetryAsymmetricPolyhedron()
603 sops[3].order = 1; in findSymmetryAsymmetricPolyhedron()
604 sops[3].power = 1; in findSymmetryAsymmetricPolyhedron()
605 sops[4].type = REFLECTION; in findSymmetryAsymmetricPolyhedron()
606 sops[4].order = 1; in findSymmetryAsymmetricPolyhedron()
607 sops[4].power = 1; in findSymmetryAsymmetricPolyhedron()
608 sops[5].type = REFLECTION; in findSymmetryAsymmetricPolyhedron()
609 sops[5].order = 1; in findSymmetryAsymmetricPolyhedron()
610 sops[5].power = 1; in findSymmetryAsymmetricPolyhedron()
611 sops[6].type = INVERSION; in findSymmetryAsymmetricPolyhedron()
612 sops[6].order = 1; in findSymmetryAsymmetricPolyhedron()
613 sops[6].power = 1; in findSymmetryAsymmetricPolyhedron()
614 sops[6].v[0] = sops[6].v[1] = sops[6].v[2] = 0; in findSymmetryAsymmetricPolyhedron()
619 *rsops = sops; in findSymmetryAsymmetricPolyhedron()
623 free(sops); in findSymmetryAsymmetricPolyhedron()
633 msym_symmetry_operation_t *sops = NULL; in findSymmetrySpherical() local
639 sops = malloc(sopsl*sizeof(msym_symmetry_operation_t)); in findSymmetrySpherical()
640 vcopy(t,sops[0].v); in findSymmetrySpherical()
641 vnorm(sops[0].v); in findSymmetrySpherical()
642 sops[0].type = PROPER_ROTATION; in findSymmetrySpherical()
643 sops[0].order = 0; in findSymmetrySpherical()
644 sops[0].power = 1; in findSymmetrySpherical()
648 *rsops = sops; in findSymmetrySpherical()
672 msym_symmetry_operation_t *sops = malloc(sizeof(msym_symmetry_operation_t[120])); in findSymmetryCubic() local
714 if(!findSymmetryOperation(&sopc, sops, sopsl, thresholds)){ in findSymmetryCubic()
717 copySymmetryOperation(&sops[sopsl], &sopc); in findSymmetryCubic()
718 (ac[4])[nc[4]++] = &sops[sopsl++]; in findSymmetryCubic()
722 if(!findSymmetryOperation(&sopc, sops, sopsl, thresholds)){ in findSymmetryCubic()
723 copySymmetryOperation(&sops[sopsl], &sopc); in findSymmetryCubic()
724 (ac[2])[nc[2]++] = &sops[sopsl++]; in findSymmetryCubic()
729 if(!findSymmetryOperation(&sopsigma, sops, sopsl, thresholds)){ in findSymmetryCubic()
732 copySymmetryOperation(&sops[sopsl], &sopsigma); in findSymmetryCubic()
733 sigma[nsigma++] = &sops[sopsl++]; in findSymmetryCubic()
743 if(!findSymmetryOperation(&sopc, sops, sopsl, thresholds)){ in findSymmetryCubic()
746 copySymmetryOperation(&sops[sopsl], &sopc); in findSymmetryCubic()
747 (ac[2])[nc[2]++] = &sops[sopsl++]; in findSymmetryCubic()
755 if(!findSymmetryOperation(&sopsigma, sops, sopsl, thresholds)){ in findSymmetryCubic()
759 copySymmetryOperation(&sops[sopsl], &sopsigma); in findSymmetryCubic()
760 sigma[nsigma++] = &sops[sopsl++]; in findSymmetryCubic()
770 for(msym_symmetry_operation_t *sopi = sops; sopi < (sops + sopsl);sopi++){ in findSymmetryCubic()
771 for(msym_symmetry_operation_t *sopj = sops; sopj < (sops + sopsl); sopj++){ in findSymmetryCubic()
774 copySymmetryOperation(&sops[sopsl], sopj); in findSymmetryCubic()
775 vcrossnorm(sopi->v, sopj->v, sops[sopsl].v); in findSymmetryCubic()
776 if(!findSymmetryOperation(&sops[sopsl], sops, sopsl,thresholds)){ in findSymmetryCubic()
777 (ac[2])[nc[2]++] = &sops[sopsl++]; in findSymmetryCubic()
781 copySymmetryOperation(&sops[sopsl], sopj); in findSymmetryCubic()
782 applySymmetryOperation(sopi,sops[sopsl].v,sops[sopsl].v); in findSymmetryCubic()
783 if(!findSymmetryOperation(&sops[sopsl], sops, sopsl,thresholds)){ in findSymmetryCubic()
785 sigma[nsigma++] = &sops[sopsl++]; in findSymmetryCubic()
787 (ac[sopj->order])[nc[sopj->order]++] = &sops[sopsl++]; in findSymmetryCubic()
817 copySymmetryOperation(&sops[sopsl++], &sopinversion); in findSymmetryCubic()
837 sops[sopsl].type = REFLECTION; in findSymmetryCubic()
838 vadd(sigma[i]->v, sigma[j]->v, sops[sopsl].v); in findSymmetryCubic()
839 vnorm(sops[sopsl].v); in findSymmetryCubic()
840 if(!findSymmetryOperation(&(sops[sopsl]), sops, sopsl,thresholds)){ in findSymmetryCubic()
841 sigma[nsigma+gsigma++] = &(sops[sopsl++]); in findSymmetryCubic()
846 sops[sopsl].type = REFLECTION; in findSymmetryCubic()
847 vsub(sigma[i]->v, sigma[j]->v, sops[sopsl].v); in findSymmetryCubic()
848 vnorm(sops[sopsl].v); in findSymmetryCubic()
849 if(!findSymmetryOperation(&(sops[sopsl]), sops, sopsl,thresholds)){ in findSymmetryCubic()
850 sigma[nsigma+gsigma++] = &(sops[sopsl++]); in findSymmetryCubic()
879 sops[sopsl].type = PROPER_ROTATION; in findSymmetryCubic()
880 sops[sopsl].order = 4; in findSymmetryCubic()
881 sops[sopsl].power = 1; in findSymmetryCubic()
882 vcopy((ac[2])[i]->v,sops[sopsl].v); in findSymmetryCubic()
883 if(!findSymmetryOperation(&sops[sopsl], sops, sopsl,thresholds)){ in findSymmetryCubic()
884 (ac[4])[nc[4]] = &(sops[sopsl++]); in findSymmetryCubic()
918 vadd(c3b[0]->v,c3b[1]->v,sops[sopsl].v); in findSymmetryCubic()
919 vadd(sops[sopsl].v,c3b[2]->v,sops[sopsl].v); in findSymmetryCubic()
920 vnorm(sops[sopsl].v); in findSymmetryCubic()
921 sops[sopsl].type = PROPER_ROTATION; in findSymmetryCubic()
922 sops[sopsl].order = 3; in findSymmetryCubic()
923 sops[sopsl].power = 1; in findSymmetryCubic()
924 (ac[3])[nc[3]] = &(sops[sopsl++]); in findSymmetryCubic()
928 vadd(c3b[0]->v,c3b[1]->v,sops[sopsl].v); in findSymmetryCubic()
929 vsub(sops[sopsl].v,c3b[2]->v,sops[sopsl].v); in findSymmetryCubic()
930 vnorm(sops[sopsl].v); in findSymmetryCubic()
931 sops[sopsl].type = PROPER_ROTATION; in findSymmetryCubic()
932 sops[sopsl].order = 3; in findSymmetryCubic()
933 sops[sopsl].power = 1; in findSymmetryCubic()
934 (ac[3])[nc[3]] = &(sops[sopsl++]); in findSymmetryCubic()
938 vadd(c3b[0]->v,c3b[2]->v,sops[sopsl].v); in findSymmetryCubic()
939 vsub(sops[sopsl].v,c3b[1]->v,sops[sopsl].v); in findSymmetryCubic()
940 vnorm(sops[sopsl].v); in findSymmetryCubic()
941 sops[sopsl].type = PROPER_ROTATION; in findSymmetryCubic()
942 sops[sopsl].order = 3; in findSymmetryCubic()
943 sops[sopsl].power = 1; in findSymmetryCubic()
944 (ac[3])[nc[3]] = &(sops[sopsl++]); in findSymmetryCubic()
948 vsub(c3b[0]->v,c3b[1]->v,sops[sopsl].v); in findSymmetryCubic()
949 vsub(sops[sopsl].v,c3b[2]->v,sops[sopsl].v); in findSymmetryCubic()
950 vnorm(sops[sopsl].v); in findSymmetryCubic()
951 sops[sopsl].type = PROPER_ROTATION; in findSymmetryCubic()
952 sops[sopsl].order = 3; in findSymmetryCubic()
953 sops[sopsl].power = 1; in findSymmetryCubic()
954 (ac[3])[nc[3]] = &(sops[sopsl++]); in findSymmetryCubic()
970 vcrossnorm(cb[i]->v, cb[j]->v, sops[sopsl].v); in findSymmetryCubic()
971 sops[sopsl].type = PROPER_ROTATION; in findSymmetryCubic()
972 sops[sopsl].order = k; in findSymmetryCubic()
973 sops[sopsl].power = 1; in findSymmetryCubic()
974 if(!findSymmetryOperation(&(sops[sopsl]), sops, sopsl,thresholds)){ in findSymmetryCubic()
975 (ac[k])[nc[k]] = &(sops[sopsl]); in findSymmetryCubic()
989 if(sops[i].type == PROPER_ROTATION){ in findSymmetryCubic()
991 if(!inversion && sops[i].order == 2){ in findSymmetryCubic()
992 vcopy(sops[i].v, sops[sopsl].v); in findSymmetryCubic()
993 sops[sopsl].type = IMPROPER_ROTATION; in findSymmetryCubic()
994 sops[sopsl].order = 4; in findSymmetryCubic()
995 sops[sopsl].power = 1; in findSymmetryCubic()
999 else if (inversion && sops[i].order > 2) { in findSymmetryCubic()
1000 vcopy(sops[i].v, sops[sopsl].v); in findSymmetryCubic()
1001 sops[sopsl].type = IMPROPER_ROTATION; in findSymmetryCubic()
1002 sops[sopsl].order = sops[i].order + (sops[i].order % 2)*sops[i].order; in findSymmetryCubic()
1003 sops[sopsl].power = 1; in findSymmetryCubic()
1018 *rsops = sops; in findSymmetryCubic()
1029 free(sops); in findSymmetryCubic()
1036 msym_error_t reduceSymmetry(int sopsl, msym_symmetry_operation_t sops[sopsl], msym_thresholds_t *th… in reduceSymmetry()
1053 …if(!((sops[i].type == PROPER_ROTATION && sops[i].order == 0) || sops[i].type == INVERSION || sops[… in reduceSymmetry()
1055 inv[1] = inv[1] || sops[i].type == INVERSION; in reduceSymmetry()
1056 if(sops[i].type == PROPER_ROTATION && sops[i].order == 0) cinf[1] = &sops[i]; in reduceSymmetry()
1141 …if(MSYM_SUCCESS != (ret = filterSymmetryOperations(sopsl,sops,thresholds,&rsopsl,&rsops))) goto er… in reduceSymmetry()
1153 if(sops[i].type == IMPROPER_ROTATION){ in reduceSymmetry()
1154 add = vparallel(sops[i].v,v,thresholds->angle); in reduceSymmetry()
1155 } else if(sops[i].type == PROPER_ROTATION){ in reduceSymmetry()
1156 if(sops[i].order != 2){ in reduceSymmetry()
1157 add = vparallel(sops[i].v,v,thresholds->angle); in reduceSymmetry()
1159 …add = vparallel(sops[i].v,v,thresholds->angle) || (vperpendicular(sops[i].v,cinf[0]->v,thresholds-… in reduceSymmetry()
1161 } else if(sops[i].type == REFLECTION){ in reduceSymmetry()
1162 add = vperpendicular(sops[i].v,v,thresholds->angle); in reduceSymmetry()
1166 copySymmetryOperation(&rsops[rsopsl], &sops[i]); in reduceSymmetry()
1170 …if(MSYM_SUCCESS != (ret = filterSymmetryOperations(sopsl,sops,thresholds,&rsopsl,&rsops))) goto er… in reduceSymmetry()
1173 … msym_symmetry_operation_t *fsop = findSymmetryOperation(&rsops[i], sops, sopsl,thresholds); in reduceSymmetry()
1202 …if(MSYM_SUCCESS != (ret = filterSymmetryOperations(sopsl,sops,thresholds, &rsopsl,&rsops))) goto e… in reduceSymmetry()
1212 msym_error_t filterSymmetryOperations(int sopsl, msym_symmetry_operation_t sops[sopsl], msym_thresh… in filterSymmetryOperations()
1218 msym_symmetry_operation_t *fsop = findSymmetryOperation(&rsops[i], sops, sopsl,thresholds); in filterSymmetryOperations()