iCAx开思工具箱

标题: 请教一个生成曲面的程序 [打印本页]

作者: maomao521 时间: 2006-5-9 10:53
标题: 请教一个生成曲面的程序
请教一个程序，该程序是.通过给定点生成B样条曲面。原来的例子是给定十六个点，我改成了36个点，改了相应得参数，编译连接都没有错误，但是到UG里面执行时出现下面错误，显示内部错误：内存访问冲突，恳请高手指点，
#include <stdio.h>
#include <stdlib.h>
#include <uf.h>
#include <uf_part.h>
#include <uf_defs.h>
#include <uf_obj.h>
#include <uf_modl.h>

#define UF_CALL(X) (report( __FILE__, __LINE__, #X, (X)))
static int report( char *file, int line, char *call, int irc)
{
  if (irc)
  {
   char messg[133];
   printf("%s, line %d:  %s\n", file, line, call);
   (UF_get_fail_message(irc, messg)) ?
   printf(" returned a %d\n", irc) :
   printf(" returned error %d:  %s\n", irc, messg);
  }
  return(irc);
}
static void do_ugopen_api(void)
{
  char *prtnam = "create_bsurf";
  /* Points for B-surface  */
  double points[36][3] = {
            0.000000,0.000000,0.005914,
            0.000000,0.200000,0.005074,
            0.000000,0.400000,0.002990,
            0.000000,0.600000,0.001533,
            0.000000,0.800000,0.000828,
            0.000000,1.000000,0.000577,
            0.200000,0.000000,0.003775,
            0.200000,0.200000,0.004451,
            0.200000,0.400000,0.004584,
            0.200000,0.600000,0.005139,
            0.200000,0.800000,0.006252,
            0.200000,1.000000,0.007520,
            0.400000,0.000000,-0.000929,
            0.400000,0.200000,0.001315,
            0.400000,0.400000,0.006451,
            0.400000,0.600000,0.017956,
            0.400000,0.800000,0.038411,
            0.400000,1.000000,0.060684,
            0.600000,0.000000,-0.009034,
            0.600000,0.200000,-0.006401,
            0.600000,0.400000,0.006438,
            0.600000,0.600000,0.041146,
            0.600000,0.800000,0.112490,
            0.600000,1.000000,0.195991,
            0.800000,0.000000,0.008290,
            0.800000,0.200000,0.005289,
            0.800000,0.400000,0.014894,
            0.800000,0.600000,0.031797,
            0.800000,0.800000,0.052069,
            0.800000,1.000000,0.071069,
            1.000000,0.000000,0.159674,
            1.000000,0.200000,0.076535,
            1.000000,0.400000,0.032681,
            1.000000,0.600000,0.013937,
            1.000000,0.800000,0.006826,
            1.000000,1.000000,0.004618,

/*             3.5, 2.5, 0.0,
   4.0, 2.6, 0.2,
   4.5, 2.6, 0.2,
   5.0, 2.5, 0.0,
   3.3, 3.0, 0.2,
   4.1, 3.2, 0.3,
   4.5, 3.2, 0.3,
   5.0, 3.0, 0.2,
   3.0, 3.5, 0.2,
   3.8, 3.6, 0.3,
   4.5, 3.6, 0.3,
   5.0, 3.5, 0.2,
   2.6, 4.0, 0.0,
   3.3, 4.0, 0.2,
   4.5, 4.0, 0.2,
   5.0, 4.0, 0.0 */
   };

  int create_mode = 1;
  int u_closed_status = 0;
  int v_closed_status = 0;
  int u_degree = 3;
  int v_degree = 3;
  int num_rows = 6;
  int pts_per_row[6] = {6, 6, 6, 6};
  int ii, jj, kk, indx;
  tag_t part_tag;
  tag_t bsurf_obj_id = NULL_TAG;
  UF_MODL_bsurf_row_info_t *pts_info_per_row; /* pts info for each row */
  UF_MODL_bsurface_t bsurf;
  /* Open a new part */
  UF_CALL( UF_PART_new(prtnam, METRIC, &part_tag) );
  pts_info_per_row = ( UF_MODL_bsurf_row_info_t *)malloc( num_rows *
                     sizeof(UF_MODL_bsurf_row_info_t) );
  /* allocate and load point and other information for each row of points */
  indx = 0;
  for (ii=0; ii<num_rows; ii++)
  {
pts_info_per_row[ii].num_points = pts_per_row[ii];  /* can vary */
pts_info_per_row[ii].points = (double *)malloc( pts_per_row[ii] * 3 *
                                 sizeof(double) );
pts_info_per_row[ii].weight = (double *)malloc( pts_per_row[ii] *
                                 sizeof(double) );
/* load up point and weight info */
for (jj=0, kk = 0; kk<pts_per_row[ii]; kk++, jj+=3, indx++)
{
   pts_info_per_row[ii].points[jj] = points[indx][0];
   pts_info_per_row[ii].points[jj+1] = points[indx][1];
   pts_info_per_row[ii].points[jj+2] = points[indx][2];
   pts_info_per_row[ii].weight[kk] = 1.0;
}
  }
  /* create B-surface */
  UF_CALL(UF_MODL_create_bsurf_thru_pts( create_mode, u_closed_status,
                                       v_closed_status, u_degree,
                                       v_degree, num_rows, pts_info_per_row,
                                 &bsurf_obj_id));

  if ( bsurf_obj_id == NULL_TAG )
printf("*** ERROR CREATING B-SURFACE ***\n");

  /* free allocated memory */
  for (ii=0; ii<num_rows; ii++)
  {
free(pts_info_per_row[ii].points);
free(pts_info_per_row[ii].weight);
  }
  free(pts_info_per_row);
  /* ask the B-surface information */
  UF_CALL(UF_MODL_ask_bsurf( bsurf_obj_id, &bsurf ));

  /* validate the surface */
  if ( (bsurf.order_u != u_degree+1) || (bsurf.order_v != v_degree+1) )
printf("*** ERROR ASKING B-SURFACE ***\n");

  UF_CALL(UF_MODL_free_bsurf_data(&bsurf));
}
/*ARGSUSED*/
void ufusr(char *param, int *retcode, int param_len)
{
  if (!UF_CALL(UF_initialize()))
  {
do_ugopen_api();
UF_CALL(UF_terminate());
  }
}
int ufusr_ask_unload(void)
{
  return (UF_UNLOAD_IMMEDIATELY);
}

欢迎光临 iCAx开思工具箱 (https://t.icax.org/)