source: data/branches/Shader_HS18/programs/Example/GLSL150/SampleFieldVS.glsl @ 12091

#version 150

// Ogre port of Nvidia's IsoSurf.cg file
// Modified code follows. See http://developer.download.nvidia.com/SDK/10/opengl/samples.html for original
//
// Cg port of Yury Uralsky's metaball FX shader
//
// Authors: Simon Green and Yury Urlasky
// Email: sdkfeedback@nvidia.com
//
// Copyright (c) NVIDIA Corporation. All rights reserved.
////////////////////////////////////////////////////////////////////////////////////////////////////

uniform float IsoValue = 1.0;
uniform mat4 WorldViewProj;
uniform mat4 origWorldViewIT;

// Size of the sampling grid
const ivec3 SizeMask = ivec3( 63, 63, 63 );
const ivec3 SizeShift = ivec3( 0, 6, 12 );
uniform vec4 Metaballs[2] = vec4[](vec4( -0.5, 0, 0, 0.2 ), vec4( 0.6, 0, 0, 0.1 ));

in vec4 vertex;

out VertexData {
    vec3 N;
    vec2 Field;
} VertexOut;

// Metaball function
// Returns metaball function value in .w and its gradient in .xyz

vec4 Metaball(vec3 Pos, vec3 Center, float RadiusSq)
{
        const float epsilon = 0.001;

        vec3 Dist = Pos - Center;
        float InvDistSq = 1 / (dot(Dist, Dist) + epsilon);

        vec4 o;
        o.xyz = -2 * RadiusSq * InvDistSq * InvDistSq * Dist;
        o.w = RadiusSq * InvDistSq;
        return o;
}

void main()
{
        // Sum up contributions from all metaballs
/*      vec4 Field;
        for (int i = 0; i<2; i++) 
        {
                Field += Metaball(vertex.xyz, Metaballs[i].xyz, Metaballs[i].w);
        }
        
        mat3 WorldViewIT = mat3(origWorldViewIT[0].xyz, origWorldViewIT[1].xyz, origWorldViewIT[2].xyz);
        
        // Transform position and normals
        gl_Position = WorldViewProj * vertex;
        VertexOut.N = WorldViewIT * Field.xyz;  // we want normals in world space
        VertexOut.Field.x = Field.w;
        
        // Generate in-out flags
        VertexOut.Field.y = (Field.w < 1.0) ? 1 : 0;
*/
        vec3 Pos;

        // Generate sampling point position based on its index
        Pos.x = float((gl_VertexID >> SizeShift.x) & SizeMask.x) / (SizeMask.x + 1);
        Pos.y = float((gl_VertexID >> SizeShift.y) & SizeMask.y) / (SizeMask.y + 1);
        Pos.z = float((gl_VertexID >> SizeShift.z) & SizeMask.z) / (SizeMask.z + 1);
        Pos = Pos*2 - 1;

        // Sum up contributions from all metaballs
        vec4 Field;
        for (int i = 0; i<2; i++)
                Field += Metaball(Pos, Metaballs[i].xyz, Metaballs[i].w);

        // Transform position and normals
        gl_Position = WorldViewProj * vec4(Pos, 1);
        VertexOut.N = mat3(origWorldViewIT) * Field.xyz;        // we want normals in world space
        VertexOut.Field.x = Field.w;

        // Generate in-out flags
        VertexOut.Field.y = (Field.w < IsoValue) ? 1 : 0;

//      o.Color = (Field*0.5+0.5) * (Field.w / 10.0);
}