Changeset 8659 for issm/trunk

Timestamp:

06/17/11 16:18:24 (14 years ago)

Author:

Mathieu Morlighem

Message:

better PatchToVec for P0

File:

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

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

@@ -7 +7 @@
 %if the patch is P0, we need to transform the vector a little bit first
 if size(Patch.value,2)==1,
-        Patch.value=[Patch.value Patch.value Patch.value];
+        vec=Patch.value(Patch.element);
+else
+        connectivity=sparse(Patch.index(:),1,1);
+        value       =sparse(Patch.index(:),1,Patch.value(:));
+        vec=full(value./connectivity);
 end
-connectivity=sparse(Patch.index(:),1,1);
-value       =sparse(Patch.index(:),1,Patch.value(:));
-vec=full(value./connectivity);