Fonc_Num_Computed * OP_Naire_Not_Comp::compute (const Arg_Fonc_Num_Comp & arg) { Fonc_Num_Computed * fc = _f.compute(arg); Tjs_El_User.ElAssert ( _accept_dim_in(fc->idim_out()), EEM0 << " Unexpected dimension for operator " << _name << "\n" << "| dim = " << fc->idim_out() ); INT dout = _dim_out(fc->idim_out()); if (fc->integral()) { if( _real) return new OP_un_Ndim_Comp_RI(arg,fc,arg.flux(),dout,_fri); else return new OP_un_Ndim_Comp_II(arg,fc,arg.flux(),dout,_fii); } else { if ( _real) return new OP_un_Ndim_Comp_RR(arg,fc,arg.flux(),dout,_frr); else return new OP_un_Ndim_Comp_IR(arg,fc,arg.flux(),dout,_fir); } }