再追記修正 調べたら http://www.opengl-tutorial.org/jp/beginners-tutorials/tutorial-3-matrices/ こういうことだった。
float型だと手計算と一致せず、double型だと手計算と一致したのでdouble型にしました。
でも、三回目の計算がまったくあたらないので注意です。
後、計算が遅いです。
でも、平行移動のベクトルは
1 0 0 Tx x x+Tx*0 x
0 1 0 Ty y y+Ty*0 y
0 0 1 Tz * z ＝ z+Tz*0 = z
0 0 0 1 0 0 0 ベクトルのままではないのか （引用 ゲームを動かす技術と発想R ）
数値は、なんでも言いわけではなさそうだった。で、現状方程式で解いてる。（3Dグラフィックスのための数学では普通に４＊４の計算をしている。）
double a[16]={
1.000000,0.000000,0.000000,0.000000,
0.000000,1.000000,0.000000,0.000000,
0.000000,0.000000,1.000000,0.000000,
0.000000,0.000000,0.000000,1.000000};
double ｂ[16]={
1.000000,0.000000,0.000000,0.000000,
0.000000,1.000000,0.000000,0.000000,
0.000000,0.000000,1.000000,0.000000,
0.000000,0.000000,0.000000,1.000000};で数値は何でもいいとして
a×bの
4×4の行列(間違えたら嫌なので強調)を算出する簡単な方法がありましたら、教えてください。
-
気になる質問をクリップする
クリップした質問は、後からいつでもマイページで確認できます。
またクリップした質問に回答があった際、通知やメールを受け取ることができます。
クリップを取り消します
-
良い質問の評価を上げる
以下のような質問は評価を上げましょう
- 質問内容が明確
- 自分も答えを知りたい
- 質問者以外のユーザにも役立つ
評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。
質問の評価を上げたことを取り消します
-
評価を下げられる数の上限に達しました
評価を下げることができません
- 1日5回まで評価を下げられます
- 1日に1ユーザに対して2回まで評価を下げられます
質問の評価を下げる
teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。
- プログラミングに関係のない質問
- やってほしいことだけを記載した丸投げの質問
- 問題・課題が含まれていない質問
- 意図的に内容が抹消された質問
- 過去に投稿した質問と同じ内容の質問
- 広告と受け取られるような投稿
評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。
質問の評価を下げたことを取り消します
この機能は開放されていません
評価を下げる条件を満たしてません
質問の評価を下げる機能の利用条件
この機能を利用するためには、以下の事項を行う必要があります。
- 質問回答など一定の行動
-
メールアドレスの認証
メールアドレスの認証
-
質問評価に関するヘルプページの閲覧
質問評価に関するヘルプページの閲覧
+2
Frame_Scene_Root[0]とFrame_body[0]がfloat型にできなかったんだが、元の値が変わってしまうのですがどうすべきでしょうか。
なに言ってるかわからん。
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
+1
評価する人は自前に ４＊４のコードを作成し書き込んでいけっていったんだから、書き込んでいけよ。
<?php
$intA = [
3, 4, 5, 5,
6, 3, 4, 3,
3, 4, 5, 7,
9, 2, 3, 4
];
$intB = [
3, 4, 5, 5,
6, 3, 4, 3,
3, 4, 5, 7,
9, 2, 3, 4
];
$result = [];
for( $i = 0; $i < 16; $i++ ){
$result[] = $intA[$i] * $intB[$i];
}
for( $i = 0; $i < 4; $i++ ){
for( $j = 0; $j < 4; $j++ ){
echo $result[ ($i * 4) + $j ], ' ';
}
echo '<br />';
}もうすぐ定期試験2科目始まるのでレスは遅れますが悪しからず。
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
+1
a×bの 4×4の行列(間違えたら嫌なので強調)を算出する簡単な方法がありましたら
何らかのライブラリを使う，という方法．（マジレス）
一体何のために行列計算しているのかわからないけれども，その行列を使う目的によっては何らかのライブラリを導入している可能性があり，そうであって且つそこに行列を扱う機能が提供されているならばそれを使うのが手っ取り早い．
OpenCVとかEigenなら，行列用の型にoperator *が提供されているから，
A * Bとか書けばそれで済む．
OpenGLにはglMultMatrixなんてのがあるみたい（掛ける方向とか行列データの並び順とかの面が若干面倒そうかもだが）．
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
checkベストアンサー
0
簡単もなにもベクトルの内積を4x4回繰り返す以外なんかあるん?
[追記] 愚直な実装でよければ:
#include <stdio.h>
#define at(t,y,x) t[y*4+x]
int main() {
float a[16] = {
1.000000,0.000000,0.000000,0.000000,
0.000000,1.000000,0.000000,0.000000,
0.000000,0.000000,1.000000,0.000000,
0.000000,0.000000,0.000000,1.000000 };
float b[16] = {
1.000000,0.000000,0.000000,0.000000,
0.000000,1.000000,0.000000,0.000000,
0.000000,0.000000,1.000000,0.000000,
0.000000,0.000000,0.000000,1.000000 };
float c[16];
int row, col;
for ( row = 0; row < 4; ++row ) {
for ( col =0; col < 4; ++col ) {
float sum = 0.0f;
int t;
for ( t = 0; t <4; ++t ) {
sum += at(a,row,t) * at(b,t,col);
}
at(c,row,col) = sum;
}
}
// できたかな?
for ( row = 0; row < 4; ++row) {
for ( col = 0; col < 4; ++col) {
printf("%f ", at(c,row,col));
}
puts("");
}
return 0;
}投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
0
[別解] 関数 multiplies に仕立ててみた
#include <stdio.h>
#define at(t,y,x) t[y*4+x]
void multiplies(float* c, const float* a, const float* b) {
int row, col;
for ( row = 0; row < 4; ++row ) {
for ( col =0; col < 4; ++col ) {
float sum = 0.0f;
int t;
for ( t = 0; t <4; ++t ) {
sum += at(a,row,t) * at(b,t,col);
}
at(c,row,col) = sum;
}
}
}
int main() {
float a[16] = {
1.000000,0.000000,0.000000,0.000000,
0.000000,1.000000,0.000000,0.000000,
0.000000,0.000000,1.000000,0.000000,
0.000000,0.000000,0.000000,1.000000 };
float b[16] = {
1.000000,0.000000,0.000000,0.000000,
0.000000,1.000000,0.000000,0.000000,
0.000000,0.000000,1.000000,0.000000,
0.000000,0.000000,0.000000,1.000000 };
float c[16];
int row, col;
multiplies(c, a, b);
// できたかな?
for ( row = 0; row < 4; ++row) {
for ( col = 0; col < 4; ++col) {
printf("%f ", at(c,row,col));
}
puts("");
}
return 0;
}投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
0
解決したんじゃなかったんか?
a[4][4] と b[4][4] の積を c[4][4] に求める:
#include <stdio.h>
void multiplies(float c[4][4], const float a[4][4], const float b[4][4]) {
int row, col;
for ( row = 0; row < 4; ++row ) {
for ( col =0; col < 4; ++col ) {
float sum = 0.0f;
int t;
for ( t = 0; t <4; ++t ) {
sum += a[row][t] * b[t][col];
}
c[row][col] = sum;
}
}
}
int main() {
float a[4][4] = {
{ 1.000000,0.000000,0.000000,0.000000 },
{ 0.000000,1.000000,0.000000,0.000000 },
{ 0.000000,0.000000,1.000000,0.000000 },
{ 0.000000,0.000000,0.000000,1.000000 },
};
float b[4][4] = {
{ 1.000000,0.000000,0.000000,0.000000 },
{ 0.000000,1.000000,0.000000,0.000000 },
{ 0.000000,0.000000,1.000000,0.000000 },
{ 0.000000,0.000000,0.000000,1.000000 },
};
float c[4][4];
int row, col;
multiplies(c, a, b);
// できたかな?
for ( row = 0; row < 4; ++row) {
for ( col = 0; col < 4; ++col) {
printf("%f ", c[row][col]);
}
puts("");
}
return 0;
}投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
-5
引数で使った &Frame_Scene_Root[0],&Frame_body[0]が、逆だったかもしれません。逆の可能性が高いという感じです。
あと関数を何回も計算するのに、引数にconst を使うのは再利用の関点から、有り得ないことに気がつきました。自分が言いだしっぺです。すいません。
問題はここか。右から計算？
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
-6
手計算したんだけど、double型にしないと精度が駄目でした。
double Frame_Scene_Root[16]={
1.000000,0.000000,0.000000,0.000000,
0.000000,1.000000,0.000000,0.000000,
0.000000,0.000000,1.000000,0.000000,
0.000000,0.000000,0.000000,1.000000};
double Frame_body[16]={
1.278853,0.000000,-0.000000,0.000000,
0.000000,0.000000,1.123165,0.000000,
0.000000,-1.470235,0.000000,0.000000,
0.135977,2.027985,133.967667,1.000000
};
double Frame[16]={
1.000000,-0.000000,-0.000000,0.000000,
-0.000000,1.000000,0.000000,0.000000,
-0.000000,0.000000,1.000000,0.000000,
-0.142114,0.000023,-49.556850,1.000000};
double temp[16];
void Ranks( double *temp ,double* a, double *b)
{
temp[0]=a[0]*b[0]+a[1]*b[4]+a[2]*b[8]+a[3]*b[12];//ok
temp[1]=a[0]*b[1]+a[1]*b[5]+a[2]*b[9]+a[3]*b[13];//ok
temp[2]=a[0]*b[2]+a[1]*b[6]+a[2]*b[10]+a[3]*b[14];//ok
temp[3]=a[0]*b[3]+a[1]*b[7]+a[2]*b[11]+a[3]*b[15];//ok
temp[4]=a[4]*b[0]+a[5]*b[4]+a[6]*b[8]+a[7]*b[12];//ok
temp[5]=a[4]*b[1]+a[5]*b[5]+a[6]*b[9]+a[7]*b[13];//ok
temp[6]=a[4]*b[2]+a[5]*b[6]+a[6]*b[10]+a[7]*b[14];//ok
temp[7]=a[4]*b[3]+a[5]*b[7]+a[6]*b[11]+a[7]*b[15];//ok
temp[8]=a[8]*b[0]+a[9]*b[4]+a[10]*b[8]+a[11]*b[12];//0k
temp[9]=a[8]*b[1]+a[9]*b[5]+a[10]*b[9]+a[11]*b[13];//0k
temp[10]=a[8]*b[2]+a[9]*b[6]+a[10]*b[10]+a[11]*b[14];//0k
temp[11]=a[8]*b[3]+a[9]*b[7]+a[10]*b[11]+a[11]*b[15];//0k
temp[12]=a[12]*b[0]+a[13]*b[4]+a[14]*b[8]+a[15]*b[12];//0k
temp[13]=a[12]*b[1]+a[13]*b[5]+a[14]*b[9]+a[15]*b[13];
temp[14]=a[12]*b[2]+a[13]*b[6]+a[14]*b[10]+a[15]*b[14];
temp[15]=a[12]*b[3]+a[13]*b[7]+a[14]*b[11]+a[15]*b[15];
}
Ranks(&temp[0],&Frame_Scene_Root[0],&Frame_body[0]);ロジックがあってるかは別として書いてみたんだけどさ、
Frame_Scene_Root[0]とFrame_body[0]がdouble型にできなかったんだが、元の値が変わってしまうのですがどうすべきでしょうか。どこかを、constすればいいのか教えてください
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
-11
const char * char * const char * const
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
15分調べてもわからないことは、teratailで質問しよう！
- ただいまの回答率 88.92%
- 質問をまとめることで、思考を整理して素早く解決
- テンプレート機能で、簡単に質問をまとめられる
質問への追記・修正、ベストアンサー選択の依頼
Zuishin
2020/08/03 08:01
> 作らないやつは、他人をあれこれ評価できる権利ないっていってんだよ。
「作らない」はこの文脈では「回答しない」と読めますが、誰がそんなことを決めましたか？
質問者にはそれを決める権利はありません。
決める権利は運営が持っています。
現状、回答しなくても評価できるシステムになっていますし、それは合理的です。評価によって質問・回答がブラッシュアップされ、資産価値が増すからです。
むしろ回答しなければ評価できないシステムの方が不合理です。
episteme
2020/08/03 08:57
質問者は最も「作らないやつ」であるが。
2020/08/03 17:40
複数のユーザーから「問題・課題が含まれていない質問」という意見がありました
teratailでは、漠然とした興味から票を募るような質問や、意見の主張をすることを目的とした投稿は推奨していません。
「編集」ボタンから編集を行い、質問の意図や解決したい課題を明確に記述していただくと回答が得られやすくなります。
Zuishin
2020/08/09 10:54
解決済みなのになぜまだやってるのかわかりませんが、どうせ編集するなら評価するひと云々というところもついでに削除したらいいと思います。
txty
2020/08/09 11:02
まだやってた。だってわけわかんなくて解けないし。
＞どうせ編集するなら評価するひと云々というところもついでに削除。
そんな気になりますか。
Zuishin
2020/08/09 11:05
そりゃ通知が来たらどこが変わったか読むし、その時にゴミがいつまでも残ってたら気になるでしょう。
解けないなら解決済みにしなければ良かったと思いますが、なぜあの回答をもらって解けないのかわかりません。そのまま答えじゃないですか。
Zuishin
2020/08/09 11:06
で、人に尋ねる気がないなら編集せず一人で解けばいいんですよ。尋ねる気があるなら余計な憎まれ口は叩かず素直に聞けばいいと思いますね。
txty
2020/08/09 11:22 編集
正論ですが、外した途端にDOS攻撃みたいになったら嫌なので、お守りみたいなものです。ここでは保守？みたいなものでいます。
Zuishin
2020/08/09 11:27
その文言があるから書かないというおとなしい人は DOS 攻撃をしません。
txty
2020/08/09 11:32
どうだかわからないから、とりあえずはずしません。すいません。
Zuishin
2020/08/09 11:36
はずさないことで、この質問を見た人からあなたへの対応が親切でなくなることはあると思いますが、それもどうだかわからないと言うなら好きにすればいいと思います。ただ、それなら回答は期待しないことです。一人でどうぞ。
あまりしつこく無意味な通知を送ってくるなら、暑くてイライラしがちな季節なので、送られた人の態度もあまり良くない方向へと変化するかもしれませんね。
txty
2020/08/09 11:41
はい、あとはひとりでって方向にします。もう当初、欲しかったものは手に入ったので静かにしてます。
2020/08/09 14:33
複数のユーザーから「やってほしいことだけを記載した丸投げの質問」という意見がありました
「質問を編集する」ボタンから編集を行い、調査したこと・試したことを記入していただくと、回答が得られやすくなります。