James Tursa

Last seen: Today

3 978 total contributions since 2007

Interested in external applications, mex & engine applications with Fortran, C, C++, etc, particularly in the areas of speed improvement and memory management efficiencies.

Contact

James Tursa's Badges

View details...

Answered
Why can't complex integer arithmetic be done?
I think the main reason for disallowing integer array arithmetic in MATLAB is because of the ambiguity caused by the max(min( ))...

ungefär 8 timmar ago | 0

Answered
How can I split C/C++ Code in multiple files and still use MEX
Use the -c option to compile the B and C stuff to produce only object files without any linking: mex -c B.c C.c Then compile y...

3 dagar ago | 1

| accepted

Answered
Matrix multiplication. Trying to get scalar solution
Those should be matrix multiples, not element-wise multiplies. Remove the dots. E.g., top=nMode * M * u0; bot=nMode * M * nM...

5 dagar ago | 0

| accepted

Answered
Adapt RK4 ODE Solver from first order to 2nd order
Try this for starters yt = zeros(numel(y0),length(h:tfinal)); %Memory Allocation And then maybe this is what is intended for i...

7 dagar ago | 0

Answered
Error using vertcat Dimensions of arrays being concatenated are not consistent.
You have a mismatch in the number of arguments. E.g., this [V,y] = ode45(@(~,y) fun(del_h1,del_cp1,cp1),Vspan,y0); should be ...

8 dagar ago | 0

Answered
Replacing matrix n column values
E.g., temp = A(:,70:80); temp(temp<0) = 0; A(:,70:80) = temp;

8 dagar ago | 1

| accepted

Answered
Solving an overdetermined system of linear equations
Since you can simply plug in those values to verify that they satisfy all four equations, the conclusion would be that you have ...

8 dagar ago | 0

| accepted

Answered
How to create a cell array with binary numbers?
Your pop variable has 1000x11 = 11000 elements. Your for-loop only assigns values to 1000 of those elements (the first column)....

8 dagar ago | 0

Answered
Vector sum of forces
sqrt(sum(F.*F,2))

9 dagar ago | 0

Answered
invalid use of operator
You are missing a ( after the second "if", and you are missing an "end" statement. That being said, my suggestion would be to u...

11 dagar ago | 0

Answered
Need Help Solving Two Coupled ODEs
The derivative function needs to have the input arguments specified. In your current code above, you have no inputs specified. ...

15 dagar ago | 0

| accepted

Answered
Versions of Euler Methods
methodA has a fundamental flaw: y(i+1) = y(i)+ h*func( t(i) + h/2, (y(i)+y(i+1))/2 ); You can't use y(i+1) on the right hand s...

15 dagar ago | 1

Answered
Problem with matrix differential equation resolve in Matlab with ode45
Normally I would have expected ode45( ) to pass q as a 4x1 column vector. Also, I don't see anything in your code solving for q...

18 dagar ago | 0

| accepted

Answered
Question regarding quaternion conventions; particularly with respect to point vs frame rotations
See this post for comments regarding "point" vs "frame" convention: https://www.mathworks.com/matlabcentral/answers/524471-what...

18 dagar ago | 0

Answered
Using vectorization to shorten run time
E.g., [P,D] = meshgrid(p,d); % generate all possible combinations X = [P(:) D(:)]; % combine them into side-by-side column vec...

22 dagar ago | 1

| accepted

Answered
why does dimension of input variable change in ode45
The dimension changes to a column vector because that is the way ode45( ) operates. Your derivative function needs to be able t...

22 dagar ago | 1

| accepted

Answered
Maatlab crash with .mex function
I had to make some assumptions about how the matrix data areas are supposed to be loaded. We can try this first to see if these...

22 dagar ago | 0

Answered
How to calculate standard deviation in double for loop?
The short answer is don't do that. This style of programming has many problems. See this link for alternatives: https://www.m...

23 dagar ago | 0

Answered
Add row to an empty array if a condition is fulfilled
No need for an explicit loop. Just use logical indexing: result = m(strcmp(m(:,3),'A'),:);

23 dagar ago | 0

| accepted

Answered
reproduce random numbers for the third dimension
If you don't want the random numbers to change from iteration to iteration, then don't regenerate them inside the loop. Calcula...

23 dagar ago | 0

| accepted

Answered
Using an Array as an input
The error message actually tells you exactly what to do. Use element-wise operator .^ instead of the matrix operator ^ Energy ...

23 dagar ago | 1

Answered
Runge-Kutta 2
The + 4.*x1 should be - 4.*x1 in your derivative function. Also, using both X1 and x1 and X2 and x2, even though MATLAB is case...

23 dagar ago | 0

| accepted

Answered
3D array of rotation matrices multiplied on 2D array of vectors
To answer your side question of how to transpose 2D pages of an array, you can use the permute function to swap the first two di...

23 dagar ago | 0

Answered
Passing structs/objects to functions
MATLAB typically passes shared data copies of input arguments to functions. That means creating a separate mxArray header for t...

25 dagar ago | 0

Answered
How can I solve a system of ODE's, in which the equations have a dependency on time, as well as the other ODEs?
Just code it up as it appears using a 3-element state vector, where the elements are defined as follows y(1) = x y(2) = y y(3...

ungefär en månad ago | 0

Answered
How to improve the speed of large matrix addition?
If you are generating a new matrix C each time, then you are essentially doing two things: (1) Copying all the elements of A in...

ungefär en månad ago | 0

Answered
Dealing with very small or very large numbers
I guess the obvious suggestion is to work in units of nanometers or maybe picometers instead of meters. Can you do that?

ungefär en månad ago | 0

| accepted

Answered
How to use if statement in function?
Add this code to the top of your function: if( N <= 0 ) error('N needs to be > 0') end That being said, your loop probab...

ungefär en månad ago | 1

| accepted

Answered
How to build a quaternion from a normal vector?
If you have the Aerospace Toolbox, you could do the following n = 1x3 unit normal vector a = angle to rotate in radians v = 1...

ungefär en månad ago | 0

Answered
Runge-Kutta 2
Normally one would plot the position, not the velocity. So I would have expected you to plot the "1" index of your solutions. ...

ungefär en månad ago | 0