From maechler at stat.math.ethz.ch  Sat Apr  5 16:16:05 2008
From: maechler at stat.math.ethz.ch (Martin Maechler)
Date: Sat, 5 Apr 2008 16:16:05 +0200
Subject: [Rd] optim: why is REPORT not used in SANN?
In-Reply-To: <47DD17AF.5030606@TU-Dresden.de>
References: <47DD17AF.5030606@TU-Dresden.de>
Message-ID: <18423.35237.59390.18811@cmath-5.math.ethz.ch>

>>>>> "TP" == Thomas Petzoldt <Thomas.Petzoldt at tu-dresden.de>
>>>>>     on Sun, 16 Mar 2008 13:50:55 +0100 writes:

    TP> Hello, I wonder why the control parameter REPORT is not
    TP> supported by method SANN. Looking into optim.c I found
    TP> an internal constant:

    TP> #define STEPS 100

    TP> ... and decreasing this to 10 helped me fine-tuning the
    TP> annealing parameters in an actual problem.

    TP> Is there any reason why not passing nREPORT to samin and
    TP> setting something like:

    TP> STEPS = nREPORT / tmax

Sorry to reply late (but then, rather than never ..).

You ask for reasons... I see/guess :

- the  SANN  method also was contributed from "outside" 
  (as ?optim mentions); and the original authors may not have
  seen a use for such more flexible monitoring.

- the R core members are probably not using 'samin' very often

- If there is a need you can write the function you are
  optimizing in a way that it prints info.

- Nobody has contributed a well-tested patch against R-devel to
  both code and documentation
  which would implement your proposal ___ BACK COMPATIBLY __
  (i.e. the default for SANN should remain to print every 100th;
   and this differs from the default for BFGS where the default
  'REPORT' leads to output every 10th eval).

Regards,
Martin

    TP> Thomas P.



    TP> -- Thomas Petzoldt Technische Universitaet Dresden
    TP> Institut fuer Hydrobiologie
    TP> thomas.petzoldt at tu-dresden.de 01062 Dresden
    TP> http://tu-dresden.de/hydrobiologie/ GERMANY



[Rd] optim: why is REPORT not used in SANN?
Martin Maechler