www.neuron.yale.edu

Code: Select all

Translating vecevent.mod into vecevent.c
Notice: VERBATIM blocks are not thread safe

Code: Select all

:  Vector stream of events

NEURON {
	ARTIFICIAL_CELL VecStim
}

ASSIGNED {
	index
	etime (ms)
	space
}

INITIAL {
	index = 0
	element()
	if (index > 0) {
		net_send(etime - t, 1)
	}
}

NET_RECEIVE (w) {
	if (flag == 1) {
		net_event(t)
		element()
		if (index > 0) {
			net_send(etime - t, 1)
		}
	}
}

VERBATIM
extern double* vector_vec();
extern int vector_capacity();
extern void* vector_arg();
ENDVERBATIM

PROCEDURE element() {
VERBATIM	
  { void* vv; int i, size; double* px;
	i = (int)index;
	if (i >= 0) {
		vv = *((void**)(&space));
		if (vv) {
			size = vector_capacity(vv);
			px = vector_vec(vv);
			if (i < size) {
				etime = px[i];
				index += 1.;
			}else{
				index = -1.;
			}
		}else{
			index = -1.;
		}
	}
  }
ENDVERBATIM
}

PROCEDURE play() {
VERBATIM
	void** vv;
	vv = (void**)(&space);
	*vv = (void*)0;
	if (ifarg(1)) {
		*vv = vector_arg(1);
	}
ENDVERBATIM
}

Code: Select all

NEURON {
	POINT_PROCESS Pr
	RANGE P, P0, random, f, tau_F, d1, tau_D1, d2, tau_D2, F, D1, D2, tlast
    THREADSAFE
    POINTER randObjPtr
}

PARAMETER {
    : the (1) is needed for the range limits to be effective
	P0 = 0.5 (1) < 0, 1 >               : basal release probability
    : these values are from Fig.3 in Varela et al. 1997
    f = 0.917 (1) < 0, 1e9 >            : additive facilitation per spike
    tau_F = 94 (ms) < 1e-9, 1e9 >       : rate of decay back to baseline following facilitation
    d1 = 0.416 (1) < 0, 1 >             : multiplicative fast depression per spike
    tau_D1 = 380 (ms) < 1e-9, 1e9 >     : rate of decay back to baseline following fast depression
    d2 = 0.975 (1) < 0, 1 >             : multiplicative slow depression per spike
    tau_D2 = 9200 (ms) < 1e-9, 1e9 >    : rate of decay back to baseline following slow depression
}

ASSIGNED {
	P				        : instantaneous release probability
    randObjPtr              : pointer to a hoc random number generator Random.uniform(0,1)
    random                  : individual instance of random number
    F                       : current level of facilitation
    D1                      : current level of fast depression
    D2                      : current level of slow depression
    tlast (ms)              : time of last spike
}

INITIAL {
	P = P0
    random = 1
    F = 1
    D1 = 1
    D2 = 1
}

NET_RECEIVE(weight) {
    INITIAL {
        tlast = t
    }
    F = 1 + (F-1)*exp(-(t - tlast)/tau_F)
    D1 = 1 - (1-D1)*exp(-(t - tlast)/tau_D1)
:    D2 = 1 - (1-D2)*exp(-(t - tlast)/tau_D2)
:    if (P0*F*D1*D2 > 1) {
    if (P0*F*D1 > 1) {
        P = 1
    } else {
:        P = P0*F*D1*D2
       P = P0*F*D1
    }
    random = randGen()
    if (random <= P) {
        net_event(t)
    }
    tlast = t
    F = F + f
    D1 = D1 * d1
:    D2 = D2 * d2
}

VERBATIM
double nrn_random_pick(void* r);
void* nrn_random_arg(int argpos);
ENDVERBATIM

FUNCTION randGen() {
VERBATIM
   if (_p_randObjPtr) {
      /*
      :Supports separate independent but reproducible streams for
      : each instance. However, the corresponding hoc Random
      : distribution MUST be set to Random.uniform(0,1)
      */
      _lrandGen = nrn_random_pick(_p_randObjPtr);
   }else{
      hoc_execerror("Random object ref not set correctly for randObjPtr"," only via hoc Random");
   }
ENDVERBATIM
}

PROCEDURE setRandObjRef() {
VERBATIM
   void** pv4 = (void**)(&_p_randObjPtr);
   if (ifarg(1)) {
      *pv4 = nrn_random_arg(1);
   }else{
      *pv4 = (void*)0;
   }
ENDVERBATIM
}

Code: Select all

vecevent.c:282:23: error: parameter named '_p' is missing
        Prop *_prop; double* _p; Datum* _ppvar; Datum* _thread;
                             ^
vecevent.c:282:34: error: parameter named '_ppvar' is missing
        Prop *_prop; double* _p; Datum* _ppvar; Datum* _thread;
                                        ^
vecevent.c:282:49: error: parameter named '_thread' is missing
        Prop *_prop; double* _p; Datum* _ppvar; Datum* _thread;
                                                       ^
vecevent.c:281:8: warning: type specifier missing, defaults to 'int' [-Wimplicit-int]
static _destructor(_prop)
~~~~~~ ^
vecevent.c:281:8: error: conflicting types for '_destructor'
vecevent.c:129:14: note: previous declaration is here
 static void _destructor(Prop*);
             ^
vecevent.c:284:2: error: use of undeclared identifier '_thread'
        _thread = (Datum*)0;
        ^
vecevent.c:285:2: error: use of undeclared identifier '_p'
        _p = prop->param; _ppvar = _prop->dparam;
        ^
vecevent.c:285:7: error: use of undeclared identifier 'prop'
        _p = prop->param; _ppvar = _prop->dparam;
             ^
vecevent.c:285:20: error: use of undeclared identifier '_ppvar'
        _p = prop->param; _ppvar = _prop->dparam;
                          ^
vecevent.c:290:21: error: use of undeclared identifier '_ppvar'
        void* vv = (void*)(_p_ptr);  
                           ^
vecevent.c:44:16: note: expanded from macro '_p_ptr'
#define _p_ptr  _ppvar[2]._pval
                ^
1 warning and 9 errors generated.