Lapicque.cc

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) NumeriCor GmbH
// SPDX-License-Identifier: LicenseRef-APL-1.1
 
#include "Lapicque.h"
#include <iostream>
 
Lapicque:: Lapicque( opencarp::timer_manager *to, float pmin, float pmax, int np, float tol, float delay ):
          _tmo(to),_delay(delay),_tol(tol)
{
  _widths.resize(np);
  _thresh.resize(np);
  float dlog = _widths.size()>1 ? (log(pmax)-log(pmin))/(_widths.size()-1) : 0;
  for( int i=0; i<_widths.size(); i++ ) {
    _widths[i] =  pmin*exp(dlog*i);
  }
  _ap = new Bisect( 0, INFINITY, tol );
}
 
 
bool
Lapicque::check( bool fire, float &curr, int doppind, int stimind ) 
{
  if( _ap->update( curr, fire ) ){
     std::cout << _widths[_currw] << "   " << curr << std::endl;
    _thresh[_currw++] = curr;
    if( _currw == _thresh.size() )
      return true;
    else {
      delete _ap;
      _ap = new Bisect( 0, INFINITY, _tol );
    }
  }
  delete _tmo->timers[doppind];
  double ttime = _tmo->time+_widths[_currw]+_delay;
  _tmo->initialize_singlestep_timer( ttime, 1, doppind, "dopple" );
  delete _tmo->timers[stimind];
  _tmo->initialize_singlestep_timer( _tmo->time, _widths[_currw], stimind, "STIM_TIMER" );
  return false;
}
 
 
void
Lapicque:: print( std::ostream &ofs )
{
  ofs << "# stimDuration(ms) stimStrength(pA/pF)" << std::endl;
  for( int i=0; i<_widths.size(); i++ )
    ofs << _widths[i] << " " << _thresh[i] << std::endl;
}
 
 
bool
Lapicque::init( int doppind, int stimind )
{
  if( _init ) return true;
  _init = true;
  double ttime = _tmo->time+_widths[0]+_delay;
  delete _tmo->timers[doppind];
  _tmo->initialize_singlestep_timer( ttime, 0, doppind, "dopple" );
  delete _tmo->timers[stimind];
  _tmo->initialize_singlestep_timer( _tmo->time, _widths[0], stimind, "STIM_TIMER" );
  return false;
}