创建一个晶体生长效果着色器，模拟晶体形成的动态过程，包含晶格结构和折射效果，提供三种晶体类型可选。

uniform float uTime;
uniform float uGrowthSpeed;
uniform float uCrystalPurity;
uniform float uRefractionIndex;

void main() {
    vUv = uv;
    vPosition = position;
    vNormal = normal;

    float growth = sin(position.x * 3.0 + uTime * uGrowthSpeed) *
                  cos(position.y * 4.0 + uTime * uGrowthSpeed * 0.8) *
                  sin(position.z * 3.0 + uTime * uGrowthSpeed * 0.6);

    float purity = sin(position.x * 5.0 + uTime * uCrystalPurity * 2.0) *
                  cos(position.z * 5.0 + uTime * uCrystalPurity * 2.0);

    float refraction = sin(position.x * uRefractionIndex * 6.0) *
                      cos(position.y * uRefractionIndex * 8.0);

    vec3 newPosition = position + normal * (growth + purity * 0.5 + refraction * 0.3) * 0.15;

    gl_Position = projectionMatrix * modelViewMatrix * vec4(newPosition, 1.0);
}

float crystalGrowth(vec2 uv, float speed) {
    float growthField = 0.0;
    growthField += sin(uv.x * 4.0 + uTime * 0.3) *
                  cos(uv.y * 6.0 + uTime * 0.5);
    growthField += sin(uv.x * 6.0 + uTime * 0.7) *
                  cos(uv.y * 8.0 + uTime * 0.9) * 0.7;
    growthField += sin(uv.x * 8.0 + uTime * 1.1) *
                  cos(uv.y * 10.0 + uTime * 1.3) * 0.4;
    return growthField * speed;
}

float crystalLattice(vec2 uv, float density) {
    float latticeField = 0.0;
    for(int i = 0; i < 20; i++) {
        vec2 latticePos = vec2(
            noise(vec2(float(i) * 0.123, 0.456)),
            noise(vec2(float(i) * 0.456, 0.789))
        );
        float dist = distance(uv, latticePos);
        float latticeSize = noise(vec2(float(i) * 0.789, 0.123)) * 0.04 + 0.02;
        float latticeLife = mod(uTime * 1.5 + float(i) * 8.0, 6.28);
        float latticeIntensity = sin(latticeLife) * 0.7 + 0.3;
        latticeField += exp(-dist * 45.0) * latticeIntensity * density;
    }
    return latticeField;
}