Uma Arquitetura de Suporte a Interações 3D ... - DCA - Unicamp

sampler sVPos; // Mapa de posições dos vértices
sampler sVN; // Mapa de normais aos vértices
sampler sAdjIdx; // Mapa de índices de adjacência
sampler sAdj; // Mapa de relações de adjacência
sampler sP1K1; // Mapa com direção e curvatura principal 1
sampler sP2K2; // Mapa com direção e curvatura principal 2
/////////////////////////////////////////////////////////////////////
// Realiza a decomposição LDL’ de uma matriz NxN simétrica positiva
// definida e resolve Ax=B.
/////////////////////////////////////////////////////////////////////
float3 ldltsv( float3x3 w, float3 m, float N ) {
// Decompõe matriz e sobrescreve matriz triangular inferior.
float3 d;
float2 v;
int i;
for( i = 0; i < N; i++ )
{
}
for( int k = 0; k < i; k++ ) v[k] = w[i][k] * d[k];
for( int j = i; j < N; j++ )
{
}
float sum = w[i][j];
for( int k = 0; k < i; k++ )
sum -= v[k] * w[j][k];
if( i == j )
else
d[i] = 1.0f / sum;
w[j][i] = sum;
// Resolve Ax=B após a decomposição LDL’
for( i = 0; i < N; i++ )
{
float sum = m[i];
for( int k = 0; k < i; k++ )
sum -= w[i][k] * m[k];
179