A fix in energy.c. Normals are summed now.

Miha Fošnarič

2016-07-11 02d65c3c9baa0e76744ce001c55ec699ab3b956f

 src/main.c

@@ -30,6 +30,9 @@
   ts_tape *tape;
   ts_uint start_iteration=0;
   force_from_tape=0;
   /* Area and volume for constant area and constant volume are initialized to be zero */
   A0=0;
   V0=0;
   /* create lock file */
   createPidFile("ts_trisurf",".lock",0);
   parse_args(argv,argc); // sets global variable command_line_args (defined in io.h)
@@ -55,7 +58,6 @@
            }
         fclose(fd);
         start_iteration++;
         printf("nucleus coords: %.17e %.17e %.17e\n",vesicle->nucleus_center[0], vesicle->nucleus_center[1], vesicle->nucleus_center[2]);
      }
      else
         ts_fprintf(stdout,"No .status file. The iteration count will start from 0");
@@ -123,7 +125,9 @@
      }
      }
   }
         //printf("nucleus coords: %.17e %.17e %.17e\n",vesicle->nucleus_center[0], vesicle->nucleus_center[1], vesicle->nucleus_center[2]);

         //write_vertex_xml_file(vesicle,1000);
   run_simulation(vesicle, tape->mcsweeps, tape->inititer, tape->iterations, start_iteration);
   write_master_xml_file(command_line_args.output_fullfilename);
   write_dout_fcompat_file(vesicle,"dout");