Skip to main content

Computing many quadratic forms in numpy

Suppose I have a collection of square matrices Sigma1, Sigma2, ... . These are stored in a numpy array Sigma such that Sigma[0, :, :] = Sigma1, Sigma[1, :, :] = Sigma2, etc.

I have a corresponding set of vectors v1, v2, ..., stored in a matrix v such that v[0, :] = v1, v[1, :] = v2, etc.

I am trying to produce a vector a such that a[0] is v1.T @ Sigma1 @ v1 (a scalar value), a[1] is v2.T @ Sigma2 @ v2, etc.

I'm a bit bewildered by the documentation for einsum tensordot and the like, which seem promising but which I have been unable to coerce into giving me the desired answer.

For example:

Sigma = np.empty(4, 3, 3)
for i in range(4):
    Z = np.random.randn(10, 3)
    Sigma[i, :, :] = Z.T @ Z

y = np.random.randn(4, 3)

I want a function foo(Sigma y) that returns a vector equivalent to

bar = np.empty(4)
for i in range 4:
    bar[i] = y[i, :].T @ Sigma[i, :, :] @ y[i, :]
return bar


source https://stackoverflow.com/questions/70385907/computing-many-quadratic-forms-in-numpy

Comments

Popular posts from this blog

Why is my reports service not connecting?

I am trying to pull some data from a Postgres database using Node.js and node-postures but I can't figure out why my service isn't connecting. my routes/index.js file: const express = require('express'); const router = express.Router(); const ordersCountController = require('../controllers/ordersCountController'); const ordersController = require('../controllers/ordersController'); const weeklyReportsController = require('../controllers/weeklyReportsController'); router.get('/orders_count', ordersCountController); router.get('/orders', ordersController); router.get('/weekly_reports', weeklyReportsController); module.exports = router; My controllers/weeklyReportsController.js file: const weeklyReportsService = require('../services/weeklyReportsService'); const weeklyReportsController = async (req, res) => { try { const data = await weeklyReportsService; res.json({data}) console...

How to show number of registered users in Laravel based on usertype?

i'm trying to display data from the database in the admin dashboard i used this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count(); echo $users; ?> and i have successfully get the correct data from the database but what if i want to display a specific data for example in this user table there is "usertype" that specify if the user is normal user or admin i want to user the same code above but to display a specific usertype i tried this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count()->WHERE usertype =admin; echo $users; ?> but it didn't work, what am i doing wrong? source https://stackoverflow.com/questions/68199726/how-to-show-number-of-registered-users-in-laravel-based-on-usertype

The number of nodes and edges I calculated manually and by a method provided by networkx library dont match

I am trying to analyse a graph using networkx library in Python. I have implemented a function that returns the nodes and edges of that graph. After that I counted the number of nodes and edges by implementing the len() function for nodes and edges, which gave me 17 for both nodes and edges. However, when I used the method graph.number_of_nodes() of networkx library it gave me 53. My code: import pandas as pd import networkx as nx import matplotlib.pyplot as plt header_list = ["a","b","w"] E =pd.read_csv("aves-barn-swallow-contact-network.edges",sep=" ",header= None,names=header_list) # print(E.head()) # prints first few rows of the csv file G = nx.from_pandas_edgelist(E,"a","b",["w"]) nx.draw(G) plt.show() class GraphAnalyse: #finding nodes and edges def nodeEdge(self,Graph): nodes = Graph.nodes() edges = Graph.edges() ...