Changeset 8480 for trunk/pois/src/poisAddPenalty.c

Timestamp:

Aug 22, 2006, 2:48:04 PM (20 years ago)

Author:

Paul Price

Message:

Adding GPL stuff, for distribution

File:

• trunk/pois/src/poisAddPenalty.c (modified) (3 diffs)

trunk/pois/src/poisAddPenalty.c

@@ +1 @@
+// Copyright (C) 2006  Paul A. Price (price@ifa.hawaii.edu)
+//
+// This program is free software; you can redistribute it and/or modify
+// it under the terms of the GNU General Public License as published by
+// the Free Software Foundation; either version 2 of the License, or
+// (at your option) any later version.
+//
+// This program is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+
 #include <stdio.h>
 #include <assert.h>
@@ -4 +20 @@
 #include "pois.h"
 
-void poisAddPenalty(psImage *matrix,    // The matrix
-                    psVector *vector,   // The vector
-                    float penalty,      // Penalty value
-                    const psArray *kernels, // Kernel basis functions
-                    const poisConfig *config // Configuration
+void poisAddPenalty(psImage *matrix,    // The matrix
+                    psVector *vector,   // The vector
+                    float penalty,      // Penalty value
+                    const psArray *kernels, // Kernel basis functions
+                    const poisConfig *config // Configuration
     )
 {
@@ -25 +41 @@
 
     for (int j = 0; j < numKernels; j++) {
-        poisKernelBasis *kernel = kernels->data[j];
-        if (kernel->xOrder == 0 && kernel->yOrder == 0) {
-            for (int k = j + 1; k < numKernels; k++) {
-                poisKernelBasis *compare = kernels->data[k];
-                if (compare->xOrder == 0 && compare->yOrder == 0 &&
-                    (compare->u == kernel->u + 1 || compare->u == kernel->u - 1 ||
-                     compare->v == kernel->v + 1 || compare->v == kernel->v - 1)) {
-                    matrix->data.F64[j][k] -= penalty;
-                    matrix->data.F64[k][j] -= penalty;
-                }
-            }
-        }
-            
-        if (abs(kernel->u) == config->xKernel) {
-            if (abs(kernel->v) == config->yKernel) {
-                vector->data.F64[j] -= 3.0 * penalty;
-            } else {
-                vector->data.F64[j] -= 5.0 * penalty;
-            }
-        } else {
-            if (abs(kernel->v) == config->yKernel) {
-                vector->data.F64[j] -= 5.0 * penalty;
-            } else {
-                vector->data.F64[j] -= 8.0 * penalty;
-            }
-        }
-        
+        poisKernelBasis *kernel = kernels->data[j];
+        if (kernel->xOrder == 0 && kernel->yOrder == 0) {
+            for (int k = j + 1; k < numKernels; k++) {
+                poisKernelBasis *compare = kernels->data[k];
+                if (compare->xOrder == 0 && compare->yOrder == 0 &&
+                    (compare->u == kernel->u + 1 || compare->u == kernel->u - 1 ||
+                     compare->v == kernel->v + 1 || compare->v == kernel->v - 1)) {
+                    matrix->data.F64[j][k] -= penalty;
+                    matrix->data.F64[k][j] -= penalty;
+                }
+            }
+        }
+
+        if (abs(kernel->u) == config->xKernel) {
+            if (abs(kernel->v) == config->yKernel) {
+                vector->data.F64[j] -= 3.0 * penalty;
+            } else {
+                vector->data.F64[j] -= 5.0 * penalty;
+            }
+        } else {
+            if (abs(kernel->v) == config->yKernel) {
+                vector->data.F64[j] -= 5.0 * penalty;
+            } else {
+                vector->data.F64[j] -= 8.0 * penalty;
+            }
+        }
+
     }