数学デッサン教室 – ページ 2 – プログラミング de 落描き

2021年4月14日2021年4月14日

Change ring size. Draw by openFrameworks

[ Video ]

[ About ]
Change ring size.

[ Source ]

#pragma once
#include "ofMain.h"

class ofApp : public ofBaseApp {

public:
	void setup();
	void update();
	void draw();

	void keyPressed(int key) {};
	void keyReleased(int key) {};
	void mouseMoved(int x, int y) {};
	void mouseDragged(int x, int y, int button) {};
	void mousePressed(int x, int y, int button) {};
	void mouseReleased(int x, int y, int button) {};
	void mouseEntered(int x, int y) {};
	void mouseExited(int x, int y) {};
	void windowResized(int w, int h) {};
	void dragEvent(ofDragInfo dragInfo) {};
	void gotMessage(ofMessage msg) {};

	glm::vec3 make_point(float R, float r, float u, float v);

	ofEasyCam cam;
	ofMesh face, frame;

	float r;
};

#pragma once

#include "ofMain.h"

class ofApp : public ofBaseApp {

public:

void setup();

void update();

void draw();

void keyPressed(int key) {};

void keyReleased(int key) {};

void mouseMoved(int x, int y) {};

void mouseDragged(int x, int y, int button) {};

void mousePressed(int x, int y, int button) {};

void mouseReleased(int x, int y, int button) {};

void mouseEntered(int x, int y) {};

void mouseExited(int x, int y) {};

void windowResized(int w, int h) {};

void dragEvent(ofDragInfo dragInfo) {};

void gotMessage(ofMessage msg) {};

glm::vec3 make_point(float R, float r, float u, float v);

ofEasyCam cam;

ofMesh face, frame;

float r;

};

#include "ofApp.h"	

//--------------------------------------------------------------
void ofApp::setup() {

	ofSetFrameRate(60);
	ofSetWindowTitle("openFrameworks");

	ofBackground(239);
	ofSetLineWidth(1.5);
	ofEnableDepthTest();

	this->frame.setMode(ofPrimitiveMode::OF_PRIMITIVE_LINES);
}

//--------------------------------------------------------------
void ofApp::update() {

	ofSeedRandom(39);

	this->face.clear();
	this->frame.clear();

	float R = 300;
	this->r = ofMap(ofNoise(ofRandom(1000), ofGetFrameNum() * 0.008), 0, 1, 10, 130);
	int v_span = 8;
	int u_span = 10;

	for (int v = 0; v < 360; v += v_span) {

		for (auto u = 0; u < 360; u += u_span) {

			auto noise_value = ofNoise(glm::vec4(this->make_point(R, 50, u + u_span * 0.5, v) * 0.005, ofGetFrameNum() * 0.01));

			if (noise_value > 0.6) {

				vector<glm::vec3> vertices;
				vertices.push_back(this->make_point(R, this->r, u, v - v_span * 0.5));
				vertices.push_back(this->make_point(R, this->r, u + u_span, v - v_span * 0.5));
				vertices.push_back(this->make_point(R, this->r, u + u_span, v + v_span * 0.5));
				vertices.push_back(this->make_point(R, this->r, u, v + v_span * 0.5));

				int index = face.getNumVertices();
				this->face.addVertices(vertices);

				this->face.addIndex(index + 0); face.addIndex(index + 1); face.addIndex(index + 2);
				this->face.addIndex(index + 0); face.addIndex(index + 2); face.addIndex(index + 3);

				this->frame.addVertices(vertices);

				this->frame.addIndex(index + 0); this->frame.addIndex(index + 1);
				this->frame.addIndex(index + 1); this->frame.addIndex(index + 2);
				this->frame.addIndex(index + 2); this->frame.addIndex(index + 3);
				this->frame.addIndex(index + 3); this->frame.addIndex(index + 0);
			}
		}
	}
}

//--------------------------------------------------------------
void ofApp::draw() {

	this->cam.begin();

	int number = 300;
	vector<glm::vec3> location_list(number);
	for (int i = 0; i < number; i++) {

		auto R = 300;
		auto tmp_r = ofRandom(10, this->r * 0.8);
		auto u = ofRandom(360) + ofGetFrameNum() * ofRandom(0.5, 1.0) * (ofRandom(1) < 0.5 ? 1 : -1);
		auto v = ofRandom(360) + ofGetFrameNum() * ofRandom(0.5, 1.0) * (ofRandom(1) < 0.5 ? 1 : -1);

		location_list[i] = this->make_point(R, tmp_r, u, v);
	}

	ofFill();
	ofSetColor(39);
	for (auto& location : location_list) {

		ofDrawSphere(location, 3);
		for (auto& other : location_list) {

			if (location == other) { continue; }

			auto distance = glm::distance(location, other);
			if (distance < 40) {

				ofDrawLine(location, other);
			}
		}
	}

	ofSetColor(239);
	this->face.draw();

	ofSetColor(39);
	this->frame.drawWireframe();

	this->cam.end();
}

//--------------------------------------------------------------
glm::vec3 ofApp::make_point(float R, float r, float u, float v) {

	// 数学デッサン教室 描いて楽しむ数学たち　P.31

	u *= DEG_TO_RAD;
	v *= DEG_TO_RAD;

	auto x = (R + r * cos(u)) * cos(v);
	auto y = (R + r * cos(u)) * sin(v);
	auto z = r * sin(u);

	return glm::vec3(x, y, z);
}

//--------------------------------------------------------------
int main() {

	ofSetupOpenGL(720, 720, OF_WINDOW);
	ofRunApp(new ofApp());
}

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

#include "ofApp.h"

//--------------------------------------------------------------

void ofApp::setup() {

ofSetFrameRate(60);

ofSetWindowTitle("openFrameworks");

ofBackground(239);

ofSetLineWidth(1.5);

ofEnableDepthTest();

this->frame.setMode(ofPrimitiveMode::OF_PRIMITIVE_LINES);

}

//--------------------------------------------------------------

void ofApp::update() {

ofSeedRandom(39);

this->face.clear();

this->frame.clear();

float R = 300;

this->r = ofMap(ofNoise(ofRandom(1000), ofGetFrameNum() * 0.008), 0, 1, 10, 130);

int v_span = 8;

int u_span = 10;

for (int v = 0; v < 360; v += v_span) {

for (auto u = 0; u < 360; u += u_span) {

auto noise_value = ofNoise(glm::vec4(this->make_point(R, 50, u + u_span * 0.5, v) * 0.005, ofGetFrameNum() * 0.01));

if (noise_value > 0.6) {

vector<glm::vec3> vertices;

vertices.push_back(this->make_point(R, this->r, u, v - v_span * 0.5));

vertices.push_back(this->make_point(R, this->r, u + u_span, v - v_span * 0.5));

vertices.push_back(this->make_point(R, this->r, u + u_span, v + v_span * 0.5));

vertices.push_back(this->make_point(R, this->r, u, v + v_span * 0.5));

int index = face.getNumVertices();

this->face.addVertices(vertices);

this->face.addIndex(index + 0); face.addIndex(index + 1); face.addIndex(index + 2);

this->face.addIndex(index + 0); face.addIndex(index + 2); face.addIndex(index + 3);

this->frame.addVertices(vertices);

this->frame.addIndex(index + 0); this->frame.addIndex(index + 1);

this->frame.addIndex(index + 1); this->frame.addIndex(index + 2);

this->frame.addIndex(index + 2); this->frame.addIndex(index + 3);

this->frame.addIndex(index + 3); this->frame.addIndex(index + 0);

}

//--------------------------------------------------------------

void ofApp::draw() {

this->cam.begin();

int number = 300;

vector<glm::vec3> location_list(number);

for (int i = 0; i < number; i++) {

auto R = 300;

auto tmp_r = ofRandom(10, this->r * 0.8);

auto u = ofRandom(360) + ofGetFrameNum() * ofRandom(0.5, 1.0) * (ofRandom(1) < 0.5 ? 1 : -1);

auto v = ofRandom(360) + ofGetFrameNum() * ofRandom(0.5, 1.0) * (ofRandom(1) < 0.5 ? 1 : -1);

location_list[i] = this->make_point(R, tmp_r, u, v);

}

ofFill();

ofSetColor(39);

for (auto& location : location_list) {

ofDrawSphere(location, 3);

for (auto& other : location_list) {

if (location == other) { continue; }

auto distance = glm::distance(location, other);

if (distance < 40) {

ofDrawLine(location, other);

}

ofSetColor(239);

this->face.draw();

ofSetColor(39);

this->frame.drawWireframe();

this->cam.end();

}

//--------------------------------------------------------------

glm::vec3 ofApp::make_point(float R, float r, float u, float v) {

// 数学デッサン教室描いて楽しむ数学たち　P.31

u *= DEG_TO_RAD;

v *= DEG_TO_RAD;

auto x = (R + r * cos(u)) * cos(v);

auto y = (R + r * cos(u)) * sin(v);

auto z = r * sin(u);

return glm::vec3(x, y, z);

}

//--------------------------------------------------------------

int main() {

ofSetupOpenGL(720, 720, OF_WINDOW);

ofRunApp(new ofApp());

}

[ Link ]
https://github.com/junkiyoshi/Insta20210414

2021年3月8日2021年3月7日

Vanishing torus. Draw by openFrameworks

[ Video ]

[ About ]
Vanishing torus.

[ Source ]

#pragma once
#include "ofMain.h"

class ofApp : public ofBaseApp {

public:
	void setup();
	void update();
	void draw();

	void keyPressed(int key) {};
	void keyReleased(int key) {};
	void mouseMoved(int x, int y) {};
	void mouseDragged(int x, int y, int button) {};
	void mousePressed(int x, int y, int button) {};
	void mouseReleased(int x, int y, int button) {};
	void windowResized(int w, int h) {};
	void dragEvent(ofDragInfo dragInfo) {};
	void gotMessage(ofMessage msg) {};

	glm::vec3 make_point(float R, float r, float u, float v);

	ofEasyCam cam;
	ofMesh face, frame;
};

#pragma once

#include "ofMain.h"

class ofApp : public ofBaseApp {

public:

void setup();

void update();

void draw();

void keyPressed(int key) {};

void keyReleased(int key) {};

void mouseMoved(int x, int y) {};

void mouseDragged(int x, int y, int button) {};

void mousePressed(int x, int y, int button) {};

void mouseReleased(int x, int y, int button) {};

void windowResized(int w, int h) {};

void dragEvent(ofDragInfo dragInfo) {};

void gotMessage(ofMessage msg) {};

glm::vec3 make_point(float R, float r, float u, float v);

ofEasyCam cam;

ofMesh face, frame;

};

#include "ofApp.h"

//--------------------------------------------------------------
void ofApp::setup() {

	ofSetFrameRate(30);
	ofSetWindowTitle("openFrameworks");

	ofBackground(239);
	ofSetLineWidth(3);
	ofEnableDepthTest();

	this->frame.setMode(ofPrimitiveMode::OF_PRIMITIVE_LINES);
}

//--------------------------------------------------------------
void ofApp::update() {

	ofSeedRandom(39);

	this->face.clear();
	this->frame.clear();

	float phi_deg_step = 6;
	float theta_deg_step = 6;

	int R_start = (ofGetFrameNum() * 4) % 350 + 30;
	for (int R_param = R_start; R_param < R_start + 350; R_param += 80) {

		int R = R_param % 350 + 30;
		float r = R * 0.13;

		auto rotation_y = glm::rotate(glm::mat4(), ofMap(ofNoise(ofRandom(1000), ofGetFrameNum() * 0.001), 0, 1, PI * -5, PI * 5), glm::vec3(0, 1, 0));
		auto rotation_x = glm::rotate(glm::mat4(), ofMap(ofNoise(ofRandom(1000), ofGetFrameNum() * 0.001), 0, 1, PI * -5, PI * 5), glm::vec3(1, 0, 0));

		for (float phi_deg = 0; phi_deg < 360; phi_deg += phi_deg_step) {

			for (float theta_deg = 0; theta_deg < 360; theta_deg += theta_deg_step) {

				float threshold = ofMap(R, 30, 350, 0, 1);
				auto noise_value = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg, phi_deg) * 0.005));
				if (noise_value < threshold) { continue; }

				auto noise_1 = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg - theta_deg_step, phi_deg) * 0.005));
				auto noise_2 = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg, phi_deg + phi_deg_step) * 0.005));
				auto noise_3 = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg, phi_deg - phi_deg_step) * 0.005));
				auto noise_4 = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg + theta_deg_step, phi_deg) * 0.005));

				auto index = this->face.getNumVertices();
				vector<glm::vec3> vertices;

				vertices.push_back(glm::vec3(this->make_point(R, r, theta_deg - theta_deg_step * 0.5, phi_deg - phi_deg_step * 0.5)));
				vertices.push_back(glm::vec3(this->make_point(R, r, theta_deg + theta_deg_step * 0.5, phi_deg - phi_deg_step * 0.5)));
				vertices.push_back(glm::vec3(this->make_point(R, r, theta_deg - theta_deg_step * 0.5, phi_deg + phi_deg_step * 0.5)));
				vertices.push_back(glm::vec3(this->make_point(R, r, theta_deg + theta_deg_step * 0.5, phi_deg + phi_deg_step * 0.5)));

				for (auto& vertex : vertices) {

					vertex = glm::vec4(vertex, 0) * rotation_y * rotation_x;
				}

				this->face.addVertices(vertices);

				this->face.addIndex(index + 0); this->face.addIndex(index + 1); this->face.addIndex(index + 3);
				this->face.addIndex(index + 0); this->face.addIndex(index + 3); this->face.addIndex(index + 2);

				if (noise_1 < threshold) {

					this->frame.addVertex(vertices[0]);
					this->frame.addVertex(vertices[2]);

					this->frame.addIndex(this->frame.getNumVertices() - 1);
					this->frame.addIndex(this->frame.getNumVertices() - 2);
				}

				if (noise_2 < threshold) {

					this->frame.addVertex(vertices[2]);
					this->frame.addVertex(vertices[3]);

					this->frame.addIndex(this->frame.getNumVertices() - 1);
					this->frame.addIndex(this->frame.getNumVertices() - 2);
				}

				if (noise_3 < threshold) {

					this->frame.addVertex(vertices[0]);
					this->frame.addVertex(vertices[1]);

					this->frame.addIndex(this->frame.getNumVertices() - 1);
					this->frame.addIndex(this->frame.getNumVertices() - 2);
				}

				if (noise_4 < threshold) {

					this->frame.addVertex(vertices[1]);
					this->frame.addVertex(vertices[3]);

					this->frame.addIndex(this->frame.getNumVertices() - 1);
					this->frame.addIndex(this->frame.getNumVertices() - 2);

				}
			}
		}
	}
}

//--------------------------------------------------------------
void ofApp::draw() {

	this->cam.begin();

	ofSetColor(255);
	this->frame.draw();

	ofSetColor(39);
	this->face.draw();

	this->cam.end();
}

//--------------------------------------------------------------
glm::vec3 ofApp::make_point(float R, float r, float u, float v) {

	// 数学デッサン教室 描いて楽しむ数学たち　P.31

	u *= DEG_TO_RAD;
	v *= DEG_TO_RAD;

	auto x = (R + r * cos(u)) * cos(v);
	auto y = (R + r * cos(u)) * sin(v);
	auto z = r * sin(u);

	return glm::vec3(x, y, z);
}


//--------------------------------------------------------------
int main() {

	ofSetupOpenGL(720, 720, OF_WINDOW);
	ofRunApp(new ofApp());
}

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

142

143

#include "ofApp.h"

//--------------------------------------------------------------

void ofApp::setup() {

ofSetFrameRate(30);

ofSetWindowTitle("openFrameworks");

ofBackground(239);

ofSetLineWidth(3);

ofEnableDepthTest();

this->frame.setMode(ofPrimitiveMode::OF_PRIMITIVE_LINES);

}

//--------------------------------------------------------------

void ofApp::update() {

ofSeedRandom(39);

this->face.clear();

this->frame.clear();

float phi_deg_step = 6;

float theta_deg_step = 6;

int R_start = (ofGetFrameNum() * 4) % 350 + 30;

for (int R_param = R_start; R_param < R_start + 350; R_param += 80) {

int R = R_param % 350 + 30;

float r = R * 0.13;

auto rotation_y = glm::rotate(glm::mat4(), ofMap(ofNoise(ofRandom(1000), ofGetFrameNum() * 0.001), 0, 1, PI * -5, PI * 5), glm::vec3(0, 1, 0));

auto rotation_x = glm::rotate(glm::mat4(), ofMap(ofNoise(ofRandom(1000), ofGetFrameNum() * 0.001), 0, 1, PI * -5, PI * 5), glm::vec3(1, 0, 0));

for (float phi_deg = 0; phi_deg < 360; phi_deg += phi_deg_step) {

for (float theta_deg = 0; theta_deg < 360; theta_deg += theta_deg_step) {

float threshold = ofMap(R, 30, 350, 0, 1);

auto noise_value = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg, phi_deg) * 0.005));

if (noise_value < threshold) { continue; }

auto noise_1 = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg - theta_deg_step, phi_deg) * 0.005));

auto noise_2 = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg, phi_deg + phi_deg_step) * 0.005));

auto noise_3 = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg, phi_deg - phi_deg_step) * 0.005));

auto noise_4 = ofNoise(glm::vec3(this->make_point(300, 80, theta_deg + theta_deg_step, phi_deg) * 0.005));

auto index = this->face.getNumVertices();

vector<glm::vec3> vertices;

vertices.push_back(glm::vec3(this->make_point(R, r, theta_deg - theta_deg_step * 0.5, phi_deg - phi_deg_step * 0.5)));

vertices.push_back(glm::vec3(this->make_point(R, r, theta_deg + theta_deg_step * 0.5, phi_deg - phi_deg_step * 0.5)));

vertices.push_back(glm::vec3(this->make_point(R, r, theta_deg - theta_deg_step * 0.5, phi_deg + phi_deg_step * 0.5)));

vertices.push_back(glm::vec3(this->make_point(R, r, theta_deg + theta_deg_step * 0.5, phi_deg + phi_deg_step * 0.5)));

for (auto& vertex : vertices) {

vertex = glm::vec4(vertex, 0) * rotation_y * rotation_x;

}

this->face.addVertices(vertices);

this->face.addIndex(index + 0); this->face.addIndex(index + 1); this->face.addIndex(index + 3);

this->face.addIndex(index + 0); this->face.addIndex(index + 3); this->face.addIndex(index + 2);

if (noise_1 < threshold) {

this->frame.addVertex(vertices[0]);

this->frame.addVertex(vertices[2]);

this->frame.addIndex(this->frame.getNumVertices() - 1);

this->frame.addIndex(this->frame.getNumVertices() - 2);

}

if (noise_2 < threshold) {

this->frame.addVertex(vertices[2]);

this->frame.addVertex(vertices[3]);

this->frame.addIndex(this->frame.getNumVertices() - 1);

this->frame.addIndex(this->frame.getNumVertices() - 2);

}

if (noise_3 < threshold) {

this->frame.addVertex(vertices[0]);

this->frame.addVertex(vertices[1]);

this->frame.addIndex(this->frame.getNumVertices() - 1);

this->frame.addIndex(this->frame.getNumVertices() - 2);

}

if (noise_4 < threshold) {

this->frame.addVertex(vertices[1]);

this->frame.addVertex(vertices[3]);

this->frame.addIndex(this->frame.getNumVertices() - 1);

this->frame.addIndex(this->frame.getNumVertices() - 2);

}

//--------------------------------------------------------------

void ofApp::draw() {

this->cam.begin();

ofSetColor(255);

this->frame.draw();

ofSetColor(39);

this->face.draw();

this->cam.end();

}

//--------------------------------------------------------------

glm::vec3 ofApp::make_point(float R, float r, float u, float v) {

// 数学デッサン教室描いて楽しむ数学たち　P.31

u *= DEG_TO_RAD;

v *= DEG_TO_RAD;

auto x = (R + r * cos(u)) * cos(v);

auto y = (R + r * cos(u)) * sin(v);

auto z = r * sin(u);

return glm::vec3(x, y, z);

}

//--------------------------------------------------------------

int main() {

ofSetupOpenGL(720, 720, OF_WINDOW);

ofRunApp(new ofApp());

}

[ Link ]
https://github.com/junkiyoshi/Insta20210308

2021年1月27日2021年1月27日

Shell and noise. Draw by openFrameworks

[ Video ]

[ About ]
Shell and noise.

[ Source ]

#pragma once
#include "ofMain.h"

class ofApp : public ofBaseApp {

public:
	void setup();
	void update();
	void draw();

	void keyPressed(int key) {};
	void keyReleased(int key) {};
	void mouseMoved(int x, int y) {};
	void mouseDragged(int x, int y, int button) {};
	void mousePressed(int x, int y, int button) {};
	void mouseReleased(int x, int y, int button) {};
	void windowResized(int w, int h) {};
	void dragEvent(ofDragInfo dragInfo) {};
	void gotMessage(ofMessage msg) {};

	glm::vec3 make_shell_point(float a, float b, float c, float h, float k, float u, float v, float w, int R);

	ofEasyCam cam;
	ofMesh face, line;
};

#pragma once

#include "ofMain.h"

class ofApp : public ofBaseApp {

public:

void setup();

void update();

void draw();

void keyPressed(int key) {};

void keyReleased(int key) {};

void mouseMoved(int x, int y) {};

void mouseDragged(int x, int y, int button) {};

void mousePressed(int x, int y, int button) {};

void mouseReleased(int x, int y, int button) {};

void windowResized(int w, int h) {};

void dragEvent(ofDragInfo dragInfo) {};

void gotMessage(ofMessage msg) {};

glm::vec3 make_shell_point(float a, float b, float c, float h, float k, float u, float v, float w, int R);

ofEasyCam cam;

ofMesh face, line;

};

#include "ofApp.h"

//--------------------------------------------------------------
void ofApp::setup() {

	ofSetFrameRate(60);
	ofSetWindowTitle("openFrameworks");

	ofBackground(239);
	ofSetLineWidth(1.5);
	ofEnableDepthTest();

	this->line.setMode(ofPrimitiveMode::OF_PRIMITIVE_LINES);
}

//--------------------------------------------------------------
void ofApp::update() {

	// Parameter reference by https://twitter.com/ru_sack/status/798308368007892992

	this->face.clear();
	this->line.clear();

	float a = 1.25;
	float b = 1.25;
	float c = 1.0;
	float h = 3.5;
	float k = 0;
	float w = 0.12;
	int R = 1;
	int scale = 40;

	float u_step = 0.3;
	int v_step = 15;

	int rotate_deg = ofGetFrameNum();

	for (int x = -300; x <= 300; x += 300) {

		for (int y = -300; y <= 300; y += 300) {

			rotate_deg += 40;
			float angle = rotate_deg * (float)DEG_TO_RAD;
			auto rotation = glm::rotate(glm::mat4(), angle, glm::vec3(0, 1, 0));

			for (float u = -35; u <= -1; u += u_step) {

				for (int v = 0; v <= 360; v += v_step) {

					auto noise_location = glm::vec2(x + cos(v * DEG_TO_RAD), y + sin(v * DEG_TO_RAD));
					if (ofNoise(glm::vec4(u * 0.005, noise_location * 0.5, ofGetFrameNum() * 0.01)) < 0.5) { continue; }

					vector<glm::vec3> vertices;
					vertices.push_back(this->make_shell_point(a, b, c, h, k, u + u_step * 0.5, (v + v_step * 0.5) * DEG_TO_RAD, w, R) * scale); // -4
					vertices.push_back(this->make_shell_point(a, b, c, h, k, u + u_step * 0.5, (v - v_step * 0.5) * DEG_TO_RAD, w, R) * scale); // -3
					vertices.push_back(this->make_shell_point(a, b, c, h, k, u - u_step * 0.5, (v - v_step * 0.5) * DEG_TO_RAD, w, R) * scale); // -2
					vertices.push_back(this->make_shell_point(a, b, c, h, k, u - u_step * 0.5, (v + v_step * 0.5) * DEG_TO_RAD, w, R) * scale); // -1

					for (auto& vertex : vertices) {

						vertex = glm::vec3(x, y, 0) + glm::vec4(vertex, 0) * rotation;
					}

					this->face.addVertices(vertices);

					this->face.addIndex(this->face.getNumVertices() - 1); this->face.addIndex(this->face.getNumVertices() - 2); this->face.addIndex(this->face.getNumVertices() - 3);
					this->face.addIndex(this->face.getNumVertices() - 1); this->face.addIndex(this->face.getNumVertices() - 3); this->face.addIndex(this->face.getNumVertices() - 4);

					this->line.addVertices(vertices);

					this->line.addIndex(this->line.getNumVertices() - 1); this->line.addIndex(this->line.getNumVertices() - 2);
					this->line.addIndex(this->line.getNumVertices() - 1); this->line.addIndex(this->line.getNumVertices() - 4);
					this->line.addIndex(this->line.getNumVertices() - 3); this->line.addIndex(this->line.getNumVertices() - 2);
					this->line.addIndex(this->line.getNumVertices() - 3); this->line.addIndex(this->line.getNumVertices() - 4);
				}
			}
		}
	}
}

//--------------------------------------------------------------
void ofApp::draw() {

	this->cam.begin();

	ofSetColor(239);
	this->face.draw();

	ofSetColor(39);
	this->line.draw();

	this->cam.end();
}

//--------------------------------------------------------------
glm::vec3 ofApp::make_shell_point(float a, float b, float c, float h, float k, float u, float v, float w, int R) {

	// 数学デッサン教室 描いて楽しむ数学たち 貝殻の形 shape of shell P.59

	float x = (h + a * cos(v)) * exp(w * u) * cos(c * u);
	float y = R * (h + a * cos(v)) * exp(w * u) * sin(c * u);
	float z = (k + b * sin(v)) * exp(w * u);

	return glm::vec3(x, y, z);
}

//--------------------------------------------------------------
int main() {

	ofSetupOpenGL(720, 720, OF_WINDOW);
	ofRunApp(new ofApp());
}

100

101

102

103

104

105

106

107

108

109

110

111

112

#include "ofApp.h"

//--------------------------------------------------------------

void ofApp::setup() {

ofSetFrameRate(60);

ofSetWindowTitle("openFrameworks");

ofBackground(239);

ofSetLineWidth(1.5);

ofEnableDepthTest();

this->line.setMode(ofPrimitiveMode::OF_PRIMITIVE_LINES);

}

//--------------------------------------------------------------

void ofApp::update() {

// Parameter reference by https://twitter.com/ru_sack/status/798308368007892992

this->face.clear();

this->line.clear();

float a = 1.25;

float b = 1.25;

float c = 1.0;

float h = 3.5;

float k = 0;

float w = 0.12;

int R = 1;

int scale = 40;

float u_step = 0.3;

int v_step = 15;

int rotate_deg = ofGetFrameNum();

for (int x = -300; x <= 300; x += 300) {

for (int y = -300; y <= 300; y += 300) {

rotate_deg += 40;

float angle = rotate_deg * (float)DEG_TO_RAD;

auto rotation = glm::rotate(glm::mat4(), angle, glm::vec3(0, 1, 0));

for (float u = -35; u <= -1; u += u_step) {

for (int v = 0; v <= 360; v += v_step) {

auto noise_location = glm::vec2(x + cos(v * DEG_TO_RAD), y + sin(v * DEG_TO_RAD));

if (ofNoise(glm::vec4(u * 0.005, noise_location * 0.5, ofGetFrameNum() * 0.01)) < 0.5) { continue; }

vector<glm::vec3> vertices;

vertices.push_back(this->make_shell_point(a, b, c, h, k, u + u_step * 0.5, (v + v_step * 0.5) * DEG_TO_RAD, w, R) * scale); // -4

vertices.push_back(this->make_shell_point(a, b, c, h, k, u + u_step * 0.5, (v - v_step * 0.5) * DEG_TO_RAD, w, R) * scale); // -3

vertices.push_back(this->make_shell_point(a, b, c, h, k, u - u_step * 0.5, (v - v_step * 0.5) * DEG_TO_RAD, w, R) * scale); // -2

vertices.push_back(this->make_shell_point(a, b, c, h, k, u - u_step * 0.5, (v + v_step * 0.5) * DEG_TO_RAD, w, R) * scale); // -1

for (auto& vertex : vertices) {

vertex = glm::vec3(x, y, 0) + glm::vec4(vertex, 0) * rotation;

}

this->face.addVertices(vertices);

this->face.addIndex(this->face.getNumVertices() - 1); this->face.addIndex(this->face.getNumVertices() - 2); this->face.addIndex(this->face.getNumVertices() - 3);

this->face.addIndex(this->face.getNumVertices() - 1); this->face.addIndex(this->face.getNumVertices() - 3); this->face.addIndex(this->face.getNumVertices() - 4);

this->line.addVertices(vertices);

this->line.addIndex(this->line.getNumVertices() - 1); this->line.addIndex(this->line.getNumVertices() - 2);

this->line.addIndex(this->line.getNumVertices() - 1); this->line.addIndex(this->line.getNumVertices() - 4);

this->line.addIndex(this->line.getNumVertices() - 3); this->line.addIndex(this->line.getNumVertices() - 2);

this->line.addIndex(this->line.getNumVertices() - 3); this->line.addIndex(this->line.getNumVertices() - 4);

}

//--------------------------------------------------------------

void ofApp::draw() {

this->cam.begin();

ofSetColor(239);

this->face.draw();

ofSetColor(39);

this->line.draw();

this->cam.end();

}

//--------------------------------------------------------------

glm::vec3 ofApp::make_shell_point(float a, float b, float c, float h, float k, float u, float v, float w, int R) {

// 数学デッサン教室描いて楽しむ数学たち貝殻の形 shape of shell P.59

float x = (h + a * cos(v)) * exp(w * u) * cos(c * u);

float y = R * (h + a * cos(v)) * exp(w * u) * sin(c * u);

float z = (k + b * sin(v)) * exp(w * u);

return glm::vec3(x, y, z);

}

//--------------------------------------------------------------

int main() {

ofSetupOpenGL(720, 720, OF_WINDOW);

ofRunApp(new ofApp());

}

[ Link ]
https://github.com/junkiyoshi/Insta20210127

2025年2月
月	火	水	木	金	土	日
					1	2
3	4	5	6	7	8	9
10	11	12	13	14	15	16
17	18	19	20	21	22	23
24	25	26	27	28