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);
      }

}