Thanks for rapid reply. <br>If decreasing the size of IO buffer, will it affect the speed?<br>With respect to the second way, jobs normally don't stop automatically because I am doing MD except the walltime is reached. So, cp2k may not have time to release the buffer.<br>
<br>Jun<br><br><div class="gmail_quote">On Fri, Dec 5, 2008 at 10:44 AM, Teodoro Laino <span dir="ltr"><<a href="mailto:teodor...@gmail.com">teodor...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Yes, I can give at least two possibilities:<br>
1) change the default of the machine on which you are running (you can<br>
increase/decrease the size of the IO buffer, on almost<br>
all machine starting from commodore64. I would suggest to get in touch<br>
with your system administrator for that.)<br>
2) wait that the job is finished. At that point the buffer is released<br>
and you will have your file with all the bytes at the proper place.<br>
<br>
Teo<br>
<div class="Ih2E3d"><br>
Jun Cheng wrote:<br>
> Hi Teo,<br>
><br>
</div><div class="Ih2E3d">> Thanks a lot for the help. It works well except the third individual<br>
> output file is still empty even with applying FLUSH_SHOULD_FLUSH TRUE.<br>
> It is probably because the third force_eva section is just a classic<br>
> restraint to a bond to prevent it from breaking during MD runs. So its<br>
> output may be too small to reach the size of the IO buffer in a<br>
> noticeable time. In fact, I did similar runs before and never got this<br>
> output in 12 hours (the walltime limit). Is there a way to overcome<br>
> this problem?<br>
><br>
> Cheers,<br>
> Jun<br>
><br>
><br>
> On Thu, Dec 4, 2008 at 8:51 PM, Laino Teodoro <<a href="mailto:teodor...@gmail.com">teodor...@gmail.com</a><br>
</div><div><div></div><div class="Wj3C7c">> <mailto:<a href="mailto:teodor...@gmail.com">teodor...@gmail.com</a>>> wrote:<br>
><br>
> Dear Jun,<br>
><br>
> 1) you don't see nothing in the individual .out files because the<br>
> IO is buffered.<br>
> Depending on the size of the IO buffer on the machine that you are<br>
> using it may take a WHILE<br>
> to see something in the output file.<br>
> If you want to flush you can play with this keyword:<br>
> <a href="http://cp2k.berlios.de/manual/CP2K_INPUT/GLOBAL.html#desc_FLUSH_SHOULD_FLUSH" target="_blank">http://cp2k.berlios.de/manual/CP2K_INPUT/GLOBAL.html#desc_FLUSH_SHOULD_FLUSH</a><br>
><br>
> 2) Switching to medium you are activating some print_key, in this<br>
> case DIPOLE, which for the<br>
> specific force_eval you are using (MIXED) requires that DIPOLE has<br>
> to be defined in all<br>
> force_eval underneath. To skip this problem you can disable by<br>
> hand the print_key DIPOLE:<br>
> <a href="http://cp2k.berlios.de/manual/CP2K_INPUT/FORCE_EVAL/MIXED/PRINT/DIPOLE.html" target="_blank">http://cp2k.berlios.de/manual/CP2K_INPUT/FORCE_EVAL/MIXED/PRINT/DIPOLE.html</a><br>
><br>
> .....<br>
> &DIPOLE OFF<br>
> &END<br>
> .....<br>
><br>
> This will solve your issue.<br>
> Teo<br>
><br>
> On 4 Dec 2008, at 19:15, Jun Cheng wrote:<br>
><br>
>> Hi Teo,<br>
>><br>
>> I am very sorry about this.<br>
>> The file (test.tgz) is uploaded in group files.<br>
>> Thanks.<br>
>><br>
>> Jun<br>
>><br>
>> On Thu, Dec 4, 2008 at 6:00 PM, Teodoro Laino<br>
</div></div><div><div></div><div class="Wj3C7c">>> <<a href="mailto:teodor...@gmail.com">teodor...@gmail.com</a> <mailto:<a href="mailto:teodor...@gmail.com">teodor...@gmail.com</a>>> wrote:<br>
>><br>
>><br>
>> and where is your input file?<br>
>> What you describe has not too much sense..anyway..<br>
>> Provide your input file and all files necessary to<br>
>> investigate the<br>
>> problem (i.e. run it).<br>
>> This is a general rule when asking for help, and I'm<br>
>> surprised to see<br>
>> that one has always to remind that.<br>
>> Teo<br>
>><br>
>> Jun wrote:<br>
>> > Hi all,<br>
>> ><br>
>> > I updated my cp2k yesterday, and did some tests on mixed multi<br>
>> > force_eva calcuations. I've got some problems on the outputs:<br>
>> > (i) When I used print_level low, the job got running.<br>
>> Everything<br>
>> > looked fine except for the empty output files corresponding<br>
>> to the<br>
>> > individual force_eva.(the overall output file, *ener file<br>
>> and *xyz<br>
>> > were fine).<br>
>> > (ii) Then, I tried print_low medium. I got the individual<br>
>> output<br>
>> > files, but the job stopped after initialization. The error<br>
>> message<br>
>> > is:<br>
>> > =============================================================<br>
>> > THERMOSTAT| Thermostat Info for PARTICLES<br>
>> > THERMOSTAT| Type of thermostat<br>
>> Nose-<br>
>> > Hoover-Chains<br>
>> > THERMOSTAT| Nose-Hoover-Chain<br>
>> > length 3<br>
>> > THERMOSTAT| Nose-Hoover-Chain time constant<br>
>> > [ fs] 1000.00<br>
>> > THERMOSTAT| Order of Yoshida<br>
>> > integrator 3<br>
>> > THERMOSTAT| Number of multiple time<br>
>> > steps 2<br>
>> > THERMOSTAT| Initial Potential<br>
>> > Energy 0.000000<br>
>> > THERMOSTAT| Initial Kinetic<br>
>> > Energy 0.000523<br>
>> > THERMOSTAT| End of Thermostat Info for PARTICLES<br>
>> ><br>
>> > ********************** begin of velocity initialization<br>
>> > ***********************<br>
>> > Initial<br>
>> > Temperature<br>
>> 330.00<br>
>> > K<br>
>> > COM velocity: 0.000000000000 0.000000000000<br>
>> > 0.000000000000<br>
>> > *********************** end of velocity initialization<br>
>> > ************************<br>
>> ><br>
>> ><br>
>> > *** 17:07:15 WARNING in force_env_methods:mixed_energy_forces<br>
>> > err=-300 ***<br>
>> > *** ASSERTION (cond) failed at line 1203 Error<br>
>> (0.259964E-11)<br>
>> > in ***<br>
>> > *** computing numerical derivatives larger then<br>
>> > (0.100000E-11) . ***<br>
>> > *<br>
>> > *** 17:07:15 ERRORL2 in cp_result_methods:get_result_r1<br>
>> > processor 0 ***<br>
>> > *** err=-300 Trying to access result ([DIPOLE]) which<br>
>> was never<br>
>> > stored! ***<br>
>> > *** /home/e89/e89/jc590/work/cp2k/makefiles/../src/<br>
>> > cp_result_methods.F ***<br>
>> > *** line<br>
>> > 164<br>
>> ***<br>
>> > *<br>
>> ><br>
>> > ===== Routine Calling Stack =====<br>
>> ><br>
>> > 2 qs_mol_dyn_low<br>
>> > 1 CP2K<br>
>> > CP2K| Trying to access result ([DIPOLE]) which was never<br>
>> stored! /<br>
>> ><br>
>> home/e89/e89/jc590/work/cp2k/makefiles/../src/cp_result_methods.F<br>
>> > line 164<br>
>> > CP2K| Abnormal program termination, stopped by process<br>
>> number 0<br>
>> > CP2K| Trying to access result ([DIPOLE]) which was never<br>
>> stored! /<br>
>> ><br>
>> home/e89/e89/jc590/work/cp2k/makefiles/../src/cp_result_methods.F<br>
>> > line 164<br>
>> > CP2K| Trying to access result ([DIPOLE]) which was never<br>
>> stored! /<br>
>> ><br>
>> home/e89/e89/jc590/work/cp2k/makefiles/../src/cp_result_methods.F<br>
>> > line 164<br>
>> > CP2K| Abnormal program termination, stopped by process<br>
>> number 65<br>
>> > CP2K| Abnormal program termination, stopped by process<br>
>> number 64<br>
>> > CP2K| Trying to access result ([DIPOLE]) which was never<br>
>> stored! /<br>
>> ><br>
>> home/e89/e89/jc590/work/cp2k/makefiles/../src/cp_result_methods.F<br>
>> > line 164<br>
>> > CP2K| Abnormal program termination, stopped by process<br>
>> number 14<br>
>> ><br>
>> ===================================================================<br>
>> ><br>
>> > Does anyone know what is the problem? Any reply will be<br>
>> appreciated.<br>
>> ><br>
>> > Cheers,<br>
>> > Jun<br>
>> ><br>
>> > P.S. I will upload my inputs if necessary.<br>
>> ><br>
>> > ><br>
>> ><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
><br>
><br>
><br>
><br>
><br>
> ><br>
<br>
<br>
<br>
</div></div></blockquote></div><br>