Changeset 895

Timestamp:

06/11/09 08:51:08 (16 years ago)

Author:

seroussi

Message:

compute B on element before loop on gaussian points, as in cielo

File:

• issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m (modified) (2 diffs)

issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m

@@ -278 +278 @@
         end
         if(temperature_average_is_present) temperature_average_list(i)=temperature_average_param(dof);end;
-end
+
+end
+
+%Update material if temperature is provided.
+if temperature_average_is_present,
+        temperature_average=1/3*([1 1 1]*temperature_average_list);
+        matice.B=paterson(temperature_average);
+end
+
 if (element_debug & triaelem.id==element_debugid),
         disp(sprintf('El id %i TriaElemnet input list before gaussian loop: \n',element_debugid)); 
@@ -381 +389 @@
         %disp(sprintf('Epsilon: %f %f %f\n',epsilon(1),epsilon(2),epsilon(3)));
         
-        %Update material if temperature is provided.
-        if temperature_average_is_present,
-                temperature_average=GetParameterValue(triaelem,temperature_average_list,gauss_l1l2l3);
-                matice.B=paterson(temperature_average);
-        end
-
         %Update material if flow law is specified. This will erase the previous change 
         %on B when temperature is provided.