Skip to content
Snippets Groups Projects
Select Git revision
  • e4c8f816819076455bd2c854fca272550b3e6e04
  • master default protected
  • dof-renumbering
  • gdemesy-master-patch-30528
  • eval-space-time
  • oscillating_multiharm
  • MH_movement
  • axisqu
  • write_vtu_and_ensight_formats
  • movingband
  • CP_1972_add_vtu_file_writing
  • mortar
  • fast_freq_sweep_Resolution
  • applyresolvent_again
  • marteaua-master-patch-54323
  • patch-1
  • binde-master-patch-08072
  • binde-master-patch-52461
  • BCGSL
  • resolvent
  • TreeElementsOf
  • getdp_3_5_0
  • getdp_3_4_0
  • getdp_3_3_0
  • getdp_3_2_0
  • getdp_3_1_0
  • getdp_3_0_4
  • getdp_3_0_3
  • getdp_3_0_2
  • getdp_3_0_1
  • getdp_3_0_0
  • onelab_mobile_2.1.0
  • getdp_2_11_3 protected
  • getdp_2_11_2 protected
  • getdp_2_11_1 protected
  • getdp_2_11_0 protected
  • getdp_2_10_0 protected
  • getdp_2_9_2 protected
  • getdp_2_9_1 protected
  • getdp_2_9_0 protected
  • getdp_2_8_0 protected
41 results

Pos_Formulation.cpp

Blame
  • t16.jl 1017 B
    # See the corresponding Python tutorial for detailed comments.
    
    import gmsh
    
    gmsh.initialize()
    
    gmsh.model.add("t16")
    
    gmsh.model.occ.addBox(0,0,0, 1,1,1, 1)
    gmsh.model.occ.addBox(0,0,0, 0.5,0.5,0.5, 2)
    gmsh.model.occ.cut([(3,1)], [(3,2)], 3)
    
    x = 0; y = 0.75; z = 0; r = 0.09
    
    holes = []
    for t in 1:5
        global x, z
        x += 0.166
        z += 0.166
        gmsh.model.occ.addSphere(x,y,z,r, 3 + t)
        t = (3, 3 + t)
        push!(holes, t)
    end
    
    ov = gmsh.model.occ.fragment([(3,3)], holes)
    gmsh.model.occ.synchronize()
    
    lcar1 = .1
    lcar2 = .0005
    lcar3 = .055
    
    ov = gmsh.model.getEntities(0);
    gmsh.model.mesh.setSize(ov, lcar1);
    
    ov = gmsh.model.getBoundary(holes, false, false, true);
    gmsh.model.mesh.setSize(ov, lcar3);
    
    eps = 1e-3
    ov = gmsh.model.getEntitiesInBoundingBox(0.5-eps, 0.5-eps, 0.5-eps,
                                             0.5+eps, 0.5+eps, 0.5+eps, 0)
    gmsh.model.mesh.setSize(ov, lcar2)
    
    gmsh.model.mesh.generate(3)
    
    gmsh.write("t16.msh")
    
    if !("-nopopup" in ARGS)
        gmsh.fltk.run()
    end
    
    gmsh.finalize()