ich kann diese vektoren addieren und multiplizieren mit skalaren, aber ich möchte sowas in der art machen:
v1=( x1, y1, z1 ), v2=( x2, y2, z2 ), v3=( x1*x2, y1*y2, z1*z2 )

und das geht nicht, weil ich nicht an die komponenten ran komm. aber wie gesagt geht das anscheinend. bloß wie ist die frage...

die anwendung von dem multiplies<T>() ist mir leider auch noch nicht ganz klar, zumindest was die verwendung außerhalb von partial_sum, accumulate & co angeht...

So?

Für die Fälle, in denen man das doch braucht (z.B. Skalierung), schreib dir eine Funktion ComponentwiseProduct() oder ähnlich. So hab ich es auch gemacht.

Pout<<sf_plus[1].x<<endl; (normalerweise mit .x(), .y(), .z() zugriff)
Fehler: Abfrage des Elementes »x« in »((Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*)sf_plus)->Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>::<anonymous>.Foam::DimensionedField<double, Foam::surfaceMesh>::<anonymous>.Foam::Field<double>::<anonymous>.Foam::List<double>::<anonymous>.Foam::UList<T>::operator[] [with T = double](1)«, das vom Nicht-Klassentyp »double« ist

Pout<<sf_plus[1].x()<<endl;
Fehler: Abfrage des Elementes »x« in »((Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*)sf_plus)->Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>::<anonymous>.Foam::DimensionedField<double, Foam::surfaceMesh>::<anonymous>.Foam::Field<double>::<anonymous>.Foam::List<double>::<anonymous>.Foam::UList<T>::operator[] [with T = double](1)«, das vom Nicht-Klassentyp »double« ist

also funktioniert halt alles nicht.
ich hab das ganz auch mit "multiplies" probiert:
std::multiplies <Type> mul =std::multiplies <Type> ();
Pout<<mul(sf_plus[1],sf_plus[1]);
aber das brachte als fehler:
Fehler: Umwandlung von »Foam::Tensor<double>« in nicht-skalaren Typen »Foam::Vector<double>« angefordert

mir wurde gesagt, dass es einen befehl so ähnlich wie multiply_coponents oder irgendwie abgekürzt mtply_cmp oder so geben würde, der trotzdem die multiplikation machen würde. aber außer dem multiplies hab ich nichts gefunden und der macht das leider nicht.
Kennt von euch jemand den befehl oder wie man den multiplies geeigneter einsetzt oder ist das problem echt nur mit grundlegenden änderungen der vector/tensor klasse möglich?

Klar ist, dass Du Zugriff auf die einzelnen Member brauchst, das sollte auch so ziemlich jede Vector-Klasse anbieten. Basierend darauf kannst Du eine freie Funktion implementieren, welche die Multiplikation für Dich erfüllt. Ob die Bibliothek sonst was mit multiply anbietet, weiß ich nicht, da würde ich einfach kurz die Doku durchsuchen und wenn sich nix findet halt selbst schreiben.

der rest war eher zur verdeutlichung, dass ich hier nicht eine neue vektorklasse machen will, sondern mit bestehenden zu kämpfen hab

Also nochmal: Freie Funktion mit Zugriff auf die Elemente des Vektors und selbst schreiben, falls diese FOAM-Bibliothek nichts anbietet. Sollte doch nicht so schwierig sein.

zugegebenerweiße recht versteckt. google findet für den begriff auch nur 2 treffer.

ich werd mich mal an der funktion versuchen.

das mit dem zugriff ist aber irgendwie merkwürdig. ich schreibe was für openfoam (ist n programm zur simulation von thermodynamik zeugs). An sich versteh ich nicht ganz warum der Zugriff auf die Elemente später nicht mehr geht...
wisst ihr was der grund ist? warum geht der zugriff nach dem ersten mal und danach nicht mehr?

am anfang hab ich eine Liste mit vektoren und ich kann auf jede komponente zugreifen:
Pout<<"icoFoam.C u wert: " << U.internalField()[0][0]<< endl;

dies übergeb ich dann:
fvm::div(phi, U)

template<class Type> tmp<fvMatrix<Type> >
div( const tmp<surfaceScalarField>& tflux, GeometricField<Type, fvPatchField, volMesh>& vf
)
vf ist hier die Liste mit den U-Vektoren...
und kann dann ebenfalls noch drauf zugreifen:
Pout<<"fvmDiv.C div3 u wert: " << vf.internalField()[0][0] << endl;

aber bei der nächsten übergabe gehts nicht mehr...
fv::convectionScheme<Type>::New (vf.mesh(), flux, vf.mesh().divScheme(name) )().fvmDiv(flux, vf);

template<class Type>
tmp<fvMatrix<Type> >
gaussConvectionScheme<Type>::fvmDiv( const surfaceScalarField& faceFlux, GeometricField<Type, fvPatchField, volMesh>& vf) const

und in fvmDiv geht nix mehr mit komponentenzugriff...
http://www.c-plusplus.net/forum/images/smiles/confused.gif

Am Rande: Gewöhn Dir mal an, das herzallerliebste Code-Tag ([cpp][/cpp]) zu benutzen.

Dann sehe ich in deinem Code keinen Vector mehr, sprachen wir nicht über den? Und dann solltest Du vll noch "geht nicht" genauer definieren. Ich könnte jetzt raten, aber das muss ja nicht sein.

Ansonsten hat Schaum laut Doku sehr wohl Elementzugriffe über x(), y(), z(): http://foam.sourceforge.net/doc/Doxygen/html/classFoam_1_1Vector.html

Also seh ich nicht, wieso das hier für nen Vektor nicht klappt.

template<typename T>
Foam::Vector<T> multiplyComponentWise(const Foam::Vector<T>& v1, const Foam::Vector<T>& v2)
{
    Foam::Vector<T> result(v1);
    result.x() *= v2.x();
    result.y() *= v2.y();
    result.z() *= v2.z();
    return result;
}

Aber ich glaube, ich habe dein Problem immer noch nicht so wirklich verstanden.

wenn das gehen würde, wär die multiplikation ja kein problem...