00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041 #include <stdio.h>
00042 #ifndef __UNIX_
00043 #include <stdlib.h>
00044 #endif
00045 #include "mensagem.h"
00046 #include "memvirtu.h"
00047 #include "lowparam.h"
00048 #include "lowmacro.h"
00049 #include "lowsolid.h"
00050 #include "eulerops.h"
00051 #include "curvgeom.h"
00052 #include "mancommd.h"
00053
00054 #ifdef __Extended_Version
00055 void MSD_execKEF(void)
00056 {
00057 int sn, v1, v2, fn, op;
00058
00059 while (5 != sscanf(restbuf, "%d %d %d %d %d", &sn, &op, &v1, &v2, &fn))
00060 {
00061 printf("Solid V1 V2 F\n");
00062 if (!lineins("? "))
00063 {
00064 return;
00065 }
00066 }
00067 MSD_highKEF(sn, v1, v2, fn);
00068 }
00069
00070 #endif
00071
00072 int MSD_highKEF(Id sn, Id v1, Id v2, Id fn)
00073 {
00074 SPTYPE s;
00075 FPTYPE f1;
00076 HPTYPE he1;
00077
00078 if ((s = MSD_getSolid(sn)) == SNIL)
00079 {
00080 fprintf(stderr, MEN_NaoEncontrouSolido, NAM_KEF, sn);
00081 return(ERROR);
00082 }
00083 if ((f1 = MSD_getFace(s, fn)) == FNIL)
00084 {
00085 fprintf(stderr, MEN_NaoEncontrouFace, NAM_KEF, fn, sn);
00086 return(ERROR);
00087 }
00088 if (MSD_getEdge(f1, v1, v2, &he1) == LNIL)
00089 {
00090 if (MSD_getEdge(f1, v2, v1, &he1) == LNIL)
00091 {
00092 fprintf(stderr, MEN_NaoEncontrouArestaV, NAM_KEF, v1, v2, sn);
00093 return(ERROR);
00094 }
00095 }
00096 MSD_lowKEF(he1, mate(he1));
00097 return(SUCCESS);
00098 }
00099
00100 #ifdef __Extended_Version
00101 void MSD_execKEF2(void)
00102 {
00103 int sn, v1, v2, f1, f2, op;
00104
00105 while (6 != sscanf(restbuf, "%d %d %d %d %d %d", &sn, &op, &v1, &v2, &f1, &f2))
00106 {
00107 printf("kef2: Sid V1 V2 F1 F2\n");
00108 if (!lineins("? "))
00109 {
00110 return;
00111 }
00112 }
00113 MSD_highKEF2(sn, v1, v2, f1, f2);
00114 }
00115
00116 #endif
00117
00118 int MSD_highKEF2(Id sn, Id v1, Id v2, Id f1, Id f2)
00119 {
00120 SPTYPE s;
00121 FPTYPE ff1;
00122 FPTYPE ff2;
00123 HPTYPE he1;
00124
00125 if ((s = MSD_getSolid(sn)) == SNIL)
00126 {
00127 fprintf(stderr, MEN_NaoEncontrouSolido, NAM_KEF2, sn);
00128 return(ERROR);
00129 }
00130 if ((ff1 = MSD_getFace(s, f1)) == FNIL)
00131 {
00132 fprintf(stderr, MEN_NaoEncontrouFace, NAM_KEF2, f1, sn);
00133 return(ERROR);
00134 }
00135 if ((ff2 = MSD_getFace(s, f2)) == FNIL)
00136 {
00137 fprintf(stderr, MEN_NaoEncontrouFace, NAM_KEF2, f2, sn);
00138 return(ERROR);
00139 }
00140
00141 if (MSD_getHalfedge(ff1, ff2, v1, v2, &he1) == LNIL)
00142 {
00143 fprintf(stderr, MEN_NaoEncontrouArestaV, NAM_KEF, v1, v2, sn);
00144 return(ERROR);
00145 }
00146 MSD_lowKEF(he1, mate(he1));
00147 return(SUCCESS);
00148 }
00149
00150 void MSD_lowKEF(HPTYPE he1, HPTYPE he2)
00151 {
00152 SPTYPE s;
00153 DPTYPE d;
00154 FPTYPE f1;
00155 FPTYPE f2;
00156 LPTYPE l1;
00157 LPTYPE l2;
00158 LPTYPE l;
00159 EPTYPE e;
00160 HPTYPE he;
00161 HPTYPE he1nxt;
00162 HPTYPE he2nxt;
00163 HPTYPE hetmp;
00164 VPTYPE vp1;
00165 VPTYPE vp2;
00166 Id v1, v2, v3, v4, fa, fb, fc, fd, flag;
00167
00168 assert1(HalEdg(he1) != ENIL);
00169 assert1(HalEdg(he2) != ENIL);
00170 assert1(he1 == mate(he2));
00171 assert1((l1 = HalWLoop(he1)) != (l2 = HalWLoop(he2)));
00172 assert1((f1 = LooLFace(l1)) != (f2 = LooLFace(l2)));
00173 assert1(HalVtx(he2) == HalVtx(HalNxt(he1)));
00174 assert1(HalVtx(he1) == HalVtx(HalNxt(he2)));
00175
00176 #ifdef __curved_env
00177 if (he1->curv != (CUCURV *)NIL)
00178 {
00179 lsettag(he1, (CUCURV *)NIL);
00180 }
00181 if (he2->curv != (CUCURV *)NIL)
00182 {
00183 lsettag(he2, (CUCURV *)NIL);
00184 }
00185 if (f2->surf != (SUSURF *)NIL)
00186 {
00187 lsettag(f2, (SUSURF *)NIL);
00188 }
00189 #endif
00190
00191 d = FacFShell(f1);
00192 s = SheSSolids(d);
00193 e = HalEdg(he1);
00194
00195 if (Gdoinversion)
00196 {
00197 if ((HalNxt(he2) == he2) && (HalNxt(he1) != he1))
00198 {
00199 MSD_lowMEV(he1, HalNxt(he2), -1, 0.0, 0.0, 0.0);
00200 MSD_lowKEF(HalPrv(he1), mate(HalPrv(he1)));
00201 MSD_lowKEV(he1, mate(he1));
00202 return;
00203 }
00204 fa = FacFaceNo(f1);
00205 fd = FacFaceNo(f2);
00206 for (AllLoopsFace(f2, l))
00207 {
00208 if (l2 != l)
00209 {
00210 he = LooLEdg(l);
00211 v1 = VerVertexNo(HalVtx(he));
00212 v2 = HalNxt(he) != he ? VerVertexNo(HalVtx(HalNxt(he))) : 0;
00213 flag = l == FacFLOut(f2) ? TRUE : FALSE;
00214 MSD_lowAddEulerOp(s, RMOVE, v1, v2, fa, fd, flag, 0, 0, 0, 0.0, 0.0, 0.0, 0.0);
00215 }
00216 }
00217
00218 v1 = VerVertexNo(HalVtx(he1));
00219 v2 = VerVertexNo(HalVtx(he2));
00220 v3 = VerVertexNo(HalVtx(HalNxt(HalNxt(he2))));
00221 v4 = HalNxt(he1) == he1 ? v3 : VerVertexNo(HalVtx(HalNxt(HalNxt(he1))));
00222 fb = FacFaceNo(LooLFace(HalWLoop(mate(HalNxt(he2)))));
00223 fc = HalNxt(he1) == he1 ? fb : FacFaceNo(LooLFace(HalWLoop(mate(HalNxt(he1)))));
00224 if (fb == fd)
00225 {
00226 fb = fa;
00227 }
00228 if (fc == fd)
00229 {
00230 fc = fa;
00231 }
00232 MSD_lowAddEulerOp(s, MEF, v1, v2, v3, v4, fa, fb, fc, fd, 0.0, 0.0, 0.0, 0.0);
00233 }
00234 if (Gdodirect)
00235 {
00236 v1 = VerVertexNo(HalVtx(he1));
00237 v2 = VerVertexNo(HalVtx(he2));
00238 fa = FacFaceNo(LooLFace(HalWLoop(HalNxt(he1))));
00239 fb = FacFaceNo(LooLFace(HalWLoop(HalNxt(he2))));
00240 MSD_lowAddEulerOpRedo(s, KEF, v1, v2, fa, fb, 0, 0, 0, 0, 0.0, 0.0, 0.0, 0.0);
00241 }
00242
00243 he = LooLEdg(l2);
00244 do
00245 {
00246 HalWLoop(he) = l1;
00247 } while ((he = HalNxt(he)) != LooLEdg(l2));
00248 LooLength(l1) += LooLength(l2) - 2;
00249
00250 while ((l = FacFLoops(f2)) != LNIL)
00251 {
00252 MSD_lowDelList(LOOP, (NPTYPE)l, (NPTYPE)f2);
00253 MSD_lowAddList(LOOP, (NPTYPE)l, (NPTYPE)f1);
00254 }
00255
00256 he2nxt = HalNxt(he2) == he2 ? HalNxt(he1) : HalNxt(he2);
00257 he1nxt = HalNxt(he1) == he1 ? HalNxt(he2) : HalNxt(he1);
00258
00259 HalNxt(HalPrv(he1)) = he2;
00260 HalNxt(HalPrv(he2)) = he1;
00261 he = HalPrv(he2);
00262 hetmp = HalPrv(he1);
00263 HalPrv(he2) = hetmp;
00264 HalPrv(he1) = he;
00265
00266 vp1 = HalVtx(he1);
00267 vp2 = HalVtx(he2);
00268
00269 MSD_lowDelHE(he2);
00270 MSD_lowDelHE(he1);
00271
00272 VerVEdge(vp2) = he1nxt;
00273 if (HalEdg(VerVEdge(vp2)) == ENIL)
00274 {
00275 VerVEdge(vp2) = HNIL;
00276 }
00277 VerVEdge(vp1) = he2nxt;
00278 if (HalEdg(VerVEdge(vp1)) == ENIL)
00279 {
00280 VerVEdge(vp1) = HNIL;
00281 }
00282
00283 LooLEdg(l1) = he1nxt;
00284
00285 MSD_lowDelElement(LOOP, (NPTYPE)l2, (NPTYPE)f1);
00286 MSD_lowDelElement(FACE, (NPTYPE)f2, (NPTYPE)d);
00287 MSD_lowDelElement(EDGE, (NPTYPE)e, (NPTYPE)d);
00288
00289 MSD_addEulCoef(s, 0, -1, 0, -1, 0, 0);
00290
00291 if (SetTrace)
00292 {
00293 if (TraceSolid(s))
00294 {
00295 fprintf(stderr, MEN_MSDAbortou, NAM_KEF);
00296 exit(1);
00297 }
00298 }
00299 }
00300
00301 int MSD_execStructKEF(EulerOp *op)
00302 {
00303 Id sid, id1, id2, id3, id4;
00304
00305 sid = op->solidno;
00306 id1 = op->ip1;
00307 id2 = op->ip2;
00308 id3 = op->ip3;
00309 id4 = op->ip4;
00310 return(MSD_highKEF2(sid, id1, id2, id3, id4));
00311 }
00312
00313 void MSD_printStructKEF(FILE *f, EulerOp *op)
00314 {
00315 Id sid, id1, id2, id3, id4;
00316 int com;
00317
00318 sid = op->solidno;
00319 id1 = op->ip1;
00320 id2 = op->ip2;
00321 id3 = op->ip3;
00322 id4 = op->ip4;
00323 com = op->CommandCount;
00324 fprintf(f, "kef %d %d %d %d %d %d\n", sid, com, id1, id2, id3, id4);
00325 }
00326
00327 int MSD_readStructKEF(FILE *f, EulerOp *op)
00328 {
00329 int solidno, ip1, ip2, ip3, ip4, com;
00330
00331 if (fscanf(f, "%d %d %d %d %d %d\n", &solidno, &com, &ip1, &ip2, &ip3, &ip4) != 6)
00332 {
00333 return(ERROR);
00334 }
00335 op->solidno = solidno;
00336 op->CommandCount = com;
00337 op->ip1 = ip1;
00338 op->ip2 = ip2;
00339 op->ip3 = ip3;
00340 op->ip4 = ip4;
00341 op->ip5 = 0;
00342 op->ip6 = 0;
00343 op->ip7 = 0;
00344 op->ip8 = 0;
00345 op->fp1 = 0.0;
00346 op->fp2 = 0.0;
00347 op->fp3 = 0.0;
00348 op->fp4 = 0.0;
00349 return(SUCCESS);
00350 }
00351
00352 int MSD_modifyStructKEF(EulerOp *op)
00353 {
00354 int solidno;
00355
00356 solidno = MSD_lowCheckTranslatorTable(op->solidno, SOLID);
00357 op->solidno = solidno;
00358 return(SUCCESS);
00359 }