Changeset 8668

Timestamp:

06/20/11 09:04:59 (14 years ago)

Author:

Mathieu Morlighem

Message:

Better PatchToVec for P0 interpolation'

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 need to transform the vector a little bit first
-if size(Patch.value,2)==1,
-        vec=Patch.value(Patch.element);
-else
-        connectivity=sparse(Patch.index(:),1,1);
-        value       =sparse(Patch.index(:),1,Patch.value(:));
-        vec=full(value./connectivity);
-end
+switch(Patch.interpolation(1)),
+        case P0Enum,
+                vec(Patch.element)=Patch.value;
+        case P1Enum,
+                connectivity=sparse(Patch.index(:),1,1);
+                value       =sparse(Patch.index(:),1,Patch.value(:));
+                vec=full(value./connectivity);
+        otherwise,
+                error('interpolation not supported yet');
+        end