| 24 | #include <elementAPI.h> |
| 25 | |
| 26 | int OPS_sdfResponse() |
| 27 | { |
| 28 | int numData = 1; |
| 29 | double uresidual = 0.0; |
| 30 | double umaxprev = 0.0; |
| 31 | |
| 32 | int numOptionalArgs = 0; |
| 33 | int numArgs = OPS_GetNumRemainingInputArgs(); |
| 34 | double t_end = 0.0; bool t_end_specified = false; |
| 35 | while (OPS_GetNumRemainingInputArgs() > 0) { |
| 36 | std::string type = OPS_GetString(); |
| 37 | if (type == "-tend") { |
| 38 | numOptionalArgs++; |
| 39 | if (OPS_GetNumRemainingInputArgs() > 0) { |
| 40 | numData = 1; |
| 41 | if (OPS_GetDoubleInput(&numData, &t_end) < 0) { |
| 42 | opserr << "ERROR sdfResponse - failed to read tend" << endln; |
| 43 | return 0; |
| 44 | } |
| 45 | numOptionalArgs++; |
| 46 | t_end_specified = true; |
| 47 | } |
| 48 | } |
| 49 | if (type == "-uresid" || type == "-uresidual") { |
| 50 | numOptionalArgs++; |
| 51 | if (OPS_GetNumRemainingInputArgs() > 0) { |
| 52 | numData = 1; |
| 53 | if (OPS_GetDoubleInput(&numData, &uresidual) < 0) { |
| 54 | opserr << "ERROR sdfResponse - failed to read uresidual" << endln; |
| 55 | return 0; |
| 56 | } |
| 57 | numOptionalArgs++; |
| 58 | } |
| 59 | } |
| 60 | if (type == "-umaxprev") { |
| 61 | numOptionalArgs++; |
| 62 | if (OPS_GetNumRemainingInputArgs() > 0) { |
| 63 | numData = 1; |
| 64 | if (OPS_GetDoubleInput(&numData, &umaxprev) < 0) { |
| 65 | opserr << "ERROR sdfResponse - failed to read umaxprev" << endln; |
| 66 | return 0; |
| 67 | } |
| 68 | numOptionalArgs++; |
| 69 | } |
| 70 | } |
| 71 | } |
| 72 | |
| 73 | if (numArgs > 0) { |
| 74 | OPS_ResetCurrentInputArg(-numArgs); |
| 75 | } |
| 76 | numArgs = numArgs - numOptionalArgs; |
| 77 | |
| 78 | if (numArgs < 7 || numArgs > 8) { |
| 79 | opserr << "Incorrect number of arguments to sdfResponse --"; |
| 80 | opserr << "m, zeta, k, Fy, alpha, dtF, filename, dt, <-uresidual, uresid, -umaxprev, umaxp, -tend, tend>" << endln; |
| 81 | opserr << "m, zeta, k, Fy, alpha, tsTag, dt, <-uresidual, uresid, -umaxprev, umaxp, -tend, tend>" << endln; |
| 82 | return -1; |
| 83 | } |
no test coverage detected