#ifndef BOT_ATYPE_LTYPE_RTYPE_OPNAME
#define BOT_ATYPE_LTYPE_RTYPE_OPNAME
IA_BinaryOperationsTable<IA_BaseImage<IA_Point<int>, ATYPE>*,
	IA_BaseImage<IA_Point<int>, LTYPE>,
	IA_BaseImage<IA_Point<int>, RTYPE> >
	*IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl;

void IA_fill_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl();

#endif // BOT_ATYPE_LTYPE_RTYPE_OPNAME
IA_Image<IA_Point<int>, ATYPE>
OP ( const IA_Image<IA_Point<int>, LTYPE> &lhs,
		RTYPE rhs)
{
    if (!IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl)
	IA_fill_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl();

    IA_BaseImage<IA_Point<int>, ATYPE>*(*f)( IA_BaseImage<IA_Point<int>, LTYPE> &,
		IA_BaseImage<IA_Point<int>, RTYPE> &);
    f = IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->lookup_operation
	(lhs.type(), IA_ConstI<IA_Point<int>, RTYPE>::s_type());
    IA_BaseImage<IA_Point<int>,LTYPE>	&basel = *IA_FBI<IA_Point<int>,IA_Point<int>, LTYPE,LTYPE,LTYPE>
	    ::extract_baseptr(lhs);
    IA_ConstI<IA_Point<int>,LTYPE>	baser(lhs.domain(), rhs);
    return IA_Image<IA_Point<int>,ATYPE>
	(IA_FBI<IA_Point<int>,IA_Point<int>, ATYPE,ATYPE,ATYPE>
	 ::make_Image(	 f(basel,baser)));
}

//
//
//

#ifndef vec_scalar_ATYPE_LTYPE_RTYPE_OPNAME
#define vec_scalar_ATYPE_LTYPE_RTYPE_OPNAME
static IA_BaseImage<IA_Point<int>, ATYPE> *
vec_scalar_bo_OPNAME(IA_BaseImage<IA_Point<int>, LTYPE> &lhs_,
	IA_BaseImage<IA_Point<int>, RTYPE> &rhs_)
{
    const IA_Set<IA_Point<int> > domain = lhs_.domain();
    const unsigned sz = domain.card();
    ATYPE *const vec = new ATYPE[sz];

    {
	ATYPE *ds = vec;
	IA_VectorI<IA_Point<int>, LTYPE> *lhs
		= (IA_VectorI<IA_Point<int>, LTYPE> *)&lhs_;
	const LTYPE *ls = lhs->vec;
	IA_ConstI<IA_Point<int>, RTYPE> *rhs
		= (IA_ConstI<IA_Point<int>, RTYPE> *)&rhs_;
	for (unsigned i=0; i<sz; i++) {
	    *(ds++) = ( (*ls) SUBOP (rhs->value) );
	    ls++;
	}
    }
    return new IA_VectorI<IA_Point<int>, ATYPE>(domain, vec, sz, 1);
}
#endif // vec_scalar_ATYPE_LTYPE_RTYPE_OPNAME

#ifndef scalar_scalar_ATYPE_LTYPE_RTYPE_OPNAME
#define scalar_scalar_ATYPE_LTYPE_RTYPE_OPNAME
static IA_BaseImage<IA_Point<int>, ATYPE> *
scalar_scalar_bo_OPNAME(IA_BaseImage<IA_Point<int>, LTYPE> &lhs_,
	IA_BaseImage<IA_Point<int>, RTYPE> &rhs_)
{
    const IA_Set<IA_Point<int> > domain = lhs_.domain();
	IA_ConstI<IA_Point<int>, LTYPE> *lhs
		= (IA_ConstI<IA_Point<int>, LTYPE> *)&lhs_;
	IA_ConstI<IA_Point<int>, RTYPE> *rhs
		= (IA_ConstI<IA_Point<int>, RTYPE> *)&rhs_;
    ATYPE	aval = ( lhs->value SUBOP rhs->value );
    return new IA_ConstI<IA_Point<int>, ATYPE>(domain, aval);
}
#endif // scalar_scalar_ATYPE_LTYPE_RTYPE_OPNAME

#ifndef scalar_iter_ATYPE_LTYPE_RTYPE_OPNAME
#define scalar_iter_ATYPE_LTYPE_RTYPE_OPNAME
static IA_BaseImage<IA_Point<int>, ATYPE> *
scalar_iter_bo_OPNAME(IA_BaseImage<IA_Point<int>, RTYPE> &lhs_,
	IA_BaseImage<IA_Point<int>, LTYPE> &rhs_)
{
    const IA_Set<IA_Point<int> > domain = lhs_.domain();
    const unsigned sz = domain.card();
    ATYPE *const vec = new ATYPE[sz];

    {
	ATYPE *ds = vec;
	IA_ConstI<IA_Point<int>, LTYPE> *lhs
		= (IA_ConstI<IA_Point<int>, LTYPE> *)&lhs_;
	RTYPE rtmp;
	IA_BaseIVIter<IA_Point<int>, RTYPE>	*riter = rhs_.value_iterator();
	for (unsigned i=0; i<sz; i++) {
	    riter->next(rtmp);
	    *(ds++) = ( (lhs->value) SUBOP (rtmp) );
	}
	delete riter;
    }
    return new IA_VectorI<IA_Point<int>, ATYPE>(domain, vec, sz, 1);
}
#endif // scalar_iter_ATYPE_LTYPE_RTYPE_OPNAME

#ifndef iter_scalar_ATYPE_LTYPE_RTYPE_OPNAME
#define iter_scalar_ATYPE_LTYPE_RTYPE_OPNAME
static IA_BaseImage<IA_Point<int>, ATYPE> *
iter_scalar_bo_OPNAME(IA_BaseImage<IA_Point<int>, LTYPE> &lhs_,
	IA_BaseImage<IA_Point<int>, RTYPE> &rhs_)
{
    return scalar_iter_bo_OPNAME(rhs_, lhs_);
}
#endif // iter_scalar_ATYPE_LTYPE_RTYPE_OPNAME

#ifndef FILL_BOT_ATYPE_LTYPE_RTYPE_OPNAME
#define FILL_BOT_ATYPE_LTYPE_RTYPE_OPNAME
void IA_fill_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl()
{
    if (IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl)
	return;

    IA_BinaryOperationsTable<IA_BaseImage<IA_Point<int>, ATYPE>*,
	IA_BaseImage<IA_Point<int>, LTYPE>,
	IA_BaseImage<IA_Point<int>, RTYPE> >::function f;
    f = iter_iter_bo_OPNAME;

    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl = new IA_BinaryOperationsTable<IA_BaseImage<IA_Point<int>, ATYPE>*,
	IA_BaseImage<IA_Point<int>, LTYPE>,
	IA_BaseImage<IA_Point<int>, RTYPE> >(iter_iter_bo_OPNAME);
    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->add_operation
	(IA_VectorI<IA_Point<int>, LTYPE>::s_type(),
	 IA_VectorI<IA_Point<int>, RTYPE>::s_type(),
	 vec_vec_bo_OPNAME);
    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->add_operation
	(IA_VectorI<IA_Point<int>, LTYPE>::s_type(),
	 IA_ConstI<IA_Point<int>, RTYPE>::s_type(),
	 vec_scalar_bo_OPNAME);
    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->add_operation
	(IA_VectorI<IA_Point<int>, LTYPE>::s_type(),
	 0,
	 vec_iter_bo_OPNAME);

    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->add_operation
	(IA_ConstI<IA_Point<int>, LTYPE>::s_type(),
	 IA_VectorI<IA_Point<int>, RTYPE>::s_type(),
	 scalar_vec_bo_OPNAME);
    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->add_operation
	(IA_ConstI<IA_Point<int>, LTYPE>::s_type(),
	 IA_ConstI<IA_Point<int>, RTYPE>::s_type(),
	 scalar_scalar_bo_OPNAME);
    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->add_operation
	(IA_ConstI<IA_Point<int>, LTYPE>::s_type(),
	 0,
	 scalar_iter_bo_OPNAME);

    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->add_operation
	(0,
	 IA_VectorI<IA_Point<int>, RTYPE>::s_type(),
	 iter_vec_bo_OPNAME);
    IA_Image_IP_ATYPE_LTYPE_RTYPE_OPNAME_tbl->add_operation
	(0,
	 IA_ConstI<IA_Point<int>, RTYPE>::s_type(),
	 iter_scalar_bo_OPNAME);

}
#endif // BOT_ATYPE_LTYPE_RTYPE_OPNAME
