int main(int argc, const char * argv[]) { Electrostatics e; TriangleMesh *teststl = new TriangleMesh(); teststl->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere5120.stl", MeshFileFormat::MFF_STL); teststl->write("/Users/phaedon/github/bem-laplace-simple/meshes/sphere5120.obj", MeshFileFormat::MFF_OBJ); TriangleMesh *s = new TriangleMesh(); s->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere_5mm_h0cm_s3.obj", MeshFileFormat::MFF_OBJ); e.addBubble(s); s = new TriangleMesh(); s->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere_5mm_h10cm.obj", MeshFileFormat::MFF_OBJ); e.addBubble(s); s = new TriangleMesh(); s->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere_5mm_h30cm.obj", MeshFileFormat::MFF_OBJ); e.addBubble(s); s = new TriangleMesh(); s->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere_5mm_h50cm.obj", MeshFileFormat::MFF_OBJ); e.addBubble(s); s = new TriangleMesh(); s->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere_5mm_h80cm.obj", MeshFileFormat::MFF_OBJ); e.addBubble(s); s = new TriangleMesh(); s->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere_5mm_h90cm.obj", MeshFileFormat::MFF_OBJ); e.addBubble(s); s = new TriangleMesh(); s->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere_5mm_h95cm.obj", MeshFileFormat::MFF_OBJ); e.addBubble(s); s = new TriangleMesh(); s->read("/Users/phaedon/github/bem-laplace-simple/meshes/sphere_5mm_h995mm_s3.obj", MeshFileFormat::MFF_OBJ); e.addBubble(s); TriangleMesh *plane = new TriangleMesh(); plane->read("/Users/phaedon/github/bem-laplace-simple/meshes/plane_h1_s50mm_t128.obj", MeshFileFormat::MFF_OBJ); e.setSurface(plane); std::vector<double> caps; e.capacitance(caps); for (size_t i = 0; i < caps.size(); i++) { std::cout << "capacitance of " << i << "th bubble: " << caps[i] << std::endl; } std::cout << "Finished!" << std::endl; return 0; }