变换一个三角形并用C

时间:2019-05-22 19:13:05

标签: c geometry

给出一个三角形,我必须按照以下方式将一个三角形转换为另一个三角形:相对于原点旋转给定角度的顶点,然后相对于原点应用均值。并将每个三角形保存在n + 1个矢量中,然后根据顶点a的自变量(0到2π之间的角度)对矢量进行排序,并直接插入。

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
typedef struct{
float x;
float y; }point;
typedef struct {
point a;
point b;
point c; }triangle;

      float angle(triangle *t,int i){

float a,b,c,alpha;
 a=(t[i].b.x-t[i].c.x)*(t[i].b.x-t[i].c.x)+(t[i].b.y-t[i].c.y)*(t[i].b.y-t[i].c.y);
 b=(t[i].a.x-t[i].b.x)*(t[i].a.x-t[i].b.x)+(t[i].a.y-t[i].b.y)*(t[i].a.y-t[i].b.y);
   c=(t[i].a.x-t[i].c.x)*(t[i].a.x-t[i].c.x)+(t[i].a.y-t[i].c.y)*(t[i].a.y-t[i].c.y);
  alpha=acos((b*b+c*c-a*a)/(2*b*c));
 return alpha;} 

   void move(int r,float angle,triangle *t,int i){

float alpha;
 // rotation
     alpha=angle*M_PI/180;
        t[i+1].a.x=t[i].a.x*cos(alpha)-t[i].a.y*sin(alpha);
        t[i+1].a.y=t[i].a.x*sin(alpha)+t[i].a.y*cos(alpha);
    t[i+1].b.x=t[i].b.x*cos(alpha)-t[i].b.y*sin(alpha);
    t[i+1].b.y=t[i].b.x*sin(alpha)+t[i].b.y*cos(alpha);
    t[i+1].c.x=t[i].c.x*cos(alpha)-t[i].c.y*sin(alpha);
    t[i+1].c.y=t[i].c.x*sin(alpha)+t[i].c.y*cos(alpha);

// homotecy
    t[i+1].a.x=t[i+1].a.x*r;
    t[i+1].a.y=t[i+1].a.y*r;
    t[i+1].b.x=t[i+1].b.x*r;
    t[i+1].b.y=t[i+1].b.y*r;
    t[i+1].c.x=t[i+1].c.x*r;
     t[i+1].c.y=t[i+1].c.y*r;
       return;}

void sort(triangle *t,int n){

float a,b;
int k,i;
for(i=1;i<=n-1;i++){
    a=angle(t,i);
    k=i-1;
    b=angle(t,k);
    while(a>=0 && a<=(M_PI*2) && b>a){
        t[k+1]=t[k];
        k--;
    }
    t[k+1]=t[i];
}
return;}

 int main(void){

triangle *t;
float angle;
int r,n,i;
char so[50];
FILE *s;

printf("n?\n");
scanf("%d",&n);
printf("r?\n");
scanf("%d",&r);
printf("angle?\n");
scanf("%f",&angle);
t=(triangle*)malloc((n+1)*sizeof(point));
if(t==NULL) exit(1);
printf("Triangle?\n");
scanf("%f %f %f %f %f %f",&t[0].a.x,&t[0].a.y,&t[0].b.x,&t[0].b.y,&t[0].c.x,&t[0].c.y);
for(i=0;i<n;i++){
    move(r,angle,t,i);
    sort(t,n);
}
printf("file?\n");
scanf("%s",so);
s=fopen(so,"w");
if(s==NULL) exit(1);
for(i=0;i<n;i++){
    fprintf(s,"%9.5f %9.5f %19.5f %9.5f %9.5f %9.5f\n",t[i].a.x,t[i].a.y,t[i].b.x,t[i].b.y,t[i].c.x,t[i].c.y);
}
fclose(s);  
free(t);
   return 0;}

执行时显示:

free(): invalid next size (normal)

1 个答案:

答案 0 :(得分:2)

更改以下行:

t=(triangle*)malloc((n+1)*sizeof(point));

t=(triangle*)malloc((n+1)*sizeof(triangle));甚至

t = malloc((n+1) * sizeof *t);,以避免强制转换和处理t类型。

在分配三角形的同时分配点。然后,您将您不拥有的内容写入内存并破坏了内存。这可能会导致分段错误。在您的情况下,您很幸运(或者,如果您认为您之前已经捕获了该错误,那么可能会很不幸),并且在释放分配的内存时遇到了错误。根据我的经验,释放内存时的错误几乎总是导致内存损坏,如果不小心,有时会释放两次错误。