How do I add a wave to a curved surface composed of triangle primatives in C++?

Question

I want to preface this post: This is perhaps more of a math question than a coding question.

I am developing a plant (lettuce) model which involves somewhat complex geometry. At this stage I have a surface curved in 2 dimensions but now I want to add waviness to this curved surface but am having a hard time envisioning how to do so. The surface is made of triangle primatives, the primatives take xyz vectors to encode location of vertices. I am using an API termed HELIOS to develop this procedural model of lettuce. I essentially created the surface with for loops and the sine function. Disclaimer: I do not have a strong background in geometry, computer graphics, or C++.

Here is the relevant code:

#include "Context.h"
#include "Visualizer.h"

using namespace helios;
using namespace std;


vector addLeaf(float leaf_length, float leaf_width, float leaf_bend_x, float leaf_bend_y, float rotation_z, float rotation_x, float displacement, float radius, Context* context ) {

    std::vector UUIDs;

//    float leaf_length = 10;
    float Nz = 10; //  params.s1_leaf_subdivisions ; number of times to split the dimension
    float dz = leaf_length / Nz; // length of each subdivision

//    float leaf_width = 10;
    float Ny = 10; //  params.s1_leaf_subdivisions ; number of times to split the dimension
    float dy = leaf_width / Ny; // length of each subdivision

// leaf wave
//    float A_3 = leaf_length * float(0.5); // Half waves on the leaf 10
//    float A_2 = leaf_length * float(0.1); // amplitude 0.25

    float A_3 = 1; // Half waves on the leaf 10
    float A_2 = 1; // amplitude 0.25
    float leaf_amplitude = leaf_length / float(10);

// the 2 * dx extends the sine wave a bit beyond 1/2 wavelength so base of leaves come together
    for (int i = 0; i < Nz + (2 * dz); i++) {
        for (float j = 0; j < Ny; j++) {

            float z = i * dz; //for each subdivision in z define Z coord
            float y = j * dy; //for each subdivision in y define y coord
            float x = 0; // we will also need an x coord
            float sz = dz; // the next step in z will be equal to a subdivision in z
            float sy = dy; // the next step in y will be equal to a subdivision in y

            float z_i = z * M_PI / (Nz * dz); // the second coord for z is z_i needed to define a triangle primitive
            float sz_i = (z + sz) * M_PI / (Nz * dz); //

            // this would be y_1 in sorghum model
            float y_i = (y * M_PI / (Ny * dy)) / (A_3); // the second coord for y is y_i needed to define a triangle primitive
            float sy_i = ((y + sy) * M_PI / (Ny * dy)) / (A_3);

            //waviness of leaf
            float leaf_wave_1;
            float leaf_wave_2;
            float leaf_wave_3;
            float leaf_wave_4;

            if (j == 0) {
                leaf_wave_1 = A_2 * sin(z_i);
                leaf_wave_2 = A_2 * sin(sz_i);
            } else {
                leaf_wave_1 = A_2 * sin(z_i);
                leaf_wave_2 = A_2 * sin(sz_i);
            }

            // Now define x based on z,y and add leaf bend in x and y
            x = leaf_bend_x * sin(z_i);
            x = ((x*radius + displacement + (leaf_bend_y * sin(y_i))) / 2) + leaf_wave_1;

            vec3 v0(x*radius + displacement, y, z);

            x = leaf_bend_x * sin(sz_i);
            x = ((x*radius + displacement + (leaf_bend_y * sin(y_i))) / 2) + leaf_wave_2;

            vec3 v1(x*radius + displacement, y, z + sz);

            if (j == Nz - 1) {
                leaf_wave_3 = sin(sz_i) * A_2;
                leaf_wave_4 = sin(z_i) * A_2;
            } else {
                leaf_wave_3 = sin(sz_i) * A_2;
                leaf_wave_4 = sin(z_i) * A_2;
            }

            x = leaf_bend_x * sin(sz_i);
            x = ((x*radius + displacement + (leaf_bend_y * sin(sy_i))) / 2) + leaf_wave_3 ;


            vec3 v2(x*radius + displacement, y + sy, z + sz);

            x = leaf_bend_x * sin(z_i);
            x = ((x*radius + displacement + (leaf_bend_y * sin(sy_i))) / 2) + leaf_wave_4 ;

            vec3 v3(x*radius + displacement, y + sy, z);

            // set of two triangles which form a rectangle or square as subunits of leaf
            UUIDs.push_back(context->addTriangle(v0, v1, v2, RGB::cyan));
            UUIDs.push_back(context->addTriangle(v0, v2, v3, RGB::magenta));
        }
    }

     return UUIDs;
}




// call to functions and build lettuce geometries
int main( void ){

    Context context;
    float leaf_length = 10;
    float leaf_width = 10;
    float radius = 1; // additional control leaf curvature

// add leaves one by one; 'i' here is # of leaves external to whorl
    for (int i = 0; i  < 6; i++) {

        if (i == 0)addLeaf(leaf_length, leaf_width,  0.5*leaf_length, 0.5*leaf_width, 4 * M_PI / 9*i, 0, i/5, radius, &context);
//        if (i == 1)addLeaf(leaf_length, leaf_width,  0.5*leaf_length, 0.5*leaf_width, 4 * M_PI / 9*i, -M_PI/ 20, i/5, radius, &context);
//        if (i == 2)addLeaf(leaf_length, leaf_width,  0.5*leaf_length, 0.5*leaf_width, 4 * M_PI / 9*i, -M_PI/ 10, i/5, radius, &context);
//        if (i == 3)addLeaf(leaf_length, leaf_width,  0.5*leaf_length, 0.5*leaf_width, 4 * M_PI / 9*i, -M_PI/ 9, i/5, radius, &context);
//        if (i == 4)addLeaf(leaf_length, leaf_width,  0.5*leaf_length, 0.5*leaf_width, 4 * M_PI / 9*i, -M_PI/ 7, i/5, radius, &context);
//        if (i == 5)addLeaf(leaf_length, leaf_width,  0.5*leaf_length, 0.5*leaf_width, 4 * M_PI / 9*i, -M_PI/ 5, i/5, radius, &context);

    }
    Visualizer visualizer(800);
    visualizer.buildContextGeometry(&context);
    visualizer.setLightingModel(Visualizer::LIGHTING_PHONG);
    visualizer.plotInteractive();

}

I tried to use a sine function and an additional for loop to create a series of values to add to the X coordinate of the triangles but did not obtain the result I was looking for.

How do I add a wave to a curved surface composed of triangle primatives in C++?

Answers (1)

Related Questions