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_execMEF(void)
00056 {
00057 int sn, v1, v2, f1, f2, op;
00058
00059 while (6 != sscanf(restbuf, "%d %d %d %d %d %d", &sn, &op, &v1, &v2, &f1, &f2))
00060 {
00061 printf("mef: Sid V1 V2 F1 F2\n");
00062 if (!lineins("? "))
00063 {
00064 return;
00065 }
00066 }
00067 MSD_highMEF(sn, v1, v2, f1, f2);
00068 }
00069
00070 #endif
00071
00072 int MSD_highMEF(Id sn, Id v1, Id v2, Id f1, Id f2)
00073 {
00074 SPTYPE s;
00075 FPTYPE oldface;
00076 LPTYPE oldloop;
00077 HPTYPE he1;
00078 HPTYPE he2;
00079
00080 if ((s = MSD_getSolid(sn)) == SNIL)
00081 {
00082 fprintf(stderr, MEN_NaoEncontrouSolido, NAM_MEF, sn);
00083 return(ERROR);
00084 }
00085 if ((oldface = MSD_getFace(s, f1)) == FNIL)
00086 {
00087 fprintf(stderr, MEN_NaoEncontrouFace, NAM_MEF, f1, sn);
00088 return(ERROR);
00089 }
00090 if ((oldloop = MSD_getLoop(oldface, v1, &he1)) == LNIL)
00091 {
00092 fprintf(stderr, MEN_NaoEncontrouVertice, NAM_MEF, v1, sn);
00093 return(ERROR);
00094 }
00095 if (MSD_getLoop(oldface, v2, &he2) != oldloop)
00096 {
00097 fprintf(stderr, MEN_NaoEncontrouVertice, NAM_MEF, v2, sn);
00098 return(ERROR);
00099 }
00100 MSD_lowMEF(he1, he2, f2);
00101 return(SUCCESS);
00102 }
00103
00104 #ifdef __Extended_Version
00105 void MSD_execMEF2(void)
00106 {
00107 int sn, v1, v2, v3, v4, f1, f2, f3, f4, op;
00108
00109 while (10 != sscanf(restbuf, "%d %d %d %d %d %d %d %d %d %d", &sn, &op, &v1,
00110 &v2, &v3, &v4, &f1, &f2, &f3, &f4))
00111 {
00112 printf("mef2: Sid V1 V2 V3 V4 F1 F2 F3 F4\n");
00113 if (!lineins("? "))
00114 {
00115 return;
00116 }
00117 }
00118 MSD_highMEF2(sn, v1, v2, v3, v4, f1, f2, f3, f4);
00119 }
00120
00121 #endif
00122
00123 int MSD_highMEF2(Id sn, Id v1, Id v2, Id v3, Id v4, Id f1, Id f2, Id f3, Id f4)
00124 {
00125 SPTYPE s;
00126 FPTYPE ff1;
00127 FPTYPE ff2;
00128 FPTYPE ff3;
00129 HPTYPE he1;
00130 HPTYPE he2;
00131
00132 if ((s = MSD_getSolid(sn)) == SNIL)
00133 {
00134 fprintf(stderr, MEN_NaoEncontrouSolido, NAM_MEF2, sn);
00135 return(ERROR);
00136 }
00137 if ((ff1 = MSD_getFace(s, f1)) == FNIL)
00138 {
00139 fprintf(stderr, MEN_NaoEncontrouFace, NAM_MEF2, f1, sn);
00140 return(ERROR);
00141 }
00142 if ((ff2 = MSD_getFace(s, f2)) == FNIL)
00143 {
00144 fprintf(stderr, MEN_NaoEncontrouFace, NAM_MEF2, f2, sn);
00145 return(ERROR);
00146 }
00147 if ((ff3 = MSD_getFace(s, f3)) == FNIL)
00148 {
00149 fprintf(stderr, MEN_NaoEncontrouFace, NAM_MEF2, f3, sn);
00150 return(ERROR);
00151 }
00152 if (MSD_getHalfedge(ff1, ff2, v1, v3, &he1) == LNIL)
00153 {
00154 fprintf(stderr, MEN_NaoEncontrouArestaF, NAM_MEF2, v1, v3, f1, f2);
00155 return(ERROR);
00156 }
00157 if (MSD_getHalfedge(ff1, ff3, v2, v4, &he2) == LNIL)
00158 {
00159 fprintf(stderr, MEN_NaoEncontrouArestaF, NAM_MEF2, v2, v4, f1, f3);
00160 return(ERROR);
00161 }
00162 MSD_lowMEF(he1, he2, f4);
00163 return(SUCCESS);
00164 }
00165
00166 FPTYPE MSD_lowMEF(HPTYPE he1, HPTYPE he2, Id fn)
00167 {
00168 SPTYPE s;
00169 DPTYPE d;
00170 FPTYPE newface;
00171 FPTYPE oldface;
00172 LPTYPE newloop;
00173 LPTYPE oldloop;
00174 EPTYPE newedge;
00175 HPTYPE nhe1;
00176 HPTYPE nhe2;
00177 HPTYPE he;
00178 Id v1, v2, f1;
00179 Id v3, v4, f2, f3;
00180
00181 #ifdef __curved_env
00182 short copy_tag, is_wire;
00183 vector v1;
00184 #endif
00185
00186 #ifdef __VIRTUAL_MEM
00187 MSD_checkHalfedge(he1);
00188 MSD_checkHalfedge(he2);
00189 #endif
00190
00191 oldloop = HalWLoop(he1);
00192 oldface = LooLFace(oldloop);
00193 d = FacFShell(oldface);
00194 s = SheSSolids(d);
00195
00196 assert1(HalWLoop(he1) == HalWLoop(he2));
00197 assert1(findl(oldloop, oldface) == 1);
00198 assert1(findf(d, oldface) == 1);
00199 #ifdef __curved_env
00200 copy_tag = he1 == he2 ? 0 : 1;
00201 is_wire = 0;
00202 if ((he1->prv == mate(he1)) && (he2->prv == mate(he2)) &&
00203 (s->sfaces->nextf == (FTYPE *)NIL))
00204 {
00205 is_wire = 1;
00206 }
00207 #endif
00208
00209 if (Gdoinversion)
00210 {
00211 v1 = VerVertexNo(HalVtx(he1));
00212 v2 = VerVertexNo(HalVtx(he2));
00213 f1 = FacFaceNo(oldface);
00214 MSD_lowAddEulerOp(s, KEF, v1, v2, f1, fn, 0, 0, 0, 0, 0.0, 0.0, 0.0, 0.0);
00215 }
00216
00217 if (Gdodirect)
00218 {
00219 v1 = VerVertexNo(HalVtx(he1));
00220 v2 = VerVertexNo(HalVtx(he2));
00221 v3 = VerVertexNo(HalVtx(HalNxt(he1)));
00222 v4 = VerVertexNo(HalVtx(HalNxt(he2)));
00223 f1 = FacFaceNo(oldface);
00224 f2 = HalEdg(he1) == ENIL ? f1 : FacFaceNo(LooLFace(HalWLoop(mate(he1))));
00225 f3 = HalEdg(he2) == ENIL ? f1 : FacFaceNo(LooLFace(HalWLoop(mate(he2))));
00226 MSD_lowAddEulerOpRedo(s, MEF, v1, v2, v3, v4, f1, f2, f3, fn, 0.0, 0.0, 0.0, 0.0);
00227 }
00228
00229 newface = (FPTYPE)MSD_lowNewElement(FACE, (NPTYPE)d);
00230 newloop = (LPTYPE)MSD_lowNewElement(LOOP, (NPTYPE)newface);
00231 newedge = (EPTYPE)MSD_lowNewElement(EDGE, (NPTYPE)d);
00232 EdgEdgeNo(newedge) = -1;
00233 FacFaceNo(newface) = fn;
00234 FacFLOut(newface) = newloop;
00235
00236 for (LooLength(newloop) = 0, he = he1; he != he2; he = HalNxt(he))
00237 {
00238 HalWLoop(he) = newloop;
00239 LooLength(newloop)++;
00240 LooLength(oldloop)--;
00241 }
00242
00243 nhe1 = MSD_lowAddHE(newedge, HalVtx(he2), he1, MINUS);
00244 nhe2 = MSD_lowAddHE(newedge, HalVtx(he1), he2, PLUS);
00245
00246 HalNxt(HalPrv(nhe1)) = nhe2;
00247 HalNxt(HalPrv(nhe2)) = nhe1;
00248 he = HalPrv(nhe1);
00249 HalPrv(nhe1) = HalPrv(nhe2);
00250 HalPrv(nhe2) = he;
00251
00252 LooLEdg(newloop) = nhe1;
00253 LooLEdg(oldloop) = nhe2;
00254
00255 #ifdef __curved_env
00256 if (!is_wire && (copy_tag == 1) && (oldface->surf != (SUSURF *)NIL))
00257 {
00258 lsettag(newface, oldface->surf);
00259 }
00260 if (is_wire)
00261 {
00262 faceeq(newface->flout, v1);
00263 lsettag(newface, newplane(++maxst, v1[0], v1[1], v1[2], v1[3]));
00264 faceeq(oldface->flout, v1);
00265 lsettag(oldface, newplane(++maxst, v1[0], v1[1], v1[2], v1[3]));
00266 }
00267 #endif
00268
00269 MSD_addEulCoef(s, 0, 1, 0, 1, 0, 0);
00270
00271 if (SetTrace)
00272 {
00273 if (TraceSolid(s))
00274 {
00275 fprintf(stderr, MEN_MSDAbortou, NAM_MEF);
00276 exit(1);
00277 }
00278 }
00279 return(newface);
00280 }
00281
00282 int MSD_execStructMEF(EulerOp *op)
00283 {
00284 Id sid, id1, id2, id3, id4, id5, id6, id7, id8;
00285
00286 sid = op->solidno;
00287 id1 = op->ip1;
00288 id2 = op->ip2;
00289 id3 = op->ip3;
00290 id4 = op->ip4;
00291 id5 = op->ip5;
00292 id6 = op->ip6;
00293 id7 = op->ip7;
00294 id8 = op->ip8;
00295 return(MSD_highMEF2(sid, id1, id2, id3, id4, id5, id6, id7, id8));
00296 }
00297
00298 void MSD_printStructMEF(FILE *f, EulerOp *op)
00299 {
00300 Id sid, id1, id2, id3, id4, id5, id6, id7, id8;
00301 int com;
00302
00303 sid = op->solidno;
00304 id1 = op->ip1;
00305 id2 = op->ip2;
00306 id3 = op->ip3;
00307 id4 = op->ip4;
00308 id5 = op->ip5;
00309 id6 = op->ip6;
00310 id7 = op->ip7;
00311 id8 = op->ip8;
00312 com = op->CommandCount;
00313 fprintf(f, "mef %d %d %d %d %d %d %d %d %d %d\n", sid,
00314 com, id1, id2, id3, id4, id5, id6, id7, id8);
00315 }
00316
00317 int MSD_readStructMEF(FILE *f, EulerOp *op)
00318 {
00319 int solidno, ip1, ip2, ip3, ip4, ip5, ip6, ip7, ip8, com;
00320
00321 if (fscanf(f, "%d %d %d %d %d %d %d %d %d %d\n", &solidno, &com, &ip1, &ip2, &ip3,
00322 &ip4, &ip5, &ip6, &ip7, &ip8) != 10)
00323 {
00324 return(ERROR);
00325 }
00326 op->solidno = solidno;
00327 op->CommandCount = com;
00328 op->ip1 = ip1;
00329 op->ip2 = ip2;
00330 op->ip3 = ip3;
00331 op->ip4 = ip4;
00332 op->ip5 = ip5;
00333 op->ip6 = ip6;
00334 op->ip7 = ip7;
00335 op->ip8 = ip8;
00336 op->fp1 = 0.0;
00337 op->fp2 = 0.0;
00338 op->fp3 = 0.0;
00339 op->fp4 = 0.0;
00340 return(SUCCESS);
00341 }
00342
00343 int MSD_modifyStructMEF(EulerOp *op)
00344 {
00345 int solidno;
00346
00347 solidno = MSD_lowCheckTranslatorTable(op->solidno, SOLID);
00348 op->solidno = solidno;
00349 return(SUCCESS);
00350 }