Changeset 9673

Timestamp:

09/07/11 16:21:42 (14 years ago)

Author:

Mathieu Morlighem

Message:

Simpler PatchToVec

File:

• issm/trunk/src/m/utils/Interp/PatchToVec.m (modified) (1 diff)

issm/trunk/src/m/utils/Interp/PatchToVec.m

@@ -6 +6 @@
 
 %if the patch is P0: we have element values, return an element vector
-switch(Patch.interpolation(1)),
-        case P0Enum,
+switch(size(Patch.value,2)),
+        case 1,
                 vec(Patch.element)=Patch.value;
-                %connectivity=sparse(Patch.index(:),1,1);
-                %value       =sparse(Patch.index(:),1,Patch.value(:));
-                %vec=full(value./connectivity);
-        case P1Enum,
+        case 3,
+                connectivity=sparse(Patch.index(:),1,1);
+                value       =sparse(Patch.index(:),1,Patch.value(:));
+                vec=full(value./connectivity);
+        case 6,
                 connectivity=sparse(Patch.index(:),1,1);
                 value       =sparse(Patch.index(:),1,Patch.value(:));