Tutorials - plate_hole.obj not working

[sky-yunxiao]

Describe the bug
Hello, I am using Ubuntu20.04, and run unit_tests with no problems. However, when trying to run through the tutorial plate_hole.obj, I meet run_time error

[2023-08-02 14:00:49.264] [polyfem] [info] Saving output to
[2023-08-02 14:00:49.265] [polyfem] [info] Loading mesh ...
[2023-08-02 14:00:49.266] [geogram] [info] Loading file /home/guanyx/mechanics/polyfem/build/plate_hole.obj...
[2023-08-02 14:00:49.274] [geogram] [info] (FP64) nb_v:8549 nb_e:0 nb_f:16797 nb_b:299 tri:1 dim:3
[2023-08-02 14:00:49.274] [geogram] [info] Attributes on vertices: point[3]
[2023-08-02 14:00:49.287] [polyfem] [debug] -- bbox before normalization:
[2023-08-02 14:00:49.287] [polyfem] [debug] min : 0.500183 0.5
[2023-08-02 14:00:49.287] [polyfem] [debug] max : 100.5 50.5
[2023-08-02 14:00:49.287] [polyfem] [debug] extent: 99.999817 50
[2023-08-02 14:00:49.287] [polyfem] [debug] -- bbox after normalization:
[2023-08-02 14:00:49.287] [polyfem] [debug] min : 0 0
[2023-08-02 14:00:49.287] [polyfem] [debug] max : 1 0.5000009150016745
[2023-08-02 14:00:49.287] [polyfem] [debug] extent: 1 0.5000009150016745
[2023-08-02 14:00:49.287] [polyfem] [debug] -- edge length after normalization:
[2023-08-02 14:00:49.287] [polyfem] [debug] min: 0.004206997358882737
[2023-08-02 14:00:49.287] [polyfem] [debug] max: 0.018750034312562808
[2023-08-02 14:00:49.287] [polyfem] [debug] avg: 0.00831950579693002
[2023-08-02 14:00:49.289] [polyfem] [info] mesh bb min [[0, 0]], max [[1, 0.5000009150016745]]
[2023-08-02 14:00:49.290] [polyfem] [info] took 0.02504s
[2023-08-02 14:00:49.290] [polyfem] [info] Loading obstacles...
[2023-08-02 14:00:49.290] [polyfem] [info] took 0.000399s
[2023-08-02 14:00:49.290] [polyfem] [info] simplex_count: 16797
[2023-08-02 14:00:49.290] [polyfem] [info] regular_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] regular_boundary_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] simple_singular_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] multi_singular_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] boundary_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] multi_singular_boundary_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] non_regular_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] non_regular_boundary_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] undefined_count: 0
[2023-08-02 14:00:49.290] [polyfem] [info] total count: 16797
[2023-08-02 14:00:49.290] [polyfem] [info] Building isoparametric basis...
[2023-08-02 14:00:49.307] [polyfem] [info] Counting flipped elements...
[2023-08-02 14:00:49.315] [polyfem] [info] done
[2023-08-02 14:00:49.315] [polyfem] [info] Building collision mesh...
[2023-08-02 14:00:49.317] [polyfem] [info] Done!
[2023-08-02 14:00:49.317] [polyfem] [debug] Building node mapping...
[2023-08-02 14:00:49.321] [polyfem] [debug] Done (took 0.004075s)
[2023-08-02 14:00:49.322] [polyfem] [info] hmin: 0.004206997358882737
[2023-08-02 14:00:49.322] [polyfem] [info] hmax: 0.018750034312562808
[2023-08-02 14:00:49.322] [polyfem] [info] havg: 0.00831950579693002
[2023-08-02 14:00:49.322] [polyfem] [info] n_bases 8549
[2023-08-02 14:00:49.322] [polyfem] [info] took 0.016533s
[2023-08-02 14:00:49.322] [polyfem] [info] flipped elements 0
[2023-08-02 14:00:49.322] [polyfem] [info] h: 0.018750034312562808
[2023-08-02 14:00:49.322] [polyfem] [info] n bases: 8549
[2023-08-02 14:00:49.322] [polyfem] [info] n pressure bases: 0
[2023-08-02 14:00:49.322] [polyfem] [info] Building cache...
[2023-08-02 14:00:49.338] [polyfem] [info] took 0.016645s
[2023-08-02 14:00:49.338] [polyfem] [info] Assigning rhs...
[2023-08-02 14:00:49.389] [polyfem] [info] took 0.050380999999999995s
[2023-08-02 14:00:49.389] [polyfem] [info] Solving Laplacian
[2023-08-02 14:00:49.389] [polyfem] [info] Solve using Eigen::SparseLU linear solver
[2023-08-02 14:00:49.389] [polyfem] [info] Eigen::SparseLU...
[2023-08-02 14:00:49.391] [polyfem] [info] Assembling stiffness mat...
[2023-08-02 14:00:49.397] [polyfem] [info] took 0.005912s
[2023-08-02 14:00:49.397] [polyfem] [info] sparsity: 59239/73085401
[2023-08-02 14:00:49.434] [polyfem] [debug] Solver error: 4.367620351922937e-14
[2023-08-02 14:00:49.434] [polyfem] [info] took 0.045617s
[2023-08-02 14:00:49.434] [polyfem] [info] Computing errors...
[2023-08-02 14:00:49.448] [polyfem] [info] took 0.013875s
[2023-08-02 14:00:49.448] [polyfem] [info] -- L2 error: 0.45393639236472705
[2023-08-02 14:00:49.448] [polyfem] [info] -- Lp error: 0.8016236718410377
[2023-08-02 14:00:49.448] [polyfem] [info] -- H1 error: 1.441222086165727
[2023-08-02 14:00:49.448] [polyfem] [info] -- H1 semi error: 1.3678679955824635
[2023-08-02 14:00:49.448] [polyfem] [info] -- Linf error: 1.2195435452710095
[2023-08-02 14:00:49.448] [polyfem] [info] -- grad max error: 2.107988764518685
[2023-08-02 14:00:49.448] [polyfem] [info] total time: 0.068062s
[2023-08-02 14:00:49.448] [polyfem] [info] Saving json...
[2023-08-02 14:00:49.449] [polyfem] [info] done
[2023-08-02 14:00:49.459] [polyfem] [error] Solve the problem first!
[2023-08-02 14:00:49.459] [polyfem] [error] Assemble grad not implemented by Laplacian!
terminate called after throwing an instance of 'std::runtime_error'
what(): Assemble grad not implemented by Laplacian!

Platform

• Windows
• macOS
• Linux