422:         for (int iel = 0; iel < nels; ++iel)

423:         {

424:           // Operate on electron with index iel

425:           els.setActive(iel);

426:           // Compute gradient at the current position

427:           Timers[Timer_evalGrad].get().start();

428:           PosType grad_now = wavefunction.evalGrad(els, iel);

429:           Timers[Timer_evalGrad].get().stop();

430: 

431:           // Construct trial move

432:           els.makeMove(iel, delta[iel]);

433: 

434:           // Compute gradient at the trial position

435:           Timers[Timer_ratioGrad].get().start();

436:           PosType grad_new;

437:           wavefunction.ratioGrad(els, iel, grad_new);

438:           Timers[Timer_ratioGrad].get().stop();

439: 

440:           // Accept/reject the trial move

441:           if (ur[iel] < accept) // MC

442:           {

443:             // Update position, and update temporary storage

444:             Timers[Timer_Update].get().start();

445:             wavefunction.acceptMove(els, iel);

446:             Timers[Timer_Update].get().stop();

447:             els.acceptMove(iel);

448:             my_accepted++;

449:           }

450:           else

451:           {

452:             els.rejectMove(iel);

453:             wavefunction.restore(iel);