changement locaux : resistance / calcul

ik1ap.cpp

@@ -262,7 +262,7 @@ class IK1AP01
 					case '0': case '1': case '2': case '3': case '4': case '5':
 					case '6': case '7': case '8': case '9': 
 						{
-						double new_ik1dens = ((double)(ch - '0'))*0.2 + 0.2;
+						double new_ik1dens = ((double)(ch - '0'))*0.4 + 0.4;
 						if( ik1dens != new_ik1dens )
 							{
 							ik1dens = new_ik1dens ;
@@ -499,19 +499,19 @@ class IK1AP01
 						value_between(
 							(clock_now % MILLISECONDS_TO_SECONDS(Stimperiod_ms)) ,
 							MILLISECONDS_TO_SECONDS((Stimperiod_ms - 20)),
-							MILLISECONDS_TO_SECONDS((Stimperiod_ms - 19))
+							MILLISECONDS_TO_SECONDS((Stimperiod_ms - 19.7))
 							)
 						)
 						{
 
-						injected_current = -25 * membrane_cap;
+						injected_current = -5 * membrane_cap;
 						
 						
 						if( membranePotentialValues.size()< ARRAY_SIZE &&
 							value_between(
 							(clock_now % MILLISECONDS_TO_SECONDS(Stimperiod_ms)) ,
-							MILLISECONDS_TO_SECONDS(((double)Stimperiod_ms - 19.1)),
-							MILLISECONDS_TO_SECONDS((Stimperiod_ms - 19))
+							MILLISECONDS_TO_SECONDS(((double)Stimperiod_ms - 19.9)),
+							MILLISECONDS_TO_SECONDS((Stimperiod_ms - 19.8))
 							))
 							{
 							//fprintf(stderr,"membrane_potential: %f\n",membrane_potential);
@@ -527,8 +527,8 @@ class IK1AP01
 							controlMembranePotentialValues.size() < ARRAY_SIZE &&
 							this->value_between(
 								(clock_now % MILLISECONDS_TO_SECONDS(Stimperiod_ms)) ,
-								MILLISECONDS_TO_SECONDS(((double)Stimperiod_ms - 21.5)),
-								MILLISECONDS_TO_SECONDS((Stimperiod_ms - 21))
+								MILLISECONDS_TO_SECONDS(((double)Stimperiod_ms - 20.2)),
+								MILLISECONDS_TO_SECONDS((Stimperiod_ms - 20.1))
 								)
 							)
 							{
@@ -546,7 +546,7 @@ class IK1AP01
 								double avg1 = average(membranePotentialValues);
 								double avg2 = average(controlMembranePotentialValues);	
 								
-								resistance = -1 *  (avg1-avg2)/(25.0 * membrane_cap) * 1E3 ;
+								resistance = (0-(avg1-avg2))/(5.0)*1000;
 								membranePotentialValues.clear();
 								controlMembranePotentialValues.clear();
 								//fprintf(stderr,"CLEAR!!");