@@ -88,29 +88,62 @@ vis, anim = visualize(mech_nc, storage_nc,
88
88
name= :nc ,
89
89
animation= anim)
90
90
91
- # ## MuJoCo cone
92
- color_mj = magenta;
91
+ # ## MuJoCo pyramidal cone
92
+ color_mjlc = magenta;
93
93
94
- mech_mj = get_mechanism (:box ,
94
+ mech_mjlc = get_mechanism (:box ,
95
95
timestep= timestep,
96
96
gravity= gravity,
97
97
friction_coefficient= friction_coefficient,
98
98
contact_type= :linear ,
99
- mode= :box , color= color_mj );
99
+ mode= :box , color= color_mjlc );
100
100
101
101
# ## Load
102
- initialize! (mech_mj , :box ,
102
+ initialize! (mech_mjlc , :box ,
103
103
x= x0,
104
104
q= one (UnitQuaternion),
105
105
v= v0,
106
106
ω= ω0)
107
- file = jldopen (joinpath (@__DIR__ , " ../MuJoCo_benchmark /results/cone_compare .jld2" ))
108
- storage_mj = generate_storage (mech_mj , [get_maximal_state (mech_mj ), file[" ztraj" ]. .. ])
107
+ file = jldopen (joinpath (@__DIR__ , " ../mujoco_benchmark /results/cone_compare_pyramidal .jld2" ))
108
+ storage_mjlc = generate_storage (mech_mjlc , [get_maximal_state (mech_mjlc ), file[" ztraj" ]. .. ])
109
109
110
110
# ## Visualize
111
- vis, anim = visualize (mech_mj, storage_mj ,
111
+ vis, anim = visualize (mech_mjlc, storage_mjlc ,
112
112
vis= vis,
113
- name= :mj ,
113
+ name= :mjlc ,
114
114
animation= anim)
115
115
116
+ line_mat_mjlc = LineBasicMaterial (color= color_mjlc, linewidth= 25.0 )
117
+ points_mjlc = Vector {Point{3,Float64}} ()
118
+ for (i, xt) in enumerate (storage_mjlc. x[1 ])
119
+ k = xt
120
+ push! (points_mjlc, Point (k[1 ], k[2 ], k[3 ]))
121
+ end
122
+ setobject! (vis[:path_mjlc ], MeshCat. Line (points_mjlc, line_mat_mjlc))
123
+
124
+
125
+ # ## MuJoCo elliptic cone
126
+ color_mjnc = RGBA (0 ,0 ,0 );
127
+
128
+ mech_mjnc = get_mechanism (:box ,
129
+ timestep= timestep,
130
+ gravity= gravity,
131
+ friction_coefficient= friction_coefficient,
132
+ contact_type= :linear ,
133
+ mode= :box , color= color_mjnc);
134
+
135
+ # ## Load
136
+ initialize! (mech_mjnc, :box ,
137
+ x= x0,
138
+ q= one (UnitQuaternion),
139
+ v= v0,
140
+ ω= ω0)
141
+ file = jldopen (joinpath (@__DIR__ , " ../mujoco_benchmark/results/cone_compare_elliptic.jld2" ))
142
+ storage_mjnc = generate_storage (mech_mjnc, [get_maximal_state (mech_mjnc), file[" ztraj" ]. .. ])
143
+
144
+ # ## Visualize
145
+ vis, anim = visualize (mech_mjnc, storage_mjnc,
146
+ vis= vis,
147
+ name= :mjnc ,
148
+ animation= anim)
116
149
0 commit comments