Je voudrais effectuer une PCA sur un ensemble de données composé d'environ 40 000 échantillons, chaque échantillon affichant environ 10 000 fonctionnalités. L'utilisation de la fonction Matlab princomp prend systématiquement plus d'une demi-heure, moment auquel je tue le processus. Je voudrais...
CVXOPT: http://abel.ee.ucla.edu/cvxopt/index.html OpenOpt: http://openopt.org/Welcome Quelle est la relation entre eux? Quels sont leurs avantages / inconvénients, respectivement? BTW, existe-t-il une autre bibliothèque d'optimisation convexe à usage général de haute qualité pour Python / C ++ à...
Existe-t-il une liste publiée des meilleures pratiques pour assurer la longévité du code, en vue de résultats scientifiques reproductibles? (par exemple open source, pratiques de documentation, sélection des dépendances, sélection d'une langue, machines virtuelles, etc.). Connaître des études (ou,...
La plupart des bibliothèques mathématiques ont un certain nombre de versions de fonctions logarithmes. La plupart du temps, nous supposons qu'ils sont parfaits, mais en fait, beaucoup d'entre eux n'offrent qu'un certain nombre de chiffres de précision. Pour certaines fonctions, il existe des...
Bonjour communauté scicomp, J'ai travaillé dans le domaine des algorithmes de graphes en utilisant des frameworks tels que NetworkX (Python), JUNG et YFiles (Java). J'entre maintenant dans le domaine de l'informatique parallèle et haute performance. Pour un nouveau projet, je recherche une...
Ainsi , le théorème de décomposition de Cholesky que que toute réel symétrique définie positive matrice a une décomposition de Cholesky où est une matrice triangulaire inférieure.MMMM=LL⊤M=LL⊤M= LL^\topLLL Compte tenu , nous savons déjà , il y a des algorithmes rapides pour calculer le facteur de...
Je résout pour une énorme matrice définie positive clairsemée utilisant la méthode du gradient conjugué (CG). Est-il possible de calculer le déterminant de utilisant les informations produites lors de la résolution?A AA x =
Le calcul parallèle est fréquemment modélisé à l'aide d'un taux de calcul local déterministe, d'une surcharge de latence et d'une bande passante réseau. En réalité, ceux-ci sont spatialement variables et non déterministes. Des études telles que Skinner et Kramer (2005) observent des distributions...
Voulez-vous améliorer ce post? Fournissez des réponses détaillées à cette question, y compris des citations et une explication de la raison pour laquelle votre réponse est correcte. Les réponses sans suffisamment de détails peuvent être modifiées ou supprimées. Je suis entièrement nouveau dans la...
Disons que j'ai une fonction que je souhaite intégrer sur un tétraèdre T ⊂ R 3 . Si f était arbitraire, la quadrature de Gauss serait une bonne solution, mais je sais que f est harmonique. Dans quelle mesure la quadrature de Gauss peut-elle être accélérée en utilisant ces informations?F: R3→...
Théoriquement, comment le temps pour faire une échelle de calcul de la théorie fonctionnelle de la densité (DFT) avec le nombre d'électrons? Je m'intéresse aux implémentations DFT "typiques" telles que VASP, ABINIT, etc., pas aux codes O
J'ai un système linéaire inhomogène A x = bUNEX=b Ax=b où UNEUNEA est une matrice n × nn×nn\times n réelle avec n ≤ 4n≤4n\leq 4 . L'espace nul de UNEUNEA est garanti d'être de dimension nulle, l'équation a donc un inverse unique x = A- 1bX=UNE-1bx=A^{-1} b . Puisque le résultat entre dans le côté...
Étant donné une fonction inconnue , nous pouvons évaluer sa valeur en tout point de son domaine, mais nous n'avons pas son expression. En d'autres termes, f est comme une boîte noire pour nous.f:Rd→Rf:Rd→Rf:\mathbb R^d \to \mathbb Rfff Quel est le nom du problème de trouver le minimiseur de ?...
Comment MATLAB, par exemple, calcule-t-il la SVD d'une matrice donnée? Je suppose que la réponse implique probablement de calculer les vecteurs propres et les valeurs propres de A*A'. Si tel est le cas, j'aimerais également savoir comment calcule-t-on ces
J'ai un problème physique régi par l'équation de Poisson en deux dimensions J'ai des mesures des deux composantes du gradient ∂ u / ∂ x et ∂ u / ∂ y le long d'une partie de la frontière, Γ m , je voudrais donc imposer ∂ u- ∇2u = f( x , y) ,je nΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y), \; in \; \Omega...
Au travail, je travaille essentiellement comme consultant indépendant. Pour la direction et les clients, j'ai besoin d'estimer le temps qu'il faudra pour développer un logiciel dans le cadre de ma recherche en science informatique. Cependant, mes estimations de temps sont généralement désactivées....
J'ai un ensemble de données de points en 3 dimensions pour lesquels je voudrais construire un maillage, en utilisant python. Tous les logiciels que j'ai vus exigent que vous fournissiez les bords. Existe-t-il un programme en python qui prend en entrée un ensemble de points en 3D et sort un...
En commençant par l'équation d'advection sous forme de conservation. ut= ( a ( x ) u )Xut=(une(X)u)X u_t = (a(x)u)_x où a ( x )une(X)a(x) est une vitesse qui dépend de l'espace, et uuu est une concentration d'une espèce qui est conservée. La discrétisation du flux (où le flux F= a ( x )...
Je développe actuellement une méthode de décomposition de domaine pour la solution du problème de diffusion. Fondamentalement, je résous un système de BVP Helmholtz de manière itérative. Je discrétise les équations en utilisant la méthode des éléments finis sur des maillages triangulaires ou...
On me donne une matrice symétrique, inversible, définie positive et dense. Je dois tester si où J est la matrice des uns.12×1212×1212 \times 12QQQJdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Je fais actuellement cela avec la bibliothèque de tatou mais cela...