Scientific Computation MCP

Installation Guide

Claude Desktop

Open Claude Desktop's configuration file (claude_desktop_config.json) and add the following:

• Mac/Linux:

{
  "mcpServers": {
    "numpy_mcp": {
      "command": "npx",
      "args": [
        "-y",
        "@smithery/cli@latest",
        "run",
        "@Aman-Amith-Shastry/scientific_computation_mcp",
        "--key",
        "<YOUR_SMITHERY_API_KEY>"
      ]
    }
  }
}

• Windows:

{
  "mcpServers": {
    "numpy_mcp": {
      "command": "cmd",
      "args": [
        "/c",
        "npx",
        "-y",
        "@smithery/cli@latest",
        "run",
        "@Aman-Amith-Shastry/scientific_computation_mcp",
        "--key",
        "<YOUR_SMITHERY_API_KEY>"
      ]
    }
  }
}

Or alternatively, run the following command:

npx -y @smithery/cli@latest install @Aman-Amith-Shastry/scientific_computation_mcp --client claude --key <YOUR_SMITHERY_API_KEY>

Restart Claude to load the server properly

Cursor

If you prefer to access the server through Cursor instead, then run the following command:

npx -y @smithery/cli@latest install @Aman-Amith-Shastry/scientific_computation_mcp --client cursor --key <YOUR_SMITHERY_API_KEY>

Related MCP server: Wolfram Alpha MCP Server

Components of the Server

Tools

Tensor storage

• create_tensor: Creates a new tensor based on a given name, shape, and values, and adds it to the tensor store. For the purposes of this server, tensors are vectors and matrices.

• view_tensor: Display the contents of a tensor from the store .

• delete_tensor: Deletes a tensor based on its name in the tensor store.

Linear Algebra

• add_matrices: Adds two matrices with the provided names, if compatible.

• subtract_matrices: Subtracts two matrices with the provided names, if compatible.

• multiply_matrices: Multiplies two matrices with the provided names, if compatible.

• scale_matrix: Scales a matrix of the provided name by a certain factor, in-place by default.

• matrix_inverse: Computes the inverse of the matrix with the provided name.

• transpose: Computes the transpose of the inverse of the matrix of the provided name.

• determinant: Computes the determinant of the matrix of the provided name.

• rank: Computes the rank (number of pivots) of the matrix of the provided name.

• compute_eigen: Calculates the eigenvectors and eigenvalues of the matrix of the provided name.

• qr_decompose: Computes the QR factorization of the matrix of the provided name. The columns of Q are an orthonormal basis for the image of the matrix, and R is upper triangular.

• svd_decompose: Computes the Singular Value Decomposition of the matrix of the provided name.

• find_orthonormal_basis: Finds an orthonormal basis for the matrix of the provided name. The vectors returned are all pair-wise orthogonal and are of unit length.

• change_basis: Computes the matrix of the provided name in the new basis.

Vector Calculus

• vector_project: Projects a vector in the tensor store to the specified vector in the same vector space

• vector_dot_product: Computes the dot product of two vectors in the tensor stores based on their provided names.

• vector_cross_product: Computes the cross product of two vectors in the tensor stores based on their provided names.

• gradient: Computes the gradient of a multivariable function based on the input function. Example call: gradient("x^2 + 2xyz + zy^3"). Do NOT include the function name (like f(x, y, z) = ...`).

• curl: Computes the curl of a vector field based on the input vector field. The input string must be formatted as a python list. Example call: curl("[3xy, 2z^4, 2y]"").

• divergenceComputes the divergence of a vector field based on the input vector field. The input string must be formatted as a python list. Example call: divergence("[3xy, 2z^4, 2y]"").

• laplacianComputes the laplacian of a scalar function (as the divergence of the gradient) or a vector field (where a component-wise laplacian is computed). If a scalar function is the input, it must be input in the same format as in the gradient tool. If the input is a vector field, it must be input in the same manner as the curl/divergence tools.

• directional_deriv: Computes the directional derivative of a function in a given direction u By default, the tool normalizes u before computing the directional derivative, as specified by the unit parameter.

Visualization

• plot_vector_field: Plots a vector field (specified in the same format as in the curl/divergence functions). Currently, only 3d vector fields are supported. A 2d png perspective image of the vector field is returned. By default, the bounds of the graph are from -1 to 1 on each axis.

• plot_function: Plots a function in 2d or 3d (based on the input variables), specified in the same format as in the gradient tool. Only the variables x and y can be used.