双三次bezier曲面片c代码

以下是一个简单的C代码片段，用于计算双三次Bezier曲面片上的点坐标：

#include <stdio.h>
#include <math.h>

#define Rows 4
#define Cols 4

typedef struct {
    double x, y, z;
} Point;

typedef struct {
    Point p[Rows][Cols];
} BezierSurface;

double Bernstein(int i, int n, double u)
{
    double coeff = 1.0;
    for (int j = 1; j <= i; j++) {
        coeff *= (u - 1.0) / j;
    }
    for (int j = 1; j <= n - i; j++) {
        coeff *= u / j;
    }
    return coeff;
}

Point EvaluateBezierSurface(BezierSurface surface, double u, double v)
{
    Point p;
    p.x = 0.0;
    p.y = 0.0;
    p.z = 0.0;
    for (int i = 0; i < Rows; i++) {
        for (int j = 0; j < Cols; j++) {
            double coeff = Bernstein(i, Rows - 1, u) * Bernstein(j, Cols - 1, v);
            p.x += coeff * surface.p[i][j].x;
            p.y += coeff * surface.p[i][j].y;
            p.z += coeff * surface.p[i][j].z;
        }
    }
    return p;
}

int main()
{
    BezierSurface surface;
    // 初始化控制点
    for (int i = 0; i < Rows; i++) {
        for (int j = 0; j < Cols; j++) {
            surface.p[i][j].x = i;
            surface.p[i][j].y = j;
            surface.p[i][j].z = sin(i) * cos(j);
        }
    }
    // 在u = 0.5，v = 0.5的位置计算曲面上的点
    Point p = EvaluateBezierSurface(surface, 0.5, 0.5);
    printf("x = %f, y = %f, z = %f\n", p.x, p.y, p.z);
    return 0;
}

在这个例子中，我们定义了一个`BezierSurface`类型来表示双三次Bezier曲面片，并定义了一个`Point`类型来表示三维点。我们使用`Bernstein`函数来计算Bezier曲线上的权重值，然后使用`EvaluateBezierSurface`函数来计算曲面上的点坐标。在`main`函数中，我们初始化了控制点，并在u = 0.5，v = 0.5的位置计算了曲面上的点。