~liggghts-dev/+junk/liggghts-dev

« back to all changes in this revision

Viewing changes to src/region.cpp

Committer: Anton Gladky
Date: 2011-12-03 07:21:52 UTC
Revision ID: gladky.anton@gmail.com-20111203072152-kuse4ike0ac0g8fx

New upstream version

files modified:
doc/Section_commands.html

doc/fix_insert.html

doc/fix_insert.txt

doc/pair_gran.html

src/compute_pair_gran_local.cpp

src/create_atoms.cpp

src/dump_stl.cpp

src/fix_insert.cpp

src/fix_insert.h

src/fix_insert_pack.cpp

src/fix_insert_pack.h

src/fix_insert_stream.cpp

src/fix_insert_stream.h

src/fix_particledistribution_discrete.cpp

src/fix_particledistribution_discrete.h

src/fix_pour.cpp

src/fix_rigid.h

src/fix_template_sphere.cpp

src/fix_template_sphere.h

src/fix_wall_gran_hooke_history.cpp

src/modify.cpp

src/modify.h

src/modify_liggghts.cpp

src/pair_gran.h

src/pair_gran_hooke_history.cpp

src/particleToInsert.cpp

src/particleToInsert.h

src/region.cpp

src/region_tetmesh.cpp

src/style_fix.h

src/version.h

Show diffs side-by-side

added added

removed removed

src/region.cpp

427

// impossible to calculate volume if bbox non-existent

428

if(!bboxflag) return 0;

429

430

double pos[3];

430

double pos[3],bbox_vol,region_vol;

431

int n_in = 0;

432

433

for(int i = 0; i < n_test; i++)

442

MyMPI::My_MPI_Sum_Scalar(n_in,world);

443

if(n_in == 0) error->all("Unable to calculate region volume - are you operating on a 2d region?");

444

445

double bbox_vol = (extent_xhi - extent_xlo) * (extent_yhi - extent_ylo) * (extent_zhi - extent_zlo);

446

return (static_cast<double>(n_in)/static_cast<double>(n_test) * bbox_vol);

445

bbox_vol = (extent_xhi - extent_xlo) * (extent_yhi - extent_ylo) * (extent_zhi - extent_zlo);

446

region_vol = static_cast<double>(n_in)/static_cast<double>(n_test) * bbox_vol;

447

return region_vol;

447

448

}

Older »