$ry = sin(2 * M_PI * $l / $NL);
$cx = $rx * $RC;
$cy = $ry * $RC;
for ($c = 0; $c < $NC; ++$c) {
$x = $cx + $rx * cos(2 * M_PI * $c / $NC) * $RT;
$y = $cy + $ry * cos(2 * M_PI * $c / $NC) * $RT;
$z = sin(2 * M_PI * $c / $NC) * $RT;
print_v($x, $y, $z);
}
}
for ($l = 0; $l < $NL; ++$l) {
$rx = cos(2 * M_PI * $l / $NL);
$ry = sin(2 * M_PI * $l / $NL);
for ($c = 0; $c < $NC; ++$c) {
$x = $rx * cos(2 * M_PI * $c / $NC);
$y = $ry * cos(2 * M_PI * $c / $NC);
$z = sin(2 * M_PI * $c / $NC);
print_vn($x, $y, $z);
}
}
for ($l = 0; $l < $NL; ++$l) {
$l2 = ($l + 1) % $NL;
for ($c = 0; $c < $NC; ++$c) {
$c2 = ($c + 1) % $NC;
$v1 = 1 + $l * $NC + $c;
$v2 = 1 + $l2 * $NC + $c;
$v3 = 1 + $l2 * $NC + $c2;
$v4 = 1 + $l * $NC + $c2;
echo "f " . $v1 . "//" . $v1 . " " . $v2 . "//" . $v2 . " " . $v3 . "//" . $v3 . " " . $v4 . "//" . $v4 . "\n";
}
}
$x = cos(2 * M_PI * $c / $NC) * ($r * $R) / $NR;
$y = sin(2 * M_PI * $c / $NC) * ($r * $R) / $NR;
print_v($x, $y, $z);
}
}
$x = 0.0;
$y = 0.0;
print_v($x, $y, $z);
print_vn(0, 0, -1);
for ($c = 0; $c < $NC; ++$c) {
$x = cos(2 * M_PI * $c / $NC);
$y = sin(2 * M_PI * $c / $NC);
$z = 0;
print_vn($x, $y, $z);
}
print_vn(0, 0, 1);
$n = 1;
for ($r = 1; $r <= $NR; ++$r) {
$v0 = 2 + ($r - 1) * $NC;
for ($c = 0; $c < $NC; ++$c) {
$c2 = ($c + 1) % $NC;
if ($r == 1) {
echo "f 1//" . $n . " " . ($v0 + $c2) . "//" . $n . " " . ($v0 + $c) . "//" . $n . "\n";
} else {
echo "f " . ($v0 + $c) . "//" . $n . " " . ($v0 - $NC + $c) . "//" . $n . " " . ($v0 - $NC + $c2) . "//" . $n . " " . ($v0 + $c2) . "//" . $n . "\n";
}
}
}
for ($l = 0; $l < $NL; ++$l) {
$v0 = 2 + ($NR - 1) * $NC + $l * $NC;
for ($c = 0; $c < $NC; ++$c) {
