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>