用C语言50行代码绘制一朵玫瑰花

先看看效果图

代码在这里奥

#include <stdio.h>
#include <math.h>const int max_iterations = 128;
const float stop_threshold = 0.01f;
const float grad_step = 0.01f;
const float clip_far = 10.0f;const float PI = 3.14159265359f;
const float PI2 = 6.28318530718f;
const float DEG_TO_RAD = PI / 180.0f;typedef struct { float x, y; } vec2;
typedef struct { float x, y, z; } vec3;
typedef struct { float m[9]; } mat3;const vec3 light_pos = { 20.0f, 50.0f, 20.0f };float min(float a, float b) { return a < b ? a : b; }
float max(float a, float b) { return a > b ? a : b; }
float clamp(float f, float a, float b) { return max(min(f, b), a); }
vec2 make2(float x, float y) { vec2 r = { x, y }; return r; }
vec2 add2(vec2 a, vec2 b) { vec2 r = { a.x + b.x, a.y + b.y }; return r; }
vec2 sub2(vec2 a, vec2 b) { vec2 r = { a.x - b.x, a.y - b.y }; return r; }
float dot2(vec2 a, vec2 b) { return a.x * b.x + a.y * b.y; }
float length2(vec2 v) { return sqrt(dot2(v, v)); }
vec3 make3(float x, float y, float z) { vec3 r = { x, y, z }; return r; }
vec3 add3(vec3 a, vec3 b) { vec3 r = { a.x + b.x, a.y + b.y, a.z + b.z }; return r; }
vec3 sub3(vec3 a, vec3 b) { vec3 r = { a.x - b.x, a.y - b.y, a.z - b.z }; return r; }
vec3 mul3(vec3 a, vec3 b) { vec3 r = { a.x * b.x, a.y * b.y, a.z * b.z }; return r; }
vec3 scale3(vec3 v, float s) { vec3 r = { v.x * s, v.y * s, v.z * s }; return r; }
float dot3(vec3 a, vec3 b) { return a.x * b.x + a.y * b.y + a.z * b.z; }
float length3(vec3 v) { return sqrt(dot3(v, v)); }
vec3 normalize3(vec3 v) { return scale3(v, 1.0f / length3(v)); }
vec3 mul(mat3 m, vec3 v) {return make3(m.m[0] * v.x + m.m[3] * v.y + m.m[6] * v.z,m.m[1] * v.x + m.m[4] * v.y + m.m[7] * v.z,m.m[2] * v.x + m.m[5] * v.y + m.m[8] * v.z);
}mat3 rotationXY(float x, float y) {vec2 c = { cos(x), cos(y) }, s = { sin(x), sin(y) };mat3 m = {c.y      , 0.0f, -s.y,s.y * s.x,  c.x,  c.y * s.x,s.y * c.x, -s.x,  c.y * c.x};return m;
}float opI(float d1, float d2) { return max(d1, d2); }
float opU(float d1, float d2) { return min(d1, d2); }
float opS(float d1, float d2) { return max(-d1, d2); }float sdPetal(vec3 p, float s) {p = add3(mul3(p, make3(0.8f, 1.5f, 0.8f)), make3(0.1f, 0.0f, 0.0f));vec2 q = make2(length2(make2(p.x, p.z)), p.y);float lower = length2(q) - 1.0f;lower = opS(length2(q) - 0.97f, lower);lower = opI(lower, q.y);float upper = length2(sub2(q, make2(s, 0.0f))) + 1.0f - s;upper = opS(upper, length2(sub2(q, make2(s, 0.0f))) + 0.97f - s);upper = opI(upper, -q.y);upper = opI(upper, q.x - 2.0f);float region = length3(sub3(p, make3(1.0f, 0.0f, 0.0f))) - 1.0f;return opI(opU(upper, lower), region);
}float map(vec3 p) {float d = 1000.0f, s = 2.0f;mat3 r = rotationXY(0.1f, PI2 * 0.618034f);r.m[0] *= 1.08f;  r.m[1] *= 1.08f;  r.m[2] *= 1.08f;r.m[3] *= 0.995f; r.m[4] *= 0.995f; r.m[5] *= 0.995f;r.m[6] *= 1.08f;  r.m[7] *= 1.08f;  r.m[8] *= 1.08f;for (int i = 0; i < 21; i++) {d = opU(d, sdPetal(p, s));p = mul(r, p);p = add3(p, make3(0.0, -0.02, 0.0));s *= 1.05f;}return d;
}vec3 gradient(vec3 pos) {const vec3 dx = { grad_step, 0.0, 0.0 };const vec3 dy = { 0.0, grad_step, 0.0 };const vec3 dz = { 0.0, 0.0, grad_step };return normalize3(make3(map(add3(pos, dx)) - map(sub3(pos, dx)),map(add3(pos, dy)) - map(sub3(pos, dy)),map(add3(pos, dz)) - map(sub3(pos, dz))));
}float ray_marching(vec3 origin, vec3 dir, float start, float end) {float depth = start;for (int i = 0; i < max_iterations; i++) {float dist = map(add3(origin, scale3(dir, depth)));if (dist < stop_threshold)return depth;depth += dist * 0.3;if (depth >= end)return end;}return end;
}float shading(vec3 v, vec3 n, vec3 eye) {vec3 ev = normalize3(sub3(v, eye));vec3 vl = normalize3(sub3(light_pos, v));float diffuse = dot3(vl, n) * 0.5f + 0.5f;vec3 h = normalize3(sub3(vl, ev));float rim = pow(1.0f - max(-dot3(n, ev), 0.0f), 2.0f) * 0.15f;float ao = clamp(v.y * 0.5f + 0.5f, 0.0f, 1.0f);return (diffuse + rim) * ao;
}vec3 ray_dir(float fov, vec2 pos) {vec3 r = { pos.x, pos.y, -tan((90.0f - fov * 0.5f) * DEG_TO_RAD) };return normalize3(r);
}float f(vec2 fragCoord) {vec3 dir = ray_dir(45.0f, fragCoord);vec3 eye = { 0.0f, 0.0f, 4.5f };mat3 rot = rotationXY(-1.0f, 1.0f);dir = mul(rot, dir);eye = mul(rot, eye);float depth = ray_marching(eye, dir, 0.0f, clip_far);vec3 pos = add3(eye, scale3(dir, depth));if (depth >= clip_far)return 0.0f;elsereturn shading(pos, gradient(pos), eye);
}int main() {puts("\033[91m");for (int y = 0; y < 80; y++) {for (int x = 0; x < 160; x++)putchar("  .,-:;+=*#@"[(int)(f(make2((x / 160.0f - 0.5f) * 2.0f, (y / 80.0f - 0.5f) * -2.0f)) * 12.0f)]);putchar('\n');}
}

大家需要源码的可以来小编的学习群【669325409】领取！

用C语言50行代码绘制一朵玫瑰花相关推荐

教你如何用91行代码实现一朵玫瑰花的绘制
教你如何用91行代码实现一朵玫瑰花的绘制. 文章目录教你如何用91行代码实现一朵玫瑰花的绘制. 主体: 初始位置的设定绘制花朵形状: 绘制花枝形状绘制一个绿色叶子要用到的工具: Pycharm ...
2022年第一篇C语言50行代码制作爱心电子相册
文章目录介绍效果展示源代码素材总结介绍去年看到过,今天突然想到,记录下来.给初学者留个好印象. 关于开发工具------浪漫烟花文章有介绍效果展示动态图源代码 #include&l ...
如何用50行代码构建情感分类器
选自Toward Data Science,作者:Rohith Gandhi,机器之心编译. 本文介绍了如何构建情感分类器,从介绍自然语言处理开始,一步一步讲述构建过程. 自然语言处理简介语言把人类 ...
C语言100行代码实现推箱子
1 C语言100行代码实现推箱子 1.1 概述 C语言是很好入门编程的一个语言,它拥有着很好的移植性,基本上所有的平台都支持C语言编程.有些C语言基础的你,是不是也很想做一个项目来检验一下自己的学习 ...
python爬虫实战：利用scrapy，短短50行代码下载整站短视频
近日,有朋友向我求助一件小事儿,他在一个短视频app上看到一个好玩儿的段子,想下载下来,可死活找不到下载的方法.这忙我得帮,少不得就抓包分析了一下这个app,找到了视频的下载链接,帮他解决了这个小问题 ...
python实现50行代码_50行代码实现python计算器主要功能
实现功能:计算带有括号和四则运算的式子 3*( 4+ 50 )-(( 100 + 40 )*5/2- 3*2* 2/4+9)*((( 3 + 4)-4)-4) 基本思路:使用正则表达式提取出每一层小括 ...
pyquery获取不到网页完整源代码_爬虫神器之PyQuery实用教程（二），50行代码爬取穷游网...
爬虫神器之PyQuery实用教程(二),50行代码爬取穷游网前言上篇文章 PyQuery (一) 回顾.今天来介绍具体 PyQuery 的使用方法. 穷游网目标与分析开始之前,按照之前的套路一步 ...
c语言300行代码大作业,C语言300行代码
<C语言300行代码>由会员分享,可在线阅读,更多相关<C语言300行代码(3页珍藏版)>请在人人文库网上搜索. 1.include #include #include #in ...
【PyTorch】50行代码实现GAN——PyTorch
本文来源于PyTorch中文网. 一直想了解GAN到底是个什么东西,却一直没能腾出时间来认真研究,前几日正好搜到一篇关于PyTorch实现GAN训练的文章,特将学习记录如下,本文主要包含两个部分:GA ...

用C语言50行代码绘制一朵玫瑰花

用C语言50行代码绘制一朵玫瑰花相关推荐

最新文章

热门文章