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 #include <math.h>
00043 #include "mensagem.h"
00044 #include "memvirtu.h"
00045 #include "lowparam.h"
00046 #include "lowmacro.h"
00047 #include "lowsolid.h"
00048 #include "eulerops.h"
00049 #include "mancommd.h"
00050 #include "vectorop.h"
00051 #include "curvgeom.h"
00052
00053 #ifdef __curved_env__
00054
00055 #ifdef __Extended_Version
00056 void MSD_execEvaluateCurveLine(void)
00057 {
00058 int sn, fn, v1, v2, nsegm;
00059 char sname[20];
00060
00061 while (5 != sscanf(restbuf, "%s %d %d %d %d", sname, &fn, &v1, &v2, &nsegm))
00062 {
00063 printf("EVAL_LINE: sname fn v1 v2 nsegm\n");
00064 if (!lineins("? "))
00065 {
00066 return;
00067 }
00068 }
00069 if ((sn = MSD_getSolidIdFromName(sname)) == ERROR)
00070 {
00071 fprintf(stderr, MEN_NaoEncontrouSolidNome, NAM_EvalEdge, sname);
00072 return;
00073 }
00074 MSD_highEvaluateCurveLine(sn, fn, v1, v2, nsegm);
00075 }
00076
00077 int MSD_highEvaluateCurveLine(Id sn, Id fn, Id v1, Id v2, int nsegm)
00078 {
00079 SPTYPE s;
00080 FPTYPE f;
00081 HPTYPE he;
00082
00083 if ((s = MSD_getSolid(sn)) == SNIL)
00084 {
00085 fprintf(stderr, MEN_NaoEncontrouSolido, NAM_EvalEdge, sn);
00086 return(ERROR);
00087 }
00088 if ((f = MSD_getFace(s, fn)) == FNIL)
00089 {
00090 fprintf(stderr, MEN_NaoEncontrouFace, NAM_EvalEdge, fn, sn);
00091 return(ERROR);
00092 }
00093 if (MSD_getEdge(f, v1, v2, &he) == LNIL)
00094 {
00095 fprintf(stderr, MEN_NaoEncontrouArestaV, NAM_EvalEdge, v1, v2, sn);
00096 return(ERROR);
00097 }
00098 if (nsegm < 2)
00099 {
00100 fprintf(stderr, MEN_NumeroDeSegmentosInvalido, NAM_EvalEdge);
00101 return(ERROR);
00102 }
00103 MSD_lowEvaluateCurveLine(he, nsegm);
00104 return(SUCCESS);
00105 }
00106
00107 #endif
00108
00109 void MSD_lowEvaluateCurveLine(HPTYPE he, int nsegm)
00110 {
00111 CURVE *cu;
00112 LPOINT *p;
00113 int j;
00114 real x, y, z, deltax, deltay, deltaz;
00115
00116 maxv = MSD_getMaxVertexId(SheSSolids(FacFShell(LooLFace(HalWLoop(he)))));
00117 cu = HalCurv(he);
00118 he = HalNxt(he);
00119
00120 p = cu->lpnt;
00121
00122 Mc[0][0] = -1.0;
00123 Mc[0][1] = 3.0;
00124 Mc[0][2] = -3.0;
00125 Mc[0][3] = 1.0;
00126 Mc[1][0] = 3.0;
00127 Mc[1][1] = -6.0;
00128 Mc[1][2] = 3.0;
00129 Mc[1][3] = 0.0;
00130 Mc[2][0] = -3.0;
00131 Mc[2][1] = 3.0;
00132 Mc[2][2] = 0.0;
00133 Mc[2][3] = 0.0;
00134 Mc[3][0] = 0.0;
00135 Mc[3][1] = 0.0;
00136 Mc[3][2] = 0.0;
00137 Mc[3][3] = 0.0;
00138
00139 mattranspose(Mc, Mc);
00140
00141 x = p->pnt[0];
00142 y = p->pnt[1];
00143 z = p->pnt[2];
00144
00145 deltax = p->pnxt->pnt[0] - p->pnt[0] / nsegm;
00146 deltay = p->pnxt->pnt[1] - p->pnt[1] / nsegm;
00147 deltaz = p->pnxt->pnt[2] - p->pnt[2] / nsegm;
00148
00149 for (j = 1; j <= nsegm - 1; j++)
00150 {
00151 x += deltax;
00152 y += deltay;
00153 z += deltaz;
00154
00155 MSD_lowMEV(mate(HalPrv(he)), he, ++maxv, x, y, z);
00156 MSD_lowSetCurveAtributeToEdge(HalPrv(he), cu);
00157 }
00158 }
00159
00160 #endif